Genetic Environment of cry1 Genes Indicates Their Common Origin

Mosquitocidal toxin-like islands in Bacillus thuringiensis S2160-1 revealed by complete-genome sequence and MS proteomic analysis

Here, we present the whole genome sequence of Bt S2160-1, a potential alternative to the mosquitocidal model strain, Bti. One chromosome genome and four mega-plasmids were contained in Bt S2160-1, and 13 predicted genes encoding predicted insecticidal crystal proteins were identified clustered on one plasmid pS2160-1p2 containing two pathogenic islands (PAIs) designed as PAI-1 (Cry54Ba, Cry30Ea4, Cry69Aa-like, Cry50Ba2-like, Cry4Ca1-like, Cry30Ga2, Cry71Aa-like, Cry72Aa-like, Cry70Aa-like, Cyt1Da2-like and Vpb4C1-like) and PAI-2 (Cyt1Aa-like, and Tpp80Aa1-like). The clusters appear to represent mosquitocidal toxin islands similar to pathogenicity islands. Transcription/translation of 10 of the 13 predicted genes was confirmed by whole-proteome analysis using LTQ-Orbitrap LC-MS/MS. In summary, the present study identified the existence of a mosquitocidal toxin island in Bacillus thuringiensis, and provides important genomic information for understanding the insecticidal mechanism of B. thuringiensis.

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.

Identification of cry genes in Bacillus thuringiensis by multiplex real-time PCR

Bacillus thuringiensis is an important bacterium of the group Bacillus cereus sensu lato due to its insecticidal properties. This microorganism has high genetic variability and its strains produce different Cry toxins, known as δ-endotoxins, which are mainly responsible for its toxic effect on insects that are agricultural pests or vector human diseases. Each strain can express a variety of cry genes, out of a total of 789 cry genes described so far. The detection of these genes is very important to characterize strains, as they may indicate their toxic potential. Several methods have been used to characterize B. thuringiensis strains, but one of the most common techniques is Polymerase Chain Reaction (PCR) from primers that detect the presence of cry genes. This technique has been optimized to make real-time multiplex quantitative PCR (qPCR) assays faster, more efficient, and safer, because the presence of three genes can be detected in a single reaction. In this work, a multiplex assay was developed to identify the presence of genes from the cry1A, cry1C, and cry1F families whose respective toxins are present in both bioinsecticides, and commercial transgenic plants used to control caterpillars. Specific primers were designed to identify the families of the cited genes and the system was validated with samples that were sequenced by next-generation sequencing (NGS). The system was implemented and used to characterize 214 strains. Of these, eight were submitted to conventional PCR, and the results matched, again validating the system. Thus, the application of the proposed technique allows the reliable evaluation through this system to detect the presence of the genes of the families cry1A, cry1C, and cry1F in samples of B. thuringiensis.

Elucidating the genomic history of commercially used Bacillus thuringiensis subsp. tenebrionis strain NB176

Bacillus thuringiensis subsp. tenebrionis (Btt) produces a coleopteran-specific crystal protoxin protein (Cry3Aa δ-endotoxin). After its discovery in 1982, the strain NB125 (DSM 5526) was eventually registered in 1990 to control the Colorado potato beetle (Leptinotarsa decemlineata). Gamma-irradiation of NB125 resulted in strain NB176-1 (DSM 5480) that exhibited higher cry3Aa production and became the active ingredient of the plant protection product Novodor^® FC. Here, we report a comparative genome analysis of the parental strain NB125, its derivative NB176-1 and the current commercial production strain NB176. The entire genome sequences of the parental and derivative strains were deciphered by a hybrid de novo approach using short (Illumina) and long (Nanopore) read sequencing techniques. Genome assembly revealed a chromosome of 5.4 to 5.6 Mbp and six plasmids with a size range from 14.9 to 250.5 kbp for each strain. The major differences among the original NB125 and the derivative strains NB176-1 and NB176 were an additional copy of the cry3Aa gene, which translocated to another plasmid as well as a chromosomal deletion (~ 178 kbp) in NB176. The assembled genome sequences were further analyzed in silico for the presence of virulence and antimicrobial resistance (AMR) genes.

Conjugation across Bacillus cereus and kin: A review

Horizontal gene transfer (HGT) is a major driving force in shaping bacterial communities. Key elements responsible for HGT are conjugation-like events and transmissible plasmids. Conjugative plasmids can promote their own transfer as well as that of co-resident plasmids. Bacillus cereus and relatives harbor a plethora of plasmids, including conjugative plasmids, which are at the heart of the group species differentiation and specification. Since the first report of a conjugation-like event between strains of B. cereus sensu lato (s.l.) 40 years ago, many have studied the potential of plasmid transfer across the group, especially for plasmids encoding major toxins. Over the years, more than 20 plasmids from B. cereus isolates have been reported as conjugative. However, with the increasing number of genomic data available, in silico analyses indicate that more plasmids from B. cereus s.l. genomes present self-transfer potential. B. cereus s.l. bacteria occupy diverse environmental niches, which were mimicked in laboratory conditions to study conjugation-related mechanisms. Laboratory mating conditions remain nonetheless simplistic compared to the complex interactions occurring in natural environments. Given the health, economic and ecological importance of strains of B. cereus s.l., it is of prime importance to consider the impact of conjugation within this bacterial group.

Whole Genome Sequencing Analysis of Bacillus thuringiensis GR007 Reveals Multiple Pesticidal Protein Genes

Bacillus thuringiensis (Bt) are soil ubiquitous bacteria. They produce a great variability of insecticidal proteins, where certain of these toxins are used worldwide for pest control. Through their adaptation to diverse ecosystems, certain Bt strains have acquired genetic mobile elements by horizontal transfer, harboring genes that encode for different virulent factors and pesticidal proteins (PP). Genomic characterization of Bt strains provides a valuable source of PP with potential biotechnological applications for pest control. In this work, we have sequenced the complete genome of the bacterium Bt GR007 strain that is toxic to Spodoptera frugiperda and Manduca sexta larvae. Four replicons (one circular chromosome and three megaplasmids) were identified. The two largest megaplasmids (pGR340 and pGR157) contain multiple genes that codify for pesticidal proteins: 10 cry genes (cry1Ab, cry1Bb, cry1Da, cry1Fb, cry1Hb, cry1Id, cry1Ja, cry1Ka, cry1Nb, and cry2Ad), two vip genes (vip3Af and vip3Ag), two binary toxin genes (vpa2Ac and vpb1Ca), five genes that codify for insecticidal toxin components (Tc’s), and a truncated cry1Bd-like gene. In addition, genes that codify for several virulent factors were also found in this strain. Proteomic analysis of the parasporal crystals of GR007 revealed that they are composed of eight Cry proteins. Further cloning of these genes for their individual expression in Bt acrystalliferous strain, by means of their own intrinsic promoter showed expression of seven Cry proteins. These proteins display differential toxicity against M. sexta and S. frugiperda larvae, where Cry1Bb showed to be the most active protein against S. frugiperda larvae and Cry1Ka the most active protein against M. sexta larvae.

Pan-Genome Portrait of Bacillus mycoides Provides Insights into the Species Ecology and Evolution

Bacillus mycoides is poorly known despite its frequent occurrence in a wide variety of environments. To provide direct insight into its ecology and evolutionary history, a comparative investigation of the species pan-genome and the functional gene categorization of 35 isolates obtained from soil samples from northeastern Poland was performed. The pan-genome of these isolates is composed of 20,175 genes and is characterized by a strong predominance of adaptive genes (∼83%), a significant amount of plasmid genes (∼37%), and a great contribution of prophages and insertion sequences. The pan-genome structure and phylodynamic studies had suggested a wide genomic diversity among the isolates, but no correlation between lineages and the bacillus origin was found. Nevertheless, the two B. mycoides populations, one from Białowieża National Park, the last European natural primeval forest with soil classified as organic, and the second from mineral soil samples taken in a farm in Jasienówka, a place with strong anthropogenic pressure, differ significantly in the frequency of genes encoding proteins enabling bacillus adaptation to specific stress conditions and production of a set of compounds, thus facilitating their colonization of various ecological niches. Furthermore, differences in the prevalence of essential stress sigma factors might be an important trail of this process. Due to these numerous adaptive genes, B. mycoides is able to quickly adapt to changing environmental conditions.

IMPORTANCE This research allows deeper understanding of the genetic organization of natural bacterial populations, specifically, Bacillus mycoides, a psychrotrophic member of the Bacillus cereus group that is widely distributed worldwide, especially in areas with continental cold climates. These thorough analyses made it possible to describe, for the first time, the B. mycoides pan-genome, phylogenetic relationship within this species, and the mechanisms behind the species ecology and evolutionary history. Our study indicates a set of functional properties and adaptive genes, in particular, those encoding sigma factors, associated with B. mycoides acclimatization to specific ecological niches and changing environmental conditions.

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.

A novel anti-dipteran Bacillus thuringiensis strain: Unusual Cry toxin genes in a highly dynamic plasmid environment

Bacillus thuringiensis emerged as a major bioinsecticide on the global market. It offers a valuable alternative to chemical products classically utilized to control pest insects. Despite the efficiency of several strains and products available on the market, the scientific community is always on the lookout for novel toxins that can replace or supplement the existing products. In this study, H3, a novel B. thuringiensis strain showing mosquitocidal activity, was isolated from Lebanese soil and characterized at an in vivo, genomic and proteomic levels. H3 parasporal crystal is toxic on its own but displays an unusual killing profile with a higher LC50 than the reference B. thuringiensis serovar israelensis crystal proteins. In addition, H3 has a different toxicity order: it is more toxic to Aedes albopictus and Anopheles gambiae than to Culex pipiens Whole genome sequencing and crystal analysis revealed that H3 can produce eleven novel Cry proteins, eight of which are assembled in genes with an orf1-gap-orf2 organization, where orf2 is a potential Cry4-type crystallization domain. Moreover, pH3-180, the toxin-carrying plasmid, holds a wide repertoire of mobile genetic elements that amount to ca 22% of its size., including novel insertion sequences and class II transposable elements Two other large plasmids present in H3 carry genetic determinants for the production of many interesting molecules - such as chitinase, cellulase and bacitracin - that may add up to H3 bioactive properties. This study therefore reports a novel mosquitocidal Bacillus thuringiensis strain with unusual Cry toxin genes in a rich mobile DNA environment.IMPORTANCE Bacillus thuringiensis, a soil entomopathogenic bacteria, is at the base of many sustainable eco-friendly bio-insecticides. Hence stems the need to continually characterize insecticidal toxins. H3 is an anti-dipteran B. thuringiensis strain, isolated from Lebanese soil, whose parasporal crystal contains eleven novel Cry toxins and no Cyt toxins. In addition to its individual activity, H3 showed potential as a co-formulant with classic commercialized B. thuringiensis products, to delay the emergence of resistance and to shorten the time required for killing. On a genomic level, H3 holds three large plasmids, one of which carries the toxin-coding genes, with four occurrences of the distinct orf1-gap-orf2 organization. Moreover, this plasmid is extremely rich in mobile genetic elements, unlike its two co-residents. This highlights the important underlying evolutionary traits between toxin-carrying plasmids and the adaptation of a B. thuringiensis strain to its environment and insect host spectrum.

Application of Bacillus thuringiensis strains with conjugal and mobilizing capability drives gene transmissibility within Bacillus cereus group populations in confined habitats

Background

Bacillus thuringiensis bacteria share similar genetic, physiological, and biochemical characteristics with other members of the Bacillus cereus group. Their diversity and entomopathogenic origin are related to their mobile genetic elements. However, the effects of wide-spread application of B. thuringiensis-based pesticides on genetically related B. cereus group populations present in the environment remain poorly understood.

Results

We first identified pBMB76 from B. thuringiensis tenebrionis as a new conjugative plasmid. Mixed mating experiments suggested that pBMB76 may compete with pHT73, another known conjugative plasmid. Applications of single (tenebrionis 4AA1 and kurstaki HD73 carrying pBMB76 and pHT73, respectively) and mixed (4AA1 + HD73) B. thuringiensis strains were performed in confined plot habitats (soil and leaf) over two planting seasons. In total, 684 B. cereus group isolates were randomly selected from different treatment sets, and the transmissibility and occurrence rate of potential conjugative plasmids were surveyed. Results showed that the percentage of isolates with plasmid mobility was markedly enhanced in the B. thuringiensis-sprayed groups. Furthermore, we performed multi-locus sequence typing (MLST) for a subset of 291 isolates, which indicated that the dominant sequence types in the treated habitats were identical or related to the corresponding sprayed formulations.

Conclusions

The application of B. thuringiensis strains with conjugal and mobilizing capability drove gene transmissibility within the B. cereus group populations in confined habitats and potentially modified the population structure.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-020-02047-4.

Completed Genomic Sequence of Bacillus thuringiensis HER1410 Reveals a Cry-Containing Chromosome, Two Megaplasmids, and an Integrative Plasmidial Prophage

Bacillus thuringiensis is the most used biopesticide in agriculture. Its entomopathogenic capacity stems from the possession of plasmid-borne insecticidal crystal genes (cry), traditionally used as discriminant taxonomic feature for that species. As such, crystal and plasmid identification are key to the characterization of this species. To date, about 600 B. thuringiensis genomes have been reported, but less than 5% have been completed, while the other draft genomes are incomplete, hindering full plasmid delineation. Here we present the complete genome of Bacillus thuringiensis HER1410, a strain closely related to B. thuringiensis entomocidus and a known host for a variety of Bacillus phages. The combination of short and long-read techniques allowed fully resolving the genome and delineation of three plasmids. This enabled the accurate detection of an unusual location of a unique cry gene, cry1Ba4, located in a genomic island near the chromosome replication origin. Two megaplasmids, pLUSID1 and pLUSID2 could be delineated: pLUSID1 (368 kb), a likely conjugative plasmid involved in virulence, and pLUSID2 (156 kb) potentially related to the sporulation process. A smaller plasmidial prophage pLUSID3, with a dual lifestyle whose integration within the chromosome causes the disruption of a flagellar key component. Finally, phylogenetic analysis placed this strain within a clade comprising members from the B. thuringiensis serovar thuringiensis and other serovars and with B. cereus s. s. in agreement with the intermingled taxonomy of B. cereus sensu lato group.

Diversity of Bacillus cereus sensu lato mobilome

BACKGROUND

Bacillus cereus sensu lato s.l.) is a group of bacteria displaying close phylogenetic relationships but a high ecological diversity. The three most studied species are Bacillus anthracis, Bacillus cereus sensu stricto and Bacillus thuringiensis. While some species are pathogenic to mammals or associated with food poisoning, Bacillus thuringiensis is a well-known entomopathogenic bacterium used as biopesticide worldwide. B. cereus s.l. also contains a large variety of mobile genetic elements (MGEs).

RESULTS

In this study, we detail the occurrence and plasmid vs. chromosome distribution of several MGEs in 102 complete and annotated genomes of B. cereus s.l. These MGEs include 16 Insertion Sequence (IS) families, the Tn3 family, 18 different Bacillus cereus repeats (BCRs) and 30 known group II introns.

CONCLUSIONS

Our analysis not only shows the diversity of these MGEs among strains of the same species and between different species within the B. cereus s.l. group, but also highlights the potential impact of these elements on the plasticity of the plasmid pool, and the TEs (Transposable Elements) - species relationship within B. cereus s.l.

Role of plasmid plasticity and mobile genetic elements in the entomopathogen Bacillus thuringiensis serovar israelensis

Bacillus thuringiensis is a well-known biopesticide that has been used for more than 80 years. This spore-forming bacterium belongs to the group of Bacillus cereus that also includes, among others, emetic and diarrheic pathotypes of B. cereus, the animal pathogen Bacillus anthracis and the psychrotolerant Bacillus weihenstephanensis. Bacillus thuringiensis is rather unique since it has adapted its lifestyle as an efficient pathogen of specific insect larvae. One of the peculiarities of B. thuringiensis strains is the extent of their extrachromosomal pool, with strains harbouring more than 10 distinct plasmid molecules. Among the numerous serovars of B. thuringiensis, 'israelensis' is certainly emblematic since its host spectrum is apparently restricted to dipteran insects like mosquitoes and black flies, vectors of human and animal diseases such as malaria, yellow fever, or river blindness. In this review, the putative role of the mobile gene pool of B. thuringiensis serovar israelensis in its pathogenicity and dedicated lifestyle is reviewed, with specific emphasis on the nature, diversity, and potential mobility of its constituents. Variations among the few related strains of B. thuringiensis serovar israelensis will also be reported and discussed in the scope of this specialised insect pathogen, whose lifestyle in the environment remains largely unknown.

Bacterial Sexuality at the Nanoscale

Understanding the basic mechanisms of bacterial sexuality is an important topic in current microbiology and biotechnology. While classical methods used to study gene transfer provide information on whole cell populations, nanotechnologies offer new opportunities for analyzing the behavior of individual mating partners. We introduce an innovative atomic force microscopy (AFM) platform to study and mechanically control DNA transfer between single bacteria, focusing on the large conjugative pXO16 plasmid of the Gram-positive bacterium Bacillus thuringiensis. We demonstrate that the adhesion forces between single donor and recipient cells are very strong (∼2 nN). Using a mutant plasmid, we find that these high forces are mediated by a pXO16 aggregation locus that contains two large surface protein genes. Notably, we also show that AFM can be used to mechanically induce plasmid transfer between single partners, revealing that transfer is very fast (<15 min) and triggers major cell surface changes in transconjugant cells. We anticipate that the single-cell technology developed here will enable researchers to mechanically control gene transfer among a wide range of Gram-positive and Gram-negative bacterial species and to understand the molecular forces involved. Also, the method could be useful in nanomedicine for the design of antiadhesion compounds capable of preventing intimate cell-cell contacts, therefore providing a means to control the resistance and virulence of bacterial pathogens.