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

Interspecies Horizontal Transfer and Specific Integration of the Mosquitocidal Toxin-Encoding Plasmid pTAND672-2 from Bacillus thuringiensis subsp. israelensis to Lysinibacillus sphaericus

pTAND672-2, a 144-kb resident plasmid of Bacillus thuringiensis serovar israelensis strain TAND672, was sequenced and characterized. This extrachromosomal element carries mosquitocidal toxin-, conjugation-, and recombinase-encoding genes, together with a putative arbitrium system, a genetic module recently discovered in temperate phages controlling lysogeny-lysis transition and in mobile genetic elements (MGEs) where its function remains clarified. Using conjugation experiments, pTAND672-2 is shown to be a novel integrative and conjugative element (ICE), which can horizontally transfer from B. thuringiensis serovar israelensis to Lysinibacillus sphaericus, another mosquitocidal bacterium, where it integrates into the chromosome. Its integration and circularization are reversible and involve a single-cross recombination between 33-bp specific sites, attB in the chromosome of L. sphaericus and attP in pTAND672-2. CDS143, coding for the putative tyrosine integrase Int143 distantly related to site-specific tyrosine Xer recombinases and phage integrases, can mediate the integration of pTAND672-2 to attB. The B. thuringiensis mosquito-killing genes carried by pTAND672-2 are efficiently transcribed and expressed in L. sphaericus, displaying a slight increased toxicity in this bacterium against Aedes albopictus larvae. The occurrence of pTAND672-2-like plasmids within the Bacillus cereus group was also explored and indicated that they all share a similar genetic backbone with diverse plasmid sizes, ranging from 58 to 225 kb. Interestingly, among them, the pEFR-4-4 plasmid of Bacillus paranthracis EFR-4 and p5 of B. thuringiensis BT-59 also display conjugative capability; moreover, like pTAND672-2 displays a chimeric structure between the pCH_133-e- and pBtoxis-like plasmids, pBTHD789-3 also appears to be mosaic of two plasmids. IMPORTANCE Horizontal transfer of mobile genetic elements carrying mosquitocidal toxin genes may play a driving role in the diversity of mosquitocidal bacteria. Here, the 144-kb mosquitocidal toxin-encoding plasmid pTAND672-2 is the first verified integrative and conjugative element (ICE) identified in Bacillus thuringiensis serovar israelensis. The key tyrosine integrase Int143, involved in the specific integration, is distantly related to other tyrosine recombinases. The study also reports the occurrence and potential interspecies transmission of pTAND672-2-like plasmids with varied sizes in B. thuringiensis, Bacillus paranthracis, and Bacillus wiedmannii isolates belonging to the Bacillus cereus group. This study is important for further understanding the evolution and ecology of mosquitocidal bacteria, as well as for providing new direction for the genetic engineering of biopesticides in the control of disease-transmitting mosquitoes.

Comparative genomics of Bacillus cereus sensu lato spp. biocontrol strains in correlation to in-vitro phenotypes and plant pathogen antagonistic capacity

Bacillus cereus sensu lato (Bcsl) strains are widely explored due to their capacity to antagonize a broad range of plant pathogens. These include B. cereus sp. UW85, whose antagonistic capacity is attributed to the secondary metabolite Zwittermicin A (ZwA). We recently isolated four soil and root-associated Bcsl strains (MO2, S-10, S-25, LSTW-24) that displayed different growth profiles and in-vitro antagonistic effects against three soilborne plant pathogens models: Pythium aphanidermatum (oomycete) Rhizoctonia solani (basidiomycete), and Fusarium oxysporum (ascomycete). To identify genetic mechanisms potentially responsible for the differences in growth and antagonistic phenotypes of these Bcsl strains, we sequenced and compared their genomes, and that of strain UW85 using a hybrid sequencing pipeline. Despite similarities, specific Bcsl strains had unique secondary metabolite and chitinase-encoding genes that could potentially explain observed differences in in-vitro chitinolytic potential and anti-fungal activity. Strains UW85, S-10 and S-25 contained a (~500 Kbp) mega-plasmid that harbored the ZwA biosynthetic gene cluster. The UW85 mega-plasmid contained more ABC transporters than the other two strains, whereas the S-25 mega-plasmid carried a unique cluster containing cellulose and chitin degrading genes. Collectively, comparative genomics revealed several mechanisms that can potentially explain differences in in-vitro antagonism of Bcsl strains toward fungal plant pathogens.

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.

A small bacteriophage protein determines the hierarchy over co-residential jumbo phage in Bacillus thuringiensis serovar israelensis

Bacillus thuringiensis serovar israelensis is the most widely used biopesticide against insects, including vectors of animal and human diseases. Among several extrachromosomal elements, this endospore-forming entomopathogen harbors two bacteriophages: a linear DNA replicon named GIL01 that does not integrate into the chromosome during lysogeny and a circular-jumbo prophage known as pBtic235. Here, we show that GIL01 hinders the induction of cohabiting prophage pBtic235. The GIL01-encoded small protein, gp7, which interacts with the host LexA repressor, is a global transcription regulator and represses the induction of pBtic235 after DNA damage to presumably allow GIL01 to multiply first. In a complex with host LexA in stressed cells, gp7 down-regulates the expression of more than 250 host and pBtic235 genes, many of which are involved in the cellular functions of genome maintenance, cell-wall transport, and membrane and protein stability. We show that gp7 homologs that are found exclusively in bacteriophages act in a similar fashion to enhance LexA's binding to DNA, while likely also affecting host gene expression. Our results provide evidence that GIL01 influences both its host and its co-resident bacteriophage.

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.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

Conjugation-mediated transfer of pXO16, a large plasmid from Bacillus thuringiensis sv. israelensis, across the Bacillus cereus group and its impact on host phenotype

pXO16, the 350 kb-conjugative plasmid from Bacillus thuringiensis sv. israelensis promotes its own transfer at high efficiency, triggers the transfer of mobilizable and non-mobilizable plasmids, as well as the transfer of host chromosomal loci. Naturally found in B. thuringiensis sv. israelensis, pXO16 transfers to various strains of Bacillus cereus sensu lato (s.l.) at a wide range of frequencies. Despite this host diversity, a paradox remains between the relatively large host spectrum and the natural occurrence of pXO16, so far restricted to B. thuringiensis sv. israelensis. Proposing first insights exploring this paradox, we investigated the behaviour of pXO16 amongst different members of the B. cereus group. We first looked at the transfer of pXO16 to two new host clusters of B. cereus s.l., Bacillus mycoides and Bacillus anthracis clusters. This examination brought to light the impairment of the characteristic rhizoidal phenotype of B. mycoides in presence of pXO16. We also explored the stability of pXO16 at different temperatures as some B. cereus group members are well-known for their psychro- or thermo-tolerance. This shed light on the thermo-sensitivity of the plasmid. The influence of pXO16 on its host cell growth and on swimming capacity also revealed no or limited impact on its natural host B. thuringiensis sv. israelensis. On the contrary, pXO16 affected more strongly both the growth and swimming capacity of other B. cereus s.l. hosts. This reinforced the running hypothesis of a co-evolution between pXO16 and B. thuringiensis sv. israelensis, enabling the plasmid maintenance without impairing the host strain development.

Genomic insights into the diversity of non-coding RNAs in Bacillus cereus sensu lato

Bacillus cereus sensu lato is a group of bacteria of medical and agricultural importance in different ecological niches and with controversial taxonomic relationships. Studying the composition of non-coding RNAs (ncRNAs) in several bacterial groups has been an important tool for identifying genetic information and better understanding genetic regulation towards environment adaptation. However, to date, no comparative genomics study of ncRNA has been performed in this group. Thus, this study aimed to identify and characterize the set of ncRNAs from 132 strains of Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis to obtain an overview of the diversity and distribution of these genetic elements in these species. We observed that the number of ncRNAs differs in the chromosomes of the three species, but not in the plasmids, when species or phylogenetic clusters were compared. The prevailing functional/structural category was Cis-reg and the most frequent class was Riboswitch. However, in plasmids, the class Group II intron was the most frequent. Also, nine ncRNAs were selected for validation in the strain B. thuringiensis 407 by RT-PCR, which allowed to identify the expression of the ncRNAs. The wide distribution and diversity of ncRNAs in the B. cereus group, and more intensely in B. thuringiensis, may help improve the abilities of these species to adapt to various environmental changes. Further studies should address the expression of these genetic elements in different conditions.

Bacillus cytotoxicus Genomics: Chromosomal Diversity and Plasmidome Versatility

Bacillus cytotoxicus is the thermotolerant representative of the Bacillus cereus group. This group, also known as B. cereus sensu lato, comprises both beneficial and pathogenic members and includes psychrotolerant and thermotolerant species. Bacillus cytotoxicus was originally recovered from a fatal outbreak in France in 1998. This species forms a remote cluster from the B. cereus group members and reliably contains the cytk-1 gene, coding for a cytotoxic variant of cytotoxin K. Although this species was originally thought to be homogenous, intra-species diversity has been recently described with four clades, six random amplified polymorphic DNA (RAPD) patterns, and 11 plasmids profiles. This study aimed to get new insights into the genomic diversity of B. cytotoxicus and to decipher the underlying chromosomal and plasmidial variations among six representative isolates through whole genome sequencing (WGS). Among the six sequenced strains, four fitted the previously described genomic clades A and D, while the remaining two constituted new distinct branches. As for the plasmid content of these strains, three large plasmids were putatively conjugative and three small ones potentially mobilizable, harboring coding genes for putative leaderless bacteriocins. Mobile genetic elements, such as prophages, Insertion Sequences (IS), and Bacillus cereus repeats (bcr) greatly contributed to the B. cytotoxicus diversity. As for IS elements and bcr, IS3 and bcr1 were the most abundant elements and, along with the group II intron B.c.I8, were found in all analyzed B. cytotoxicus strains. When compared to other B. cytotoxicus strains, the type-strain NVH 391-98 displayed a relatively low number of IS. Our results shed new light on the contribution of mobile genetic elements to the genome plasticity of B. cytotoxicus and their potential role in horizontal gene transfer.

Unraveling Protein Interactions between the Temperate Virus Bam35 and Its Bacillus Host Using an Integrative Yeast Two Hybrid-High Throughput Sequencing Approach

Bacillus virus Bam35 is the model Betatectivirus and member of the family Tectiviridae, which is composed of tailless, icosahedral, and membrane-containing bacteriophages. Interest in these viruses has greatly increased in recent years as they are thought to be an evolutionary link between diverse groups of prokaryotic and eukaryotic viruses. Additionally, betatectiviruses infect bacteria of the Bacillus cereus group, which are known for their applications in industry and notorious since it contains many pathogens. Here, we present the first protein–protein interactions (PPIs) network for a tectivirus–host system by studying the Bam35–Bacillus thuringiensis model using a novel approach that integrates the traditional yeast two-hybrid system and high-throughput sequencing (Y2H-HTS). We generated and thoroughly analyzed a genomic library of Bam35′s host B. thuringiensis HER1410 and screened interactions with all the viral proteins using different combinations of bait–prey couples. Initial analysis of the raw data enabled the identification of over 4000 candidate interactions, which were sequentially filtered to produce 182 high-confidence interactions that were defined as part of the core virus–host interactome. Overall, host metabolism proteins and peptidases were particularly enriched within the detected interactions, distinguishing this host–phage system from the other reported host–phage PPIs. Our approach also suggested biological roles for several Bam35 proteins of unknown function, including the membrane structural protein P25, which may be a viral hub with a role in host membrane modification during viral particle morphogenesis. This work resulted in a better understanding of the Bam35–B. thuringiensis interaction at the molecular level and holds great potential for the generalization of the Y2H-HTS approach for other virus–host models.

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.

Unlimited Cooperativity of Betatectivirus SSB, a Novel DNA Binding Protein Related to an Atypical Group of SSBs From Protein-Primed Replicating Bacterial Viruses

Bam35 and related betatectiviruses are tail-less bacteriophages that prey on members of the Bacillus cereus group. These temperate viruses replicate their linear genome by a protein-primed mechanism. In this work, we have identified and characterized the product of the viral ORF2 as a single-stranded DNA binding protein (hereafter B35SSB). B35SSB binds ssDNA with great preference over dsDNA or RNA in a sequence-independent, highly cooperative manner that results in a non-specific stimulation of DNA replication. We have also identified several aromatic and basic residues, involved in base-stacking and electrostatic interactions, respectively, that are required for effective protein-ssDNA interaction. Although SSBs are essential for DNA replication in all domains of life as well as many viruses, they are very diverse proteins. However, most SSBs share a common structural domain, named OB-fold. Protein-primed viruses could constitute an exception, as no OB-fold DNA binding protein has been reported. Based on databases searches as well as phylogenetic and structural analyses, we showed that B35SSB belongs to a novel and independent group of SSBs. This group contains proteins encoded by protein-primed viral genomes from unrelated viruses, spanning betatectiviruses and Φ29 and close podoviruses, and they share a conserved pattern of secondary structure. Sensitive searches and structural predictions indicate that B35SSB contains a conserved domain resembling a divergent OB-fold, which would constitute the first occurrence of an OB-fold-like domain in a protein-primed genome.

ORFograph: search for novel insecticidal protein genes in genomic and metagenomic assembly graphs

BACKGROUND

Since the prolonged use of insecticidal proteins has led to toxin resistance, it is important to search for novel insecticidal protein genes (IPGs) that are effective in controlling resistant insect populations. IPGs are usually encoded in the genomes of entomopathogenic bacteria, especially in large plasmids in strains of the ubiquitous soil bacteria, Bacillus thuringiensis (Bt). Since there are often multiple similar IPGs encoded by such plasmids, their assemblies are typically fragmented and many IPGs are scattered through multiple contigs. As a result, existing gene prediction tools (that analyze individual contigs) typically predict partial rather than complete IPGs, making it difficult to conduct downstream IPG engineering efforts in agricultural genomics.

METHODS

Although it is difficult to assemble IPGs in a single contig, the structure of the genome assembly graph often provides clues on how to combine multiple contigs into segments encoding a single IPG.

RESULTS

We describe ORFograph, a pipeline for predicting IPGs in assembly graphs, benchmark it on (meta)genomic datasets, and discover nearly a hundred novel IPGs. This work shows that graph-aware gene prediction tools enable the discovery of greater diversity of IPGs from (meta)genomes.

CONCLUSIONS

We demonstrated that analysis of the assembly graphs reveals novel candidate IPGs. ORFograph identified both already known genes "hidden" in assembly graphs and potential novel IPGs that evaded existing tools for IPG identification. As ORFograph is fast, one could imagine a pipeline that processes many (meta)genomic assembly graphs to identify even more novel IPGs for phenotypic testing than would previously be inaccessible by traditional gene-finding methods. While here we demonstrated the results of ORFograph only for IPGs, the proposed approach can be generalized to any class of genes. Video abstract.

Comparative Genomics of Prophages Sato and Sole Expands the Genetic Diversity Found in the Genus Betatectivirus

Tectiviruses infecting the Bacillus cereus group represent part of the bacterial "plasmid repertoire" as they behave as linear plasmids during their lysogenic cycle. Several novel tectiviruses have been recently found infecting diverse strains belonging the B. cereus lineage. Here, we report and analyze the complete genome sequences of phages Sato and Sole. The linear dsDNA genome of Sato spans 14,852 bp with 32 coding DNA sequences (CDSs), whereas the one of Sole has 14,444 bp comprising 30 CDSs. Both phage genomes contain inverted terminal repeats and no tRNAs. Genomic comparisons and phylogenetic analyses placed these two phages within the genus Betatectivirus in the family Tectiviridae. Additional comparative genomic analyses indicated that the "gene regulation-genome replication" module of phages Sato and Sole is more diverse than previously observed among other fully sequenced betatectiviruses, displaying very low sequence similarities and containing some ORFans. Interestingly, the ssDNA binding protein encoded in this genomic module in phages Sato and Sole has very little amino acid similarity with those of reference betatectiviruses. Phylogenetic analyses showed that both Sato and Sole represent novel tectivirus species, thus we propose to include them as two novel species in the genus Betatectivirus.

Preparation of Nanomaterial Wettable Powder Formulations of Antagonistic Bacteria from Phellodendron chinense and the Biological Control of Brown Leaf Spot Disease

Brown leaf spot disease caused by Nigrospora guilinensis on Phellodendron chinense occurs in a large area in Dayi County, Chengdu City, Sichuan Province, China each year. This outbreak has severely reduced the production of Chinese medicinal plants P. chinense and caused substantial economic losses. The bacterial isolate JKB05 was isolated from the healthy leaves of P. chinense, exhibited antagonistic effects against N. guilinensis and was identified as Bacillus megaterium. The following fermentation medium and conditions improved the inhibitory effect of B. megaterium JKB05 on N. guilinensis: 2% glucose, 0.1% soybean powder, 0.1% KCl, and 0.05% MgSO4; initial concentration 6 × 106 cfu/ml, and a 42-h optimal fermentation time. A composite of 0.1% nano-SiO2 JKB05 improved the thermal stability, acid-base stability and ultraviolet resistance by 16%, 12%, and 38.9%, respectively, and nano-SiO2 was added to the fermentation process. The best formula for the wettable powder was 35% kaolin, 4% polyethylene glycol, 8% Tween, and 2% humic acid. The following quality test results for the wettable powder were obtained: wetting time 87.0 s, suspension rate 80.33%, frequency of microbial contamination 0.08%, pH 7.2, fineness 95.8%, drying loss 1.47%, and storage stability ≥83.5%. A pot experiment revealed that the ability of JKB05 to prevent fungal infections on P. chinense increased considerably and achieved levels of control as high as 94%. The use of nanomaterials significantly improved the ability of biocontrol bacteria to control this disease.

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.

Potential for Bacillus thuringiensis and Other Bacterial Toxins as Biological Control Agents to Combat Dipteran Pests of Medical and Agronomic Importance

The control of dipteran pests is highly relevant to humans due to their involvement in the transmission of serious diseases including malaria, dengue fever, Chikungunya, yellow fever, zika, and filariasis; as well as their agronomic impact on numerous crops. Many bacteria are able to produce proteins that are active against insect species. These bacteria include Bacillus thuringiensis, the most widely-studied pesticidal bacterium, which synthesizes proteins that accumulate in crystals with insecticidal properties and which has been widely used in the biological control of insects from different orders, including Lepidoptera, Coleoptera, and Diptera. In this review, we summarize all the bacterial proteins, from B. thuringiensis and other entomopathogenic bacteria, which have described insecticidal activity against dipteran pests, including species of medical and agronomic importance.

Long inverted repeats around the chromosome replication terminus in the model strain Bacillus thuringiensis serovar israelensis BGSC 4Q7

Bacillus thuringiensis serovar israelensis is the most widely used natural biopesticide against mosquito larvae worldwide. Its lineage has been actively studied and a plasmid-free strain, B. thuringiensis serovar israelensis BGSC 4Q7 (4Q7), has been produced. Previous sequencing of the genome of this strain has revealed the persistent presence of a 235 kb extrachromosomal element, pBtic235, which has been shown to be an inducible prophage, although three putative chromosomal prophages have been lost. Moreover, a 492 kb region, potentially including the standard replication terminus, has also been deleted in the 4Q7 strain, indicating an absence of essential genes in this area. We reanalysed the genome coverage distribution of reads for the previously sequenced variant strain, and sequenced two independently maintained samples of the 4Q7 strain. A 553 kb area, close to the 492 kb deletion, was found to be duplicated. This duplication presumably restored the equal sizes of the replichores, and a balanced functioning of replication termination. An analysis of genome assembly graphs revealed a transient association of the host chromosome with the pBtic235 element. This association may play a functional role in the replication of the bacterial chromosome, and the termination of this process in particular. The genome-restructuring events detected may modify the genetic status of cytotoxic or haemolytic toxins, potentially influencing strain virulence. Twelve of the single-nucleotide variants identified in 4Q7 were probably due to the procedure used for strain construction or were present in the precursor of this strain. No sequence variants were found in pBtic235, but the distribution of the corresponding 4Q7 reads indicates a significant difference from counterparts in natural B. thuringiensis serovar israelensis strains, suggesting a duplication or over-replication in 4Q7. Thus, the 4Q7 strain is not a pure plasmid-less offshoot, but a highly genetically modified derivative of its natural ancestor. In addition to potentially influencing virulence, genome-restructuring events can modify the replication termination machinery. These findings have potential implications for the conclusions of virulence studies on 4Q7 as a model, but they also raise interesting fundamental questions about the functioning of the Bacillus genome.

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.

Collagen-like sequences encoded by extremophilic and extremotolerant bacteria

Collagens and collagen-like proteins are found in a wide range of organisms. The common feature of these proteins is a triple helix fold, requiring a characteristic pattern of amino acid sequences, composed of Gly-X-Y tripeptide repeats. Collagen-like proteins from bacteria are heterogeneous in terms of length and amino acid composition of their collagenous sequences. However, different bacteria live in different environments, some at extreme temperatures and conditions. This study explores the occurrence of collagen-like sequences in the genomes of different extreme condition-adapted bacteria, and investigates features that could be linked to conditions where they thrive. Our results show that proteins containing collagen-like sequences are encoded by genomes of various extremophiles. Some of these proteins contain conservative domains, characteristic of cell or endospore surface proteins, while most other proteins are unknown. The characteristics of collagenous sequences may depend on both, the phylogenetic relationship and the living conditions of the bacteria.

Characterization and Whole Genome Sequencing of AR23, a Highly Toxic Bacillus thuringiensis Strain Isolated from Lebanese Soil

The demand for sustainable and eco-friendly control methods of pests and insects is increasing worldwide. From this came the interest in Bacillus thuringiensis, an entomopathogenic bacterium capable of replacing chemical pesticides. However, the possibility of pests developing resistance to a particular strain may impair its use, and there is a need to identify novel strains of this species as potential commercial biopesticides. B. thuringiensis sv. israelensis is one of the most successful serovars, widely commercialized for its activity against black fly and mosquito larvae. In this study, we isolated, characterized, and sequenced a new Lebanese B. thuringiensis sv. israelensis isolate, strain AR23. Compared to the commercialized reference strain AM65-52 (Vectobac^®, Sumitomo), AR23 showed an increased activity against several mosquito species. The genomic analysis revealed that this strain, compared to AM65-52, possesses a simplified plasmid content and an additional functional cry4Ba coding gene that most likely accounts for the increased effectiveness of this strain in mosquito larvae killing.

Repertoire of the Bacillus thuringiensis Virulence Factors Unrelated to Major Classes of Protein Toxins and Its Role in Specificity of Host-Pathogen Interactions

Bacillus thuringiensis (Bt) is a Gram-positive soil bacteria that infects invertebrates, predominantly of Arthropoda phylum. Due to its immense host range Bt has become a leading producer of biopesticides applied both in biotechnology and agriculture. Cytotoxic effect of Bt, as well as its host specificity, are commonly attributed either to proteinaceous crystal parasporal toxins (Cry and Cyt) produced by bacteria in a stationary phase or to soluble toxins of Vip and Sip families secreted by vegetative cells. At the same time, numerous non-toxin virulence factors of Bt have been discovered, including metalloproteases, chitinases, aminopolyol antibiotics and nucleotide-mimicking moieties. These agents act at each stage of the B. thuringiensis invasion and contribute to cytotoxic properties of Bt strains enhancing toxin activity, ensuring host immune response evasion and participating in extracellular matrix degeneration. In this review we attempt to classify Bt virulence factors unrelated to major groups of protein toxins and discuss their putative role in the establishment of Bt specificity to various groups of insects.

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.