Interspersed DNA repeats bcr1-bcr18 of Bacillus cereus group bacteria form three distinct groups with different evolutionary and functional patterns

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.

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.

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.

Comparative Bioinformatics and Experimental Analysis of the Intergenic Regulatory Regions of Bacillus cereus hbl and nhe Enterotoxin Operons and the Impact of CodY on Virulence Heterogeneity

Bacillus cereus is a food contaminant with greatly varying enteropathogenic potential. Almost all known strains harbor the genes for at least one of the three enterotoxins Nhe, Hbl, and CytK. While some strains show no cytotoxicity, others have caused outbreaks, in rare cases even with lethal outcome. The reason for these differences in cytotoxicity is unknown. To gain insight into the origin of enterotoxin expression heterogeneity in different strains, the architecture and role of 5′ intergenic regions (5′ IGRs) upstream of the nhe and hbl operons was investigated. In silico comparison of 142 strains of all seven phylogenetic groups of B. cereus sensu lato proved the presence of long 5′ IGRs upstream of the nheABC and hblCDAB operons, which harbor recognition sites for several transcriptional regulators, including the virulence regulator PlcR, redox regulators ResD and Fnr, the nutrient-sensitive regulator CodY as well as the master regulator for biofilm formation SinR. By determining transcription start sites, unusually long 5′ untranslated regions (5′ UTRs) upstream of the nhe and hbl start codons were identified, which are not present upstream of cytK-1 and cytK-2. Promoter fusions lacking various parts of the nhe and hbl 5′ UTR in B. cereus INRA C3 showed that the entire 331 bp 5′ UTR of nhe is necessary for full promoter activity, while the presence of the complete 606 bp hbl 5′ UTR lowers promoter activity. Repression was caused by a 268 bp sequence directly upstream of the hbl transcription start. Luciferase activity of reporter strains containing nhe and hbl 5′ IGR lux fusions provided evidence that toxin gene transcription is upregulated by the depletion of free amino acids. Electrophoretic mobility shift assays showed that the branched-chain amino acid sensing regulator CodY binds to both nhe and hbl 5′ UTR downstream of the promoter, potentially acting as a nutrient-responsive roadblock repressor of toxin gene transcription. PlcR binding sites are highly conserved among all B. cereus sensu lato strains, indicating that this regulator does not significantly contribute to the heterogeneity in virulence potentials. The CodY recognition sites are far less conserved, perhaps conferring varying strengths of CodY binding, which might modulate toxin synthesis in a strain-specific manner.

Single-Cell (Meta-)Genomics of a Dimorphic Candidatus Thiomargarita nelsonii Reveals Genomic Plasticity

The genus Thiomargarita includes the world's largest bacteria. But as uncultured organisms, their physiology, metabolism, and basis for their gigantism are not well understood. Thus, a genomics approach, applied to a single Candidatus Thiomargarita nelsonii cell was employed to explore the genetic potential of one of these enigmatic giant bacteria. The Thiomargarita cell was obtained from an assemblage of budding Ca. T. nelsonii attached to a provannid gastropod shell from Hydrate Ridge, a methane seep offshore of Oregon, USA. Here we present a manually curated genome of Bud S10 resulting from a hybrid assembly of long Pacific Biosciences and short Illumina sequencing reads. With respect to inorganic carbon fixation and sulfur oxidation pathways, the Ca. T. nelsonii Hydrate Ridge Bud S10 genome was similar to marine sister taxa within the family Beggiatoaceae. However, the Bud S10 genome contains genes suggestive of the genetic potential for lithotrophic growth on arsenite and perhaps hydrogen. The genome also revealed that Bud S10 likely respires nitrate via two pathways: a complete denitrification pathway and a dissimilatory nitrate reduction to ammonia pathway. Both pathways have been predicted, but not previously fully elucidated, in the genomes of other large, vacuolated, sulfur-oxidizing bacteria. Surprisingly, the genome also had a high number of unusual features for a bacterium to include the largest number of metacaspases and introns ever reported in a bacterium. Also present, are a large number of other mobile genetic elements, such as insertion sequence (IS) transposable elements and miniature inverted-repeat transposable elements (MITEs). In some cases, mobile genetic elements disrupted key genes in metabolic pathways. For example, a MITE interrupts hupL, which encodes the large subunit of the hydrogenase in hydrogen oxidation. Moreover, we detected a group I intron in one of the most critical genes in the sulfur oxidation pathway, dsrA. The dsrA group I intron also carried a MITE sequence that, like the hupL MITE family, occurs broadly across the genome. The presence of a high degree of mobile elements in genes central to Thiomargarita's core metabolism has not been previously reported in free-living bacteria and suggests a highly mutable genome.

Intrastrand triplex DNA repeats in bacteria: a source of genomic instability

Repetitive nucleic acid sequences are often prone to form secondary structures distinct from B-DNA. Prominent examples of such structures are DNA triplexes. We observed that certain intrastrand triplex motifs are highly conserved and abundant in prokaryotic genomes. A systematic search of 5246 different prokaryotic plasmids and genomes for intrastrand triplex motifs was conducted and the results summarized in the ITxF database available online at http://bioinformatics.uni-konstanz.de/utils/ITxF/. Next we investigated biophysical and biochemical properties of a particular G/C-rich triplex motif (TM) that occurs in many copies in more than 260 bacterial genomes by CD and nuclear magnetic resonance spectroscopy as well as in vivo footprinting techniques. A characterization of putative properties and functions of these unusually frequent nucleic acid motifs demonstrated that the occurrence of the TM is associated with a high degree of genomic instability. TM-containing genomic loci are significantly more rearranged among closely related Escherichia coli strains compared to control sites. In addition, we found very high frequencies of TM motifs in certain Enterobacteria and Cyanobacteria that were previously described as genetically highly diverse. In conclusion we link intrastrand triplex motifs with the induction of genomic instability. We speculate that the observed instability might be an adaptive feature of these genomes that creates variation for natural selection to act upon.

Is Cytotoxin K from Bacillus cereus a bona fide enterotoxin?

Cytotoxin K (CytK) produced by Bacillus cereus s.l. has generally been considered to be associated with the foodborne diarrhoeal syndrome. Two distinct variants of CytK have been reported: CytK-1 from Bacillus cytotoxicus and CytK-2 from B. cereus. In order to determine whether CytK plays a significant role in the diarrhoeal disease, the occurrence of cytK genes was assessed among 390 B. cereus isolates with different origins including clinical and food poisoning samples and was found to be 46%. Interestingly, the cytK occurrence was slightly lower in food poisoning and clinical isolates than in environmental samples. Seventy cytK-2 positive strains (including 28 isolates from foodborne outbreaks) were then selected in order to assess their genetic diversity. A genetic dendrogram based on the cytK-2 sequences of these 70 strains and on two cytK-1 sequences from strains NVH 391-98 and 883-00 showed an important diversity. However, no strain clustering according to the origin or source of isolation was observed. These observations were confirmed by Multi-Locus Sequences Typing (MLST) based on five different loci of housekeeping genes (ccpA, recF, sucC, purF and gdpD) for which no grouping of foodborne outbreak strains could be identified. Therefore, the choice of cytK as virulence factor for the diarrhoeal pathotype does not seem to be relevant per se, even though the involvement of CytK in the diarrhoeal syndrome cannot be fully excluded. Potential synergistic effects between CytK and other virulence factors, together with their potential variable expression levels should be further investigated.

Global mRNA decay analysis at single nucleotide resolution reveals segmental and positional degradation patterns in a Gram-positive bacterium

Background

Recent years have shown a marked increase in the use of next-generation sequencing technologies for quantification of gene expression (RNA sequencing, RNA-Seq). The expression level of a gene is a function of both its rate of transcription and RNA decay, and the influence of mRNA decay rates on gene expression in genome-wide studies of Gram-positive bacteria is under-investigated.

Results

In this work, we employed RNA-Seq in a genome-wide determination of mRNA half-lives in the Gram-positive bacterium Bacillus cereus. By utilizing a newly developed normalization protocol, RNA-Seq was used successfully to determine global mRNA decay rates at the single nucleotide level. The analysis revealed positional degradation patterns, with mRNAs being degraded from both ends of the molecule, indicating that both 5' to 3' and 3' to 5' directions of RNA decay are present in B. cereus. Other operons showed segmental degradation patterns where specific ORFs within polycistrons were degraded at variable rates, underlining the importance of RNA processing in gene regulation. We determined the half-lives for more than 2,700 ORFs in B. cereus ATCC 10987, ranging from less than one minute to more than fifteen minutes, and showed that mRNA decay rate correlates globally with mRNA expression level, GC content, and functional class of the ORF.

Conclusions

To our knowledge, this study presents the first global analysis of mRNA decay in a bacterium at single nucleotide resolution. We provide a proof of principle for using RNA-Seq in bacterial mRNA decay analysis, revealing RNA processing patterns at the single nucleotide level.

Genome of alkaliphilic Bacillus pseudofirmus OF4 reveals adaptations that support the ability to grow in an external pH range from 7.5 to 11.4

Bacillus pseudofirmus OF4 is an extreme but facultative alkaliphile that grows non-fermentatively in a pH range from 7.5 to above 11.4 and can withstand large sudden increases in external pH. It is a model organism for studies of bioenergetics at high pH, at which energy demands are higher than at neutral pH because both cytoplasmic pH homeostasis and ATP synthesis require more energy. The alkaliphile also tolerates a cytoplasmic pH > 9.0 at external pH values at which the pH homeostasis capacity is exceeded, and manages other stresses that are exacerbated at alkaline pH, e.g. sodium, oxidative and cell wall stresses. The genome of B. pseudofirmus OF4 includes two plasmids that are lost from some mutants without viability loss. The plasmids may provide a reservoir of mobile elements that promote adaptive chromosomal rearrangements under particular environmental conditions. The genome also reveals a more acidic pI profile for proteins exposed on the outer surface than found in neutralophiles. A large array of transporters and regulatory genes are predicted to protect the alkaliphile from its overlapping stresses. In addition, unanticipated metabolic versatility was observed, which could ensure requisite energy for alkaliphily under diverse conditions.

Impact of small repeat sequences on bacterial genome evolution

Intergenic regions of prokaryotic genomes carry multiple copies of terminal inverted repeat (TIR) sequences, the nonautonomous miniature inverted-repeat transposable element (MITE). In addition, there are the repetitive extragenic palindromic (REP) sequences that fold into a small stem loop rich in G–C bonding. And the clustered regularly interspaced short palindromic repeats (CRISPRs) display similar small stem loops but are an integral part of a complex genetic element. Other classes of repeats such as the REP2 element do not have TIRs but show other signatures. With the current availability of a large number of whole-genome sequences, many new repeat elements have been discovered. These sequences display diverse properties. Some show an intimate linkage to integrons, and at least one encodes a small RNA. Many repeats are found fused with chromosomal open reading frames, and some are located within protein coding sequences. Small repeat units appear to work hand in hand with the transcriptional and/or post-transcriptional apparatus of the cell. Functionally, they are multifaceted, and this can range from the control of gene expression, the facilitation of host/pathogen interactions, or stimulation of the mammalian immune system. The CRISPR complex displays dramatic functions such as an acquired immune system that defends against invading viruses and plasmids. Evolutionarily, mobile repeat elements may have influenced a cycle of active versus inactive genes in ancestral organisms, and some repeats are concentrated in regions of the chromosome where there is significant genomic plasticity. Changes in the abundance of genomic repeats during the evolution of an organism may have resulted in a benefit to the cell or posed a disadvantage, and some present day species may reflect a purification process. The diverse structure, eclectic functions, and evolutionary aspects of repeat elements are described.