Multiple-locus sequence typing analysis of Bacillus cereus and Bacillus thuringiensis reveals separate clustering and a distinct population structure of psychrotrophic strains

New insights into Bacillus cytotoxicus sources, screening, toxicity, and persistence in food production facilities

Bacillus cytotoxicus is a thermotolerant member of the Bacillus cereus group. It has been linked to rare, but at times fatal cases of diarrheal disease and might be missed at routine diagnostic screening temperatures commonly used for the B. cereus group. The pathogen is mostly found on dehydrated foods containing potato starch or insects. How it enters the food chain or whether it persists in food producing environments is largely unknown. Increased consumption of insects and convenience foods in Europe and the lack of information on the persistence of B. cytotoxicus in food environments and its virulence demand for further characterization. In this study, we aimed to obtain a better understanding of i) the food sources of B. cytotoxicus, ii) screening temperatures needed for its isolation from food matrices, iii) cytotoxicity of the organism, and iv) its ecological niche and potential epidemiological links. To this end, 112 food samples were collected, with a focus on foods exhibiting low water activity. The samples were screened for B. cytotoxicus at 42 °C and at 50 °C. Presumptive isolates were characterized by cytK-1 toxin gene PCR for differentiation of B. cytotoxicus from other B. cereus group members. Vero cell cytotoxicity assays were performed, and selected isolates were sequenced. Our results show that screening at 42 °C might be insufficient for detecting B. cytotoxicus in foods that harbor other less thermophilic Bacillus species. When screening at 50 °C, B. cytotoxicus was detected in 23% of the food samples (n = 26 isolates). The highest prevalence was detected in mashed potato products (82%) and potato flakes (67%). In contrast, a wide range of products not containing any potato ingredients did not yield B. cytotoxicus isolates. All B. cytotoxicus isolates exhibited either low or no detectable cytotoxicity. WGS analysis revealed that a highly toxic isolate is closely related to the French outbreak strain NVH 391-98. In addition, we could show that two isolates sampled 5 years apart from the same production facility only differed by seven SNPs, making it likely that B. cytotoxicus is able to persist in production facilities over a long time. Interestingly, the reoccurring strain possessed an additional plasmid and did not show cytotoxic potential when re-isolated after 5 years.

A Genomic Analysis of the Bacillus Bacteriophage Kirovirus kirovense Kirov and Its Ability to Preserve Milk

Bacteriophages are widely recognized as alternatives to traditional antibiotics commonly used in the treatment of bacterial infection diseases and in the food industry, as phages offer a potential solution in combating multidrug-resistant bacterial pathogens. In this study, we describe a novel bacteriophage, Kirovirus kirovense Kirov, which infects members of the Bacillus cereus group. Kirovirus kirovense Kirov is a broad-host-range phage belonging to the Caudoviricetes class. Its chromosome is a linear 165,667 bp double-stranded DNA molecule that contains two short, direct terminal repeats, each 284 bp long. According to bioinformatics predictions, the genomic DNA contains 275 protein-coding genes and 5 tRNA genes. A comparative genomic analysis suggests that Kirovirus kirovense Kirov is a novel species within the Kirovirus genus, belonging to the Andregratiavirinae subfamily. Kirovirus kirovense Kirov demonstrates the ability to preserve and decontaminate B. cereus from cow milk when present in milk at a concentration of 104 PFU/mL. After 4 h of incubation with the phage, the bacterial titer drops from 105 to less than 102 CFU/mL.

Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria

Core genome multilocus sequence typing (cgMLST) employs a strategy where the set of orthologous genes common to all members of a group of organisms are used for phylogenetic analysis of the group members. The Bacillus cereus group consists of species with pathogenicity towards insect species as well as warm-blooded animals including humans. While B. cereus is an opportunistic pathogen linked to a range of human disease conditions, including emesis and diarrhoea, Bacillus thuringiensis is an entomopathogenic species with toxicity toward insect larvae, and therefore used as a biological pesticide worldwide. Bacillus anthracis is a classical obligate pathogen causing anthrax, an acute lethal condition in herbivores as well as humans, and which is endemic in many parts of the world. The group also includes a range of additional species, and B. cereus group bacteria have been subject to analysis with a wide variety of phylogenetic typing systems. Here we present, based on analyses of 173 complete genomes from B. cereus group species available in public databases, the identification of a set of 1568 core genes which were used to create a core genome multilocus typing scheme for the group which is implemented in the PubMLST system as an open online database freely available to the community. The new cgMLST system provides unprecedented resolution over existing phylogenetic analysis schemes covering the B. cereus group.

A Comparative Analysis of the Core Proteomes within and among the Bacillus subtilis and Bacillus cereus Evolutionary Groups Reveals the Patterns of Lineage- and Species-Specific Adaptations

By integrating phylogenomic and comparative analyses of 1104 high-quality genome sequences, we identify the core proteins and the lineage-specific fingerprint proteins of the various evolutionary clusters (clades/groups/species) of the Bacillus genus. As fingerprints, we denote those core proteins of a certain lineage that are present only in that particular lineage and absent in any other Bacillus lineage. Thus, these lineage-specific fingerprints are expected to be involved in particular adaptations of that lineage. Intriguingly, with a few notable exceptions, the majority of the Bacillus species demonstrate a rather low number of species-specific fingerprints, with the majority of them being of unknown function. Therefore, species-specific adaptations are mostly attributed to highly unstable (in evolutionary terms) accessory proteomes and possibly to changes at the gene regulation level. A series of comparative analyses consistently demonstrated that the progenitor of the Cereus Clade underwent an extensive genomic expansion of chromosomal protein-coding genes. In addition, the majority (76–82%) of the B. subtilis proteins that are essential or play a significant role in sporulation have close homologs in most species of both the Subtilis and the Cereus Clades. Finally, the identification of lineage-specific fingerprints by this study may allow for the future development of highly specific vaccines, therapeutic molecules, or rapid and low-cost molecular tests for species identification.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan-genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well-studied Bacillus cereus sensu lato group, we identified clade-specific core genes (present in all clade members) and then used clade-specific allelic diversity to identify genes under purifying and diversifying selection. Clade-specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade-specific core or accessory genomes. Genes under purifying selection were often clade-specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade-specific selection, allowing ready identification of strains and their ecological niche.

The Membrane Proteome of Spores and Vegetative Cells of the Food-Borne Pathogen Bacillus cereus

Membrane proteins are fascinating since they play an important role in diverse cellular functions and constitute many drug targets. Membrane proteins are challenging to analyze. The spore, the most resistant form of known life, harbors a compressed inner membrane. This membrane acts not only as a barrier for undesired molecules but also as a scaffold for proteins involved in signal transduction and the transport of metabolites during spore germination and subsequent vegetative growth. In this study, we adapted a membrane enrichment method to study the membrane proteome of spores and cells of the food-borne pathogen Bacillus cereus using quantitative proteomics. Using bioinformatics filtering we identify and quantify 498 vegetative cell membrane proteins and 244 spore inner membrane proteins. Comparison of vegetative and spore membrane proteins showed there were 54 spore membrane-specific and 308 cell membrane-specific proteins. Functional characterization of these proteins showed that the cell membrane proteome has a far larger number of transporters, receptors and proteins related to cell division and motility. This was also reflected in the much higher expression level of many of these proteins in the cellular membrane for those proteins that were in common with the spore inner membrane. The spore inner membrane had specific expression of several germinant receptors and spore-specific proteins, but also seemed to show a preference towards the use of simple carbohydrates like glucose and fructose owing to only expressing transporters for these. These results show the differences in membrane proteome composition and show us the specific proteins necessary in the inner membrane of a dormant spore of this toxigenic spore-forming bacterium to survive adverse conditions.

Development of a database and standardized approach for rpoB sequence-based subtyping and identification of aerobic spore-forming Bacillales

Aerobic spore-forming Bacillales are a highly diverse and ubiquitous group that includes organisms that cause foodborne illnesses and food spoilage. Classical microbiological and biochemical identification of members of the order Bacillales represents a challenge due to the diversity of organisms in this group as well as the fact that the phenotypic-based taxonomic assignment of some named species in this group is not consistent with their phylogenomic characteristics. DNA-sequencing-based tools, on the other hand, can be fast and cost-effective, and can provide for a more reliable identification and characterization of Bacillales isolates. In comparison to 16S rDNA, rpoB was shown to better discriminate between Bacillales isolates and to allow for improved taxonomic assignment to the species level. However, the lack of a publicly accessible rpoB database, as well as the lack of standardized protocols for rpoB-based typing and strain identification, is a major challenge. Here, we report (i) the curation of a DNA sequence database for rpoB-based subtype classification of Bacillales isolates; (ii) the development of standardized protocols for generating rpoB sequence data, and a scheme for rpoB-based initial taxonomic identification of Bacillales isolates at the species level; and (iii) the integration of the database in a publicly accessible online platform that allows for the analysis of rpoB sequence data from uncharacterized Bacillales isolates. Specifically, we curated a database of DNA sequences for a 632-nt internal variable region within the rpoB gene from representative Bacillales reference type strains and a large number of isolates that we have previously isolated and characterized through multiple projects. As of May 21, 2021, the rpoB database contained more than 8350 rpoB sequences representing 1902 distinct rpoB allelic types that can be classified into 160 different genera. The database also includes 1129 rpoB sequences for representative Bacillales reference type strains as available on May 21, 2021 in the NCBI database. The rpoB database is integrated into the online Food Microbe Tracker platform (www.foodmicrobetracker.com) and can be queried using the integrated BLAST tool to initially subtype and taxonomically identify aerobic and facultative anaerobic spore-formers. While whole-genome sequencing is increasingly used in bacterial taxonomy, the rpoB sequence-based identification scheme described here provides a valuable tool as it allows for rapid and cost-effective initial isolate characterization, which can help to identify and characterize foodborne pathogens and food spoilage bacteria. In addition, the database and primers described here can also be adopted for metagenomics approaches that include rpoB as a target, improving discriminatory power and identification over what can be achieved using 16S rDNA as a target.

Effects of Bacillus coagulans as an adjunct starter culture on yogurt quality and storage

Bacillus coagulans has been widely studied for its probiotic properties. Therefore, identifying a strain that can be used as an adjunct starter culture for yogurt production would have commercial value. In this study, 30 B. coagulans strains were isolated from vegetable samples from 11 provinces or autonomous regions in China, and their pan-genomic and phylogenetic characteristics were analyzed. Phylogenetic analysis categorized 30 strains into 4 different subphylotypes, including subtype I (11 isolates), subtype II (7 isolates), subtype III (11 isolates), and subtype IV (1 isolate). Four B. coagulans strains (B. coagulans-70, B. coagulans-78, B. coagulans-79, and B. coagulans-100) were randomly selected from each subphylotype of the phylogenetic tree as adjunct starter cultures. Compared with the other tested strains, B. coagulans-70 showed the highest count in yogurt at the end of the manufacturing period. Comparative genome analysis indicated that the different bacterial levels of B. coagulans strains in yogurt may be associated with the abundance of genes related to carbohydrate transport and metabolism (e.g., sucrose utilization). Finally, differences in texture and volatile flavor compound profiles were observed between the yogurt samples. Compared with the other groups, the addition of B. coagulans-70 exerted a positive effect on the appearance and texture of yogurt products. Volatile analysis showed increased quantities of 2-heptanone, 2-nonanone, amyl alcohol, and 2-hydroxy-3-pentanone in the B. coagulans-70 group compared with control yogurts. These results above combined with the results of a sensory evaluation indicated that B. coagulans-70 is the most suitable strain for further use in functional dairy product development.

The Food Poisoning Toxins of Bacillus cereus

Bacillus cereus is a ubiquitous soil bacterium responsible for two types of food-associated gastrointestinal diseases. While the emetic type, a food intoxication, manifests in nausea and vomiting, food infections with enteropathogenic strains cause diarrhea and abdominal pain. Causative toxins are the cyclic dodecadepsipeptide cereulide, and the proteinaceous enterotoxins hemolysin BL (Hbl), nonhemolytic enterotoxin (Nhe) and cytotoxin K (CytK), respectively. This review covers the current knowledge on distribution and genetic organization of the toxin genes, as well as mechanisms of enterotoxin gene regulation and toxin secretion. In this context, the exceptionally high variability of toxin production between single strains is highlighted. In addition, the mode of action of the pore-forming enterotoxins and their effect on target cells is described in detail. The main focus of this review are the two tripartite enterotoxin complexes Hbl and Nhe, but the latest findings on cereulide and CytK are also presented, as well as methods for toxin detection, and the contribution of further putative virulence factors to the diarrheal disease.

Characterization of Bacillus cereus Group Isolates From Human Bacteremia by Whole-Genome Sequencing

Members of the Bacillus cereus group are spore-forming organisms commonly associated with food poisoning and intestinal infections. Moreover, some strains of the group (i.e., B. cereus sensu stricto and Bacillus thuringiensis) can cause bacteremia in humans, mainly in immunocompromised individuals. Here we performed the genetic characterization of 17 human clinical strains belonging to B. cereus group isolated from blood culture. The whole-genome sequencing (WGS) revealed that the isolates were closely related to B. cereus sensu stricto and B. thuringiensis-type strain. Multilocus sequence typing analysis performed on the draft genome revealed the genetic diversity of our isolates, which were assigned to different sequence types. Based on panC nucleotide sequence, the isolates were grouped in the phylogenetic groups III and IV. The NHE, cer, and inhA gene cluster, entA, entFM, plcA, and plcB, were the most commonly detected virulence genes. Although we did not assess the ability to generate biofilm by phenotypic tests, we verified the prevalence of biofilm associated genes using an in silico approach. A high prevalence of pur gene cluster, xerC, clpY, codY, tasA, sipW, sinI, and sigB genes, was found. Genes related to the resistance to penicillin, trimethoprim, and ceftriaxone were identified in most of the isolates. Intriguingly, the majority of these virulence and AMR genes appeared to be evenly distributed among B. cereus s.s. isolates, as well as closely related to B. thuringiensis isolates. We showed the WGS represents a good approach to rapidly characterize B. cereus group strains, being able to give useful information about genetic epidemiology, the presence of virulence and antimicrobial genes, and finally about the potential hazard related to this underestimated risk.

Divergence in environmental adaptation between terrestrial clades of the Bacillus cereus group

The Bacillus cereus group encompasses beneficial and harmful species in diverse niches and has a much debated taxonomy. Investigating whether selection has led to ecological divergence between phylogenetic clades can help understand the basis of speciation, and has implications for predicting biological safety across this group. Using three most terrestrial species in this group (B. cereus, Bacillus thuringiensis and Bacillus mycoides) we charactererized ecological specialization in terms of resource use, thermal adaptation and fitness in different environmental conditions and tested whether taxonomic species or phylogenetic clade best explained phenotypic variation. All isolates grew vigorously in protein rich media and insect cadavers, but exploitation of soil or plant derived nutrients was similarly weak for all. For B. thuringiensis and B. mycoides, clade and taxonomic species were important predictors of relative fitness in insect infections. Fully psychrotolerant isolates could outcompete B. thuringiensis in insects at low temperature, although psychrotolerance predicted growth in artificial media better than clade. In contrast to predictions, isolates in the Bacillus anthracis clade had sub-optimal growth at 37°C. The common ecological niche in these terrestrial B. cereus species is the ability to exploit protein rich resources such as cadavers. However, selection has led to different phylogenetic groups developing different strategies for accessing this resource. Thus, clades, as well as traditional taxonomic phenotypes, predict biologically important traits.

Polyphasic characterization and identification of the bioremediation agent Bacillus sp. SFC 500-1E

Bacillus sp. SFC 500-1E is used for the effective treatment of tannery effluents since it consistently removes hexavalent chromium from diverse contaminated matrices. The aim of the present study was to complete identification of the strain through a polyphasic characterization, which included the pattern of carbohydrate utilization, fatty acids profile, multilocus sequence analysis, multiplex PCR profile and the analysis of the complete genome sequence. Morpho-physiological and biochemical characterization results and analysis of 16S rRNA sequences were not conclusive. The strain formed a monophyletic clade with B. toyonensis BCT-7112, B. thuringiensis MC28 and B. cereus Rock 1-3. However, genomic comparisons with type strains of B. cereus and B. thuringiensis showed that the isolated belonged to a different species. Results of this study highlight the relevance of the genome sequence of this strain, identified as Bacillus toyonensis SFC 500-1E, to expand knowledge of its bioremediation potential and to explore unknown decontamination activities.

Comparative genomic analysis and mosquito larvicidal activity of four Bacillus thuringiensis serovar israelensis strains

Bacillus thuringiensis serovar israelensis (Bti) is used to control insect vectors of human and animal diseases. In the present study, the toxicity of four strains of Bti, named T0124, T0131, T0137, and T0139, toward Aedes aegypti and Culex quinquefasciatus larvae was analyzed. The T0131 strain showed the highest larvicidal activity against A. aegypti (LC₅₀ = 0.015 µg/ml) and C. quinquefasciatus larvae (LC₅₀ = 0.035 µg/ml) when compared to the other strains. Furthermore, the genomic sequences of the four strains were obtained and compared. These Bti strains had chromosomes sizes of approximately 5.4 Mb with GC contents of ~35% and 5472–5477 putative coding regions. Three small plasmids (5.4, 6.8, and 7.6 kb) and three large plasmids (127, 235, and 359 kb) were found in the extrachromosomal content of all four strains. The SNP-based phylogeny revealed close relationship among isolates from this study and other Bti isolates, and SNPs analysis of the plasmids 127 kb did not reveal any mutations in δ-endotoxins genes. This newly acquired sequence data for these Bti strains may be useful in the search for novel insecticidal toxins to improve existing ones or develop new strategies for the biological control of important insect vectors of human and animal diseases.

Suboptimal Bacillus licheniformis and Bacillus weihenstephanensis Spore Incubation Conditions Increase Heterogeneity of Spore Outgrowth Time

Changes with time of a population of Bacillus weihenstephanensis KBAB4 and Bacillus licheniformis AD978 dormant spores into germinated spores and vegetative cells were followed by flow cytometry, at pH ranges of 4.7 to 7.4 and temperatures of 10°C to 37°C for B. weihenstephanensis and 18°C to 59°C for B. licheniformis Incubation conditions lower than optimal temperatures or pH led to lower proportions of dormant spores able to germinate and extended time of germination, a lower proportion of germinated spores able to outgrow, an extension of their times of outgrowth, and an increase of the heterogeneity of spore outgrowth time. A model based on the strain growth limits was proposed to quantify the impact of incubation temperature and pH on the passage through each physiological stage. The heat treatment temperature or time acted independently on spore recovery. Indeed, a treatment at 85°C for 12 min or at 95°C for 2 min did not have the same impact on spore germination and outgrowth kinetics of B. weihenstephanensis despite the fact that they both led to a 10-fold reduction of the population. Moreover, acidic sporulation pH increased the time of outgrowth 1.2-fold and lowered the proportion of spores able to germinate and outgrow 1.4-fold. Interestingly, we showed by proteomic analysis that some proteins involved in germination and outgrowth were detected at a lower abundance in spores produced at pH 5.5 than in those produced at pH 7.0, maybe at the origin of germination and outgrowth behavior of spores produced at suboptimal pH.IMPORTANCE Sporulation and incubation conditions have an impact on the numbers of spores able to recover after exposure to sublethal heat treatment. Using flow cytometry, we were able to follow at a single-cell level the changes in the physiological states of heat-stressed spores of Bacillus spp. and to discriminate between dormant spores, germinated spores, and outgrowing vegetative cells. We developed original mathematical models that describe (i) the changes with time of the proportion of cells in their different states during germination and outgrowth and (ii) the influence of temperature and pH on the kinetics of spore recovery using the growth limits of the tested strains as model parameters. We think that these models better predict spore recovery after a sublethal heat treatment, a common situation in food processing and a concern for food preservation and safety.

Population structure and toxin gene profiles of Bacillus cereus sensu lato isolated from flour products

Data on the occurrence, population structure and toxinogenic potential of Bacillus cereus sensu lato isolated from flour is essential to enable improved risk assessment. We aimed to provide data on the occurrence of B. cereus sensu lato in flour products at retail level. In addition, we screened the isolates for Bacillus thuringiensis and Bacillus cytotoxicus and determined population structure and toxin gene profiles. We screened 89 flour products for presence of B. cereus sensu lato, resulting in 75 positive samples (84%). We were able to show that the population structure of members of the B. cereus group isolated from flour is highly diverse. Isolates were assigned to panC types II (4%), III (21%), IV (39%) and V (36%). Production of parasporal crystals characteristic for Bacillus thuringiensis was detected in seven isolates assigned to panC type III, IV and V. No B. cytotoxicus were detected. Two of the isolates harbored ces encoding cereulide, which causes the emetic syndrome. Various enterotoxin genes were found, with all isolates harboring nhe, 75% of isolates harboring hbl and 51% of the isolates harboring cytK-2. Our findings suggest that toxinogenic B. cereus sensu lato are common in flour products at retail level.

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.

Bacillus wiedmannii biovar thuringiensis: A Specialized Mosquitocidal Pathogen with Plasmids from Diverse Origins

Bacillus cereus sensu lato also known as B. cereus group is composed of an ecologically diverse bacterial group with an increasing number of related species, some of which are medically or agriculturally important. Numerous efforts have been undertaken to allow presumptive differentiation of B. cereus group species from one another. FCC41 is a Bacillus sp. strain toxic against mosquito species like Aedes aegypti, Aedes (Ochlerotatus) albifasciatus, Culex pipiens, Culex quinquefasciatus, and Culex apicinus, some of them responsible for the transmission of vector-borne diseases. Here, we report the complete genome sequence of FCC41 strain, which consists of one circular chromosome and eight circular plasmids ranging in size from 8 to 490 kb. This strain harbors six crystal protein genes, including cry24Ca, two cry4-like and two cry52-like, a cry41-like parasporin gene and multiple virulence factors. The phylogenetic analysis of the whole-genome sequence of this strain with molecular approaches places this strain into the Bacillus wiedmannii cluster. However, according with phenotypical characteristics such as the mosquitocidal activity due to the presence of Cry proteins found in the parasporal body and cry genes encoded in plasmids of different sizes, indicate that this strain could be renamed as B. wiedmannii biovar thuringiensis strain FCC41.

Genes under positive selection in the core genome of pathogenic Bacillus cereus group members

In this comparative genomics study our aim was to unravel genes under positive selection in the core genome of the Bacillus cereus group. Indeed, the members of this group share close genetic relationships but display a rather large phenotypic and ecological diversity, providing a unique opportunity for studying how genomic changes reflect ecological adaptation during the divergence of a bacterial group. For this purpose, we screened ten completely sequenced genomes of four pathogenic Bacillus species, finding that 254 out of 3093 genes have codon sites with d_N/d_S (ω) values above one. These results remained unchanged after having disentangled the confounding effects of recombination and selection signature in a Bayesian framework. The presumably adaptive nucleotide polymorphisms are distributed over a wide range of biological functions, such as antibiotic resistance, DNA repair, nutrient uptake, metabolism, cell wall assembly and spore structure. Our results indicate that adaptation to animal hosts, whether as pathogens, saprophytes or symbionts, is the major driving force in the evolution of the Bacillus cereus group. Future work should seek to understand the evolutionary dynamics of both core and accessory genes in an integrative framework to ultimately unravel the key networks involved in host adaptation.

Multilocus sequence typing for phylogenetic view and vip gene diversity of Bacillus thuringiensis strains of the Assam soil of North East India

An agriculturally important insecticidal bacterium, Bacillus thuringiensis have been isolated from the soil samples of various part of Assam including the Kaziranga National Park. Previously, the isolates were characterized based on morphology, 16S rDNA sequencing, and the presence of the various classes' crystal protein gene(s). In the present study, the phylogenetic analysis of a few selected isolates was performed by an unambiguous and quick method called the multiple locus sequence typing (MLST). A known B. thuringiensis strain kurstaki 4D4 have been used as a reference strain for MLST. A total of four the MLST locus of housekeeping genes, recF, sucC, gdpD and yhfL were selected. A total of 14 unique sequence types (STs) was identified. A total number of alleles identified for the locus gdpD and sucC was 12, followed by locus yhfL was 11, however, only 6 alleles were detected for the locus recF. The phylogenetic analysis using MEGA 7.0.26 showed three major lineages. Approximately, 87% of the isolates belonged to the STs corresponding to B. thuringiensis, whereas two isolates, BA07 and BA39, were clustered to B. cereus. The isolates were also screened for the diversity of vegetative insecticidal protein (vip) genes. In all, 8 isolates showed the presence of vip1, followed by 7 isolates having vip2 and 6 isolates for vip3 genes. The expression of Vip3A proteins was analyzed by western blot analyses and expression of the Vip3A protein was observed in the isolate BA20. Thus, the phylogenetic relationship and diversity of Bt isolates from Assam soil was established based on MLST, in addition, found isolates having vip genes, which could be used for crop improvement.

Analysis of enterotoxigenic Bacillus cereus strains from dried foods using whole genome sequencing, multi-locus sequence analysis and toxin gene prevalence and distribution using endpoint PCR analysis

Bacillus cereus strains were isolated from dried foods, which included international brands of spices from South East Asia, Mexico and India purchased from several retail stores, samples of powdered infant formula (PIF), medicated fish feed and dietary supplements. The genetic diversity of 64 strains from spices and PIF was determined using a multiplex endpoint PCR assay designed to identify hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Thirteen different B. cereus toxigenic gene patterns or profiles were identified among the strains. Randomly selected B. cereus strains were sequenced and compared with reference Genomic Groups from National Center Biotechnology Information using bioinformatics tools. A comprehensive multi-loci sequence analysis (MLSA) was designed using alleles from 25 known MLST genes specifically tailored for use with whole genome assemblies. A cohort of representative genomes of strains from a few FDA regulated commodities like dry foods and medicated fish feed was used to demonstrate the utility of the 25-MLSA approach for rapid clustering and identification of Genome Groups. The analysis clustered the strains from medicated fish feed, dry foods, and dietary supplements into phylogenetically-related groups. 25-MLSA also pointed to a greater diversity of B. cereus strains from foods and feed than previously recognized. Our integrated approach of toxin gene PCR, and to our knowledge, whole genome sequencing (WGS) based sequence analysis, may be the first of its kind that demonstrates enterotoxigenic potential and genomic diversity in parallel.

Selection procedure of bioprotective cultures for their combined use with High Pressure Processing to control spore-forming bacteria in cooked ham

High Pressure Processing (HPP) and biopreservation can contribute to food safety by inactivation of bacterial contaminants. However these treatments are inefficient against bacterial endospores. Moreover, HPP can induce spore germination. The objective of this study was to select lactic acid bacteria strains to be used as bioprotective cultures, to control vegetative cells of spore-forming bacteria in ham after application of HPP. A collection of 63 strains of various origins was screened for their antagonistic activity against spore-forming Bacillus and Clostridium species and their ability to resist to HPP. Some safety requirements should also be considered prior to their introduction into the food chain. Hence, the selection steps included the assessment of biogenic amine production and antibiotic resistance. No strain produced histamine above the threshold detection level of 50 ppm. From the assessment of antibiotic resistance against nine antibiotics, 14 susceptible strains were kept. Antagonistic action of the 14 strains was then assessed by the well diffusion method against pathogenic or spoilage spore-forming species as Bacillus cereus, Clostridium sp. like botulinum, Clostridium frigidicarnis, and Clostridium algidicarnis. One Lactobacillus curvatus strain and one Lactococcus lactis strain were ultimately selected for their widest inhibitory spectrum and their potential production of bacteriocin. A Lactobacillus plantarum strain was included as control. Their resistance to HPP and ability to regrow during chilled storage was then assessed in model ham liquid medium. Treatments of pressure intensities of 400, 500, and 600 MPa, and durations of 1, 3, 6, and 10 min were applied. After treatment, cultures were incubated at 8 °C during 30 days. Inactivation curves were then fitted by using a reparameterized Weibull model whereas growth curves were modelled with a logistic model. Although the two Lactobacillus strains were more resistant than L. lactis to HPP, the latter was the only strain able to regrow following HPP. The absence of biogenic amine production of this strain after growth on diced cube cooked ham was also shown. In conclusion this L. lactis strain could be selected as representing the best candidate for a promising preservative treatment combining biopreservation and HPP to control spore-forming bacteria in cooked ham.

Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co-evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large-scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co-evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.

Rapid, High-Throughput Identification of Anthrax-Causing and Emetic Bacillus cereus Group Genome Assemblies via BTyper, a Computational Tool for Virulence-Based Classification of Bacillus cereus Group Isolates by Using Nucleotide Sequencing Data

The Bacillus cereus group comprises nine species, several of which are pathogenic. Differentiating between isolates that may cause disease and those that do not is a matter of public health and economic importance, but it can be particularly challenging due to the high genomic similarity within the group. To this end, we have developed BTyper, a computational tool that employs a combination of (i) virulence gene-based typing, (ii) multilocus sequence typing (MLST), (iii) panC clade typing, and (iv) rpoB allelic typing to rapidly classify B. cereus group isolates using nucleotide sequencing data. BTyper was applied to a set of 662 B. cereus group genome assemblies to (i) identify anthrax-associated genes in non-B. anthracis members of the B. cereus group, and (ii) identify assemblies from B. cereus group strains with emetic potential. With BTyper, the anthrax toxin genes cya, lef, and pagA were detected in 8 genomes classified by the NCBI as B. cereus that clustered into two distinct groups using k-medoids clustering, while either the B. anthracis poly-γ-d-glutamate capsule biosynthesis genes capABCDE or the hyaluronic acid capsule hasA gene was detected in an additional 16 assemblies classified as either B. cereus or Bacillus thuringiensis isolated from clinical, environmental, and food sources. The emetic toxin genes cesABCD were detected in 24 assemblies belonging to panC clades III and VI that had been isolated from food, clinical, and environmental settings. The command line version of BTyper is available at https://github.com/lmc297/BTyper. In addition, BMiner, a companion application for analyzing multiple BTyper output files in aggregate, can be found at https://github.com/lmc297/BMiner.

IMPORTANCEBacillus cereus is a foodborne pathogen that is estimated to cause tens of thousands of illnesses each year in the United States alone. Even with molecular methods, it can be difficult to distinguish nonpathogenic B. cereus group isolates from their pathogenic counterparts, including the human pathogen Bacillus anthracis, which is responsible for anthrax, as well as the insect pathogen B. thuringiensis. By using the variety of typing schemes employed by BTyper, users can rapidly classify, characterize, and assess the virulence potential of any isolate using its nucleotide sequencing data.

In defense of Bacillus thuringiensis, the safest and most successful microbial insecticide available to humanity - a response to EFSA

The Bacillus cereus group contains vertebrate pathogens such as B. anthracis and B. cereus and the invertebrate pathogen B. thuringiensis (Bt). Microbial biopesticides based on Bt are widely recognised as being among the safest and least environmentally damaging insecticidal products available. Nevertheless, a recent food-poisoning incident prompted a European Food Safety Authority review which argued that Bt poses a health risk equivalent to B. cereus, a causative agent of diarrhoea. However, a critical examination of available data, and this latest incident, provides no solid evidence that Bt causes diarrhoea. Although relatively high levels of B. cereus-like spores can occur in foods, genotyping demonstrates that these are predominantly naturally occurring strains rather than biopesticides. Moreover, MLST genotyping of >2000 isolates show that biopesticide genotypes have never been isolated from any clinical infection. MLST data demonstrate that B. cereus group is heterogeneous and formed of distinct clades with substantial differences in biology, ecology and host association. The group posing the greatest risk (the anthracis clade) is distantly related to the clade containing all biopesticides. These recent data support the long-held view that Bt and especially the strains used in Bt biopesticides are very safe for humans.

Multilocus sequence typing, biochemical and antibiotic resistance characterizations reveal diversity of North American strains of the honey bee pathogen Paenibacillus larvae

Paenibacillus larvae is a Gram positive bacterium and the causative agent of the most widespread fatal brood disease of honey bees, American foulbrood (AFB). A total of thirty-three independent Paenibacillus larvae isolates from various geographical origins in North America and five reference strains were investigated for genetic diversity using multilocus sequence typing (MLST). This technique is regarded to be a powerful tool for epidemiological studies of pathogenic bacteria and is widely used in genotyping assays. For MLST, seven housekeeping gene loci, ilvD (dihydroxy-acid dyhydrogenase), tri (triosephosphate isomerase), purH (phospharibosyl-aminoimidazolecarboxamide), recF (DNA replication and repair protein), pyrE (orotate phosphoribosyltransferase), sucC (succinyl coenzyme A synthetase β subunit) and glpF (glycerol uptake facilitator protein) were studied and applied for primer designs. Previously, ERIC type DNA fingerprinting was applied to these same isolates and the data showed that almost all represented the ERIC I type, whereas using BOX-PCR gave an indication of more diversity. All isolates were screened for resistance to four antibiotics used by U.S. beekeepers, showing extensive resistance to tetracycline and the first records of resistance to tylosin and lincomycin. Our data highlight the intraspecies relationships of P. larvae and the potential application of MLST methods in enhancing our understanding of epidemiological relationships among bacterial isolates of different origins.

The Bacterial Sequential Markov Coalescent

Bacteria can exchange and acquire new genetic material from other organisms directly and via the environment. This process, known as bacterial recombination, has a strong impact on the evolution of bacteria, for example, leading to the spread of antibiotic resistance across clades and species, and to the avoidance of clonal interference. Recombination hinders phylogenetic and transmission inference because it creates patterns of substitutions (homoplasies) inconsistent with the hypothesis of a single evolutionary tree. Bacterial recombination is typically modeled as statistically akin to gene conversion in eukaryotes, i.e., using the coalescent with gene conversion (CGC). However, this model can be very computationally demanding as it needs to account for the correlations of evolutionary histories of even distant loci. So, with the increasing popularity of whole genome sequencing, the need has emerged for a faster approach to model and simulate bacterial genome evolution. We present a new model that approximates the coalescent with gene conversion: the bacterial sequential Markov coalescent (BSMC). Our approach is based on a similar idea to the sequential Markov coalescent (SMC)-an approximation of the coalescent with crossover recombination. However, bacterial recombination poses hurdles to a sequential Markov approximation, as it leads to strong correlations and linkage disequilibrium across very distant sites in the genome. Our BSMC overcomes these difficulties, and shows a considerable reduction in computational demand compared to the exact CGC, and very similar patterns in simulated data. We implemented our BSMC model within new simulation software FastSimBac. In addition to the decreased computational demand compared to previous bacterial genome evolution simulators, FastSimBac provides more general options for evolutionary scenarios, allowing population structure with migration, speciation, population size changes, and recombination hotspots. FastSimBac is available from https://bitbucket.org/nicofmay/fastsimbac, and is distributed as open source under the terms of the GNU General Public License. Lastly, we use the BSMC within an Approximate Bayesian Computation (ABC) inference scheme, and suggest that parameters simulated under the exact CGC can correctly be recovered, further showcasing the accuracy of the BSMC. With this ABC we infer recombination rate, mutation rate, and recombination tract length of Bacillus cereus from a whole genome alignment.

Microbial communities associated with the anthropogenic, highly alkaline environment of a saline soda lime, Poland

Soda lime is a by-product of the Solvay soda process for the production of sodium carbonate from limestone and sodium chloride. Due to a high salt concentration and alkaline pH, the lime is considered as a potential habitat of haloalkaliphilic and haloalkalitolerant microbial communities. This artificial and unique environment is nutrient-poor and devoid of vegetation, due in part to semi-arid, saline and alkaline conditions. Samples taken from the surface layer of the lime and from the depth of 2 m (both having pH ~11 and EC_e up to 423 dS m⁻¹) were investigated using culture-based (culturing on alkaline medium) and culture-independent microbiological approaches (microscopic analyses and pyrosequencing). A surprisingly diverse bacterial community was discovered in this highly saline, alkaline and nutrient-poor environment, with the bacterial phyla Proteobacteria (representing 52.8% of the total bacterial community) and Firmicutes (16.6%) showing dominance. Compared to the surface layer, higher bacterial abundance and diversity values were detected in the deep zone, where more stable environmental conditions may occur. The surface layer was dominated by members of the genera Phenylobacterium, Chelativorans and Skermanella, while in the interior layer the genus Fictibacillus was dominant. The culturable aerobic, haloalkaliphilic bacteria strains isolated in this study belonged mostly to the genus Bacillus and were closely related to the species Bacillus pseudofirmus, B. cereus, B. plakortidis, B. thuringensis and B. pumilus.

Genetic and genomic diversity of NheABC locus from Bacillus strains

Non-hemolytic enterotoxin (NHE), a tri-partite, proteinaceous toxin encoded by contiguous nheA, nheB and nheC genes of Bacillus cereus sensu lato (B. cereus s.l.), is considered to be associated with the foodborne diarrheic syndrome. However, B. cereus s.l. includes a number of closely related strains, and the occurrence of NHE among them, and other members of Bacillus is unclear. Consequently, we aimed to determine the distribution and evolution of NHE within Bacillus by confirming the presence of the nheA, B and C sequences and variation within them using published data, and to analyze the genomic and genetic diversity. The phylogenetic tree of NHE proteins (NheA, NheB and NheC) from 81 different B. cereus s.l. strains was constructed. And on the genetic determinants of the NHE toxin did not bring any obvious link between the nheABC genes sequence of a strain and its virulence in the diarrhoeal pathogenesis. Analysis of the genomic diversity of the nheA, B and C loci revealed that their upstream regions were more conserved than the downstream sequences. Multilocus sequence typing schemes (MLST) based on seven concatenated housekeeping genes and nheA, B and C genes of the 75 strains were developed. The neighbor joining phylogenetic tree based on seven housekeeping genes together with nheA, B and C genes was similiar to published Bacillus phylogenetic trees. And on the genetic determinants of the NHE toxin did not bring any obvious link between the nheABC genes sequence of a strain and its virulence in the diarrhoeal pathogenesis.The results indicate that nheA, B and C genes do not affect the diversity of housekeeping genes, and this specific NHE protein does not participate in the diarrheic syndrome.

Comparative genomics of extrachromosomal elements in Bacillus thuringiensis subsp. israelensis

Bacillus thuringiensis subsp. israelensis is one of the most important microorganisms used against mosquitoes. It was intensively studied following its discovery and became a model bacterium of the B. thuringiensis species. Those studies focused on toxin genes, aggregation-associated conjugation, linear genome phages, etc. Recent announcements of genomic sequences of different strains have not been explicitly related to the biological properties studied. We report data on plasmid content analysis of four strains using ultra-high-throughput sequencing. The strains were commercial product isolates, with their putative ancestor and type B. thuringiensis subsp. israelensis strain sequenced earlier. The assembled contigs corresponding to published and novel data were assigned to plasmids described earlier in B. thuringiensis subsp. israelensis and other B. thuringiensis strains. A new 360 kb plasmid was identified, encoding multiple transporters, also found in most of the earlier sequenced strains. Our genomic data show the presence of two toxin-coding plasmids of 128 and 100 kb instead of the reported 225 kb plasmid, a co-integrate of the former two. In two of the sequenced strains, only a 100 kb plasmid was present. Some heterogeneity exists in the small plasmid content and structure between strains. These data support the perception of active plasmid exchange among B. thuringiensis subsp. israelensis strains in nature.

Ecology of Bacillaceae

Members of the family Bacillaceae are among the most robust bacteria on Earth, which is mainly due to their ability to form resistant endospores. This trait is believed to be the key factor determining the ecology of these bacteria. However, they also perform fundamental roles in soil ecology (i.e., the cycling of organic matter) and in plant health and growth stimulation (e.g., via suppression of plant pathogens and phosphate solubilization). In this review, we describe the high functional and genetic diversity that is found within the Bacillaceae (a family of low-G+C% Gram-positive spore-forming bacteria), their roles in ecology and in applied sciences related to agriculture. We then pose questions with respect to their ecological behavior, zooming in on the intricate social behavior that is becoming increasingly well characterized for some members of Bacillaceae. Such social behavior, which includes cell-to-cell signaling via quorum sensing or other mechanisms (e.g., the production of extracellular hydrolytic enzymes, toxins, antibiotics and/or surfactants) is a key determinant of their lifestyle and is also believed to drive diversification processes. It is only with a deeper understanding of cell-to-cell interactions that we will be able to understand the ecological and diversification processes of natural populations within the family Bacillaceae. Ultimately, the resulting improvements in understanding will benefit practical efforts to apply representatives of these bacteria in promoting plant growth as well as biological control of plant pathogens.

The Phylogeny of Bacillus cereus sensu lato

The three main species of the Bacillus cereus sensu lato , B. cereus , B. thuringiensis , and B. anthracis , were recognized and established by the early 1900s because they each exhibited distinct phenotypic traits. B. thuringiensis isolates and their parasporal crystal proteins have long been established as a natural pesticide and insect pathogen. B. anthracis , the etiological agent for anthrax, was used by Robert Koch in the 19th century as a model to develop the germ theory of disease, and B. cereus , a common soil organism, is also an occasional opportunistic pathogen of humans. In addition to these three historical species designations, are three less-recognized and -understood species: B. mycoides , B. weihenstephanensis , and B. pseudomycoides . All of these “species” combined comprise the Bacillus cereus sensu lato group. Despite these apparently clear phenotypic definitions, early molecular approaches to separate the first three by various DNA hybridization and 16S/23S ribosomal sequence analyses led to some “confusion” because there were limited differences to differentiate between these species. These and other results have led to frequent suggestions that a taxonomic change was warranted to reclassify this group to a single species. But the pathogenic properties of B. anthracis and the biopesticide applications of B. thuringiensis appear to “have outweighed pure taxonomic considerations” and the separate species categories are still being maintained. B. cereus sensu lato represents a classic example of a now common bacterial species taxonomic quandary.

Massive horizontal gene transfer, strictly vertical inheritance and ancient duplications differentially shape the evolution of Bacillus cereus enterotoxin operons hbl, cytK and nhe

Background

Bacillus cereus sensu lato comprises eight closely related species including the human pathogens Bacillus anthracis and Bacillus cereus. Within B. cereus sensu lato, chromosomally and plasmid-encoded toxins exist. While plasmid-mediated horizontal gene transfer of the emetic toxin, anthrax and insecticidal toxins is known, evolution of enterotoxin genes within the group has not been studied.

Results

We report draft genome assemblies of 25 strains, a phylogenetic network of 142 strains based on ANI derived from genome sequences and a phylogeny based on whole-genome SNP analysis. The data clearly support subdivision of B. cereus sensu lato into seven phylogenetic groups. While group I, V and VII represent B. pseudomycoides, B. toyonensis and B. cytotoxicus, which are distinguishable at species level (ANI border ≥ 96 %), strains ascribed to the other five species do not match phylogenic groups. The chromosomal enterotoxin operons nheABC and hblCDAB are abundant within B. cereus both isolated from infections and from the environment. While the duplicated hbl variant hbl_a is present in 22 % of all strains investigated, duplication of nheABC is extremely rare (0.02 %) and appears to be phylogenetically unstable. Distribution of toxin genes was matched to a master tree based on seven concatenated housekeeping genes, which depicts species relationships in B. cereus sensu lato as accurately as whole-genome comparisons. Comparison to the phylogeny of enterotoxin genes uncovered ample evidence for horizontal transfer of hbl, cytK and plcR, as well as frequent deletion of both toxins and duplication of hbl. No evidence for nhe deletion was found and stable horizontal transfer of nhe is rare. Therefore, evolution of B. cereus enterotoxin operons is shaped unexpectedly different for yet unknown reasons.

Conclusions

Frequent exchange of the pathogenicity factors hbl, cytK and plcR in B. cereus sensu lato appears to be an important mechanism of B. cereus virulence evolution, including so-called probiotic or non-pathogenic species, which might have consequences for risk assessment procedures. In contrast, exclusively vertical inheritance of nhe was observed, and since nhe-negative strains appear to be extremely rare, we suggest that fitness loss may be associated with deletion or horizontal transfer of the nhe operon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0529-4) contains supplementary material, which is available to authorized users.

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.

The Water Cycle, a Potential Source of the Bacterial Pathogen Bacillus cereus

The behaviour of the sporulating soil-dwelling Bacillus cereus sensu lato (B. cereus sl) which includes foodborne pathogenic strains has been extensively studied in relation to its various animal hosts. The aim of this environmental study was to investigate the water compartments (rain and soil water, as well as groundwater) closely linked to the primary B. cereus sl reservoir, for which available data are limited. B. cereus sl was present, primarily as spores, in all of the tested compartments of an agricultural site, including water from rain to groundwater through soil. During rain events, leachates collected after transfer through the soil eventually reached the groundwater and were loaded with B. cereus sl. In groundwater samples, newly introduced spores of a B. cereus model strain were able to germinate, and vegetative cells arising from this event were detected for up to 50 days. This first B. cereus sl investigation in the various types of interrelated environments suggests that the consideration of the aquatic compartment linked to soil and to climatic events should provide a better understanding of B. cereus sl ecology and thus be relevant for a more accurate risk assessment of food poisoning caused by B. cereus sl pathogenic strains.

Characteristics and phylogeny of Bacillus cereus strains isolated from Maari, a traditional West African food condiment

Maari is a spontaneously fermented food condiment made from baobab tree seeds in West African countries. This type of product is considered to be safe, being consumed by millions of people on a daily basis. However, due to the spontaneous nature of the fermentation the human pathogen Bacillus cereus occasionally occurs in Maari. This study characterizes succession patterns and pathogenic potential of B. cereus isolated from the raw materials (ash, water from a drilled well (DW) and potash), seed mash throughout fermentation (0-96h), after steam cooking and sun drying (final product) from two production sites of Maari. Aerobic mesophilic bacterial (AMB) counts in raw materials were of 10(5)cfu/ml in DW, and ranged between 6.5×10(3) and 1.2×10(4)cfu/g in potash, 10(9)-10(10)cfu/g in seed mash during fermentation and 10(7) - 10(9) after sun drying. Fifty three out of total 290 AMB isolates were identified as B. cereus sensu lato by use of ITS-PCR and grouped into 3 groups using PCR fingerprinting based on Escherichia coli phage-M13 primer (M13-PCR). As determined by panC gene sequencing, the isolates of B. cereus belonged to PanC types III and IV with potential for high cytotoxicity. Phylogenetic analysis of concatenated sequences of glpF, gmk, ilvD, pta, pur, pycA and tpi revealed that the M13-PCR group 1 isolates were related to B. cereus biovar anthracis CI, while the M13-PCR group 2 isolates were identical to cereulide (emetic toxin) producing B. cereus strains. The M13-PCR group 1 isolates harboured poly-γ-D-glutamic acid capsule biosynthesis genes capA, capB and capC showing 99-100% identity with the environmental B. cereus isolate 03BB108. Presence of cesB of the cereulide synthetase gene cluster was confirmed by PCR in M13-PCR group 2 isolates. The B. cereus harbouring the cap genes were found in potash, DW, cooking water and at 8h fermentation. The "emetic" type B. cereus were present in DW, the seed mash at 48-72h of fermentation and in the final product, while the remaining isolates (PanC type IV) were detected in ash, at 48-72h fermentation and in the final product. This work sheds light on the succession and pathogenic potential of B. cereus species in traditional West African food condiment and clarifies their phylogenetic relatedness to B. cereus biovar anthracis. Future implementation of GMP and HACCP and development of starter cultures for controlled Maari fermentations will help to ensure a safe product.

Division of labour and terminal differentiation in a novel Bacillus thuringiensis strain

A major challenge in bacterial developmental biology has been to understand the mechanisms underlying cell fate decisions. Some differentiated cell types display cooperative behaviour. Cooperation is one of the greatest mysteries of evolutionary biology and microbes have been considered as an excellent system for experimentally testing evolution theories. Bacillus thuringiensis (Bt) is a spore-forming bacterium, which is genetically closely related to B. anthracis, the agent of anthrax, and to B. cereus, an opportunistic human pathogen. The defining feature that distinguishes Bt from its relatives is its ability to produce crystal inclusions in the sporulating cells. These toxins are solubilized after ingestion and are cooperative public goods in insect hosts. In this study, we describe a Bt strain LM1212 that presents the unique ability to terminally differentiate into crystal producers and spore formers. Transcriptional analysis based on lacZ and gfp reporter genes suggested that this phenotype is the consequence of a new type of cell differentiation associated with a novel regulation mode of cry gene expression. The differentiating crystal-producer phenotype has higher spore productivity than a typical Bt strain and is better able to compete with Cry toxin null 'cheaters'. Potentially, this division of labour provides additional fitness benefits in terms of spore viability or durability of Cry toxin.

Crystalliferous Bacillus cereus group bacteria from a Maryland hardwood forest are dominated by psychrotolerant strains

Crystal-forming bacteria of the Bacillus cereus group were isolated from soil samples collected at different elevations within a mixed hardwood forest in central Maryland, and their phylogenetic relationships determined by multilocus sequence analysis. The vast majority of isolates obtained were associated with two phylogenetic groups known to be psychrotolerant, with very few isolates representing phylogenetic groups more typically associated with Bacillus thuringiensis. Isolates from the psychrotolerant groups were found to grow on solid media at 7°C. Isolates of 11 highly related, novel sequence types (STs) from the psychrotolerant group that includes Bacillus weihenstephanensis were generally found at higher elevations, and were not associated with soils near streams. Isolates of two related STs from the second psychrotolerant group were nearly always found at the bottoms of ravines near streams, in areas abundant in earthworm castings.

Eco-genetic structure of Bacillus cereus sensu lato populations from different environments in northeastern Poland

The Bacillus cereus group, which includes entomopathogens and etiologic agents of foodborne illness or anthrax, persists in various environments. The basis of their ecological diversification remains largely undescribed. Here we present the genetic structure and phylogeny of 273 soil B. cereus s.l. isolates from diverse habitats in northeastern Poland, with samplings acquired from the last European natural forest (Białowieża National Park), the largest marshes in Europe (Biebrza National Park), and a farm. In multi-locus sequence typing (MLST), despite negative selection in seven housekeeping loci, the isolates exhibited high genetic diversity (325 alleles), mostly resulting from mutation events, and represented 148 sequencing types (131 STs new and 17 STs already described) grouped into 19 complexes corresponding with bacterial clones, and 80 singletons. Phylogenetic analyses showed that 74% of the isolates clustered with B. cereus s.l. environmental references (clade III), while only 11 and 15%, respectively, grouped with isolates of clinical origin (clade I), and B. cereus ATCC 14579 and reference B. thuringiensis (clade II). Predominantly within clade III, we found lineages adapted to low temperature (thermal ecotypes), while putative toxigenic isolates (cytK-positive) were scattered in all clades of the marsh and farm samplings. The occurrence of 92% of STs in bacilli originating from one habitat, and the description of new STs for 78% of the isolates, strongly indicate the existence of specific genotypes within the natural B. cereus s.l. populations. In contrast to the human-associated B. cereus s.l. that exhibit a significant level of similarity, the environmental isolates appear more complex. Thus we propose dividing B. cereus s.l. into two groups, the first including environmental isolates, and the second covering those that are of clinical relevance.

Phylogenetic Analysis ofBacillus cereusIsolates from Severe Systemic Infections Using Multilocus Sequence Typing Scheme

Bacillus cereus strains from cases of severe or lethal systemic infections, including respiratory symptoms cases, were analyzed using multilocus sequence typing scheme of B. cereus MLST database. The isolates were evenly distributed between the two main clades, and 60% of them had allele profiles new to the database. Half of the collection's strains clustered in a lineage neighboring Bacillus anthracis phylogenetic origin. Strains from lethal cases with respiratory symptoms were allocated in both main clades. This is the first report of strains causing respiratory symptoms to be identified as genetically distant from B. anthracis. The phylogenetic location of the presented here strains was compared with all previously submitted to the database isolates from systemic infections, and were found to appear in the same clusters where clinical isolates from other studies had been assigned. It seems that the pathogenic strains are forming clusters on the phylogenetic tree.

Extended genetic analysis of Brazilian isolates of Bacillus cereus and Bacillus thuringiensis

Multiple locus sequence typing (MLST) was undertaken to extend the genetic characterization of 29 isolates of Bacillus cereus and Bacillus thuringiensis previously characterized in terms of presence/absence of sequences encoding virulence factors and via variable number tandem repeat (VNTR). Additional analysis involved polymerase chain reaction for the presence of sequences (be, cytK, inA, pag, lef, cya and cap), encoding putative virulence factors, not investigated in the earlier study. MLST analysis ascribed novel and unique sequence types to each of the isolates. A phylogenetic tree was constructed from a single sequence of 2,838 bp of concatenated loci sequences. The strains were not monophyletic by analysis of any specific housekeeping gene or virulence characteristic. No clear association in relation to source of isolation or to genotypic profile based on the presence or absence of putative virulence genes could be identified. Comparison of VNTR profiling with MLST data suggested a correlation between these two methods of genetic analysis. In common with the majority of previous studies, MLST was unable to provide clarification of the basis for pathogenicity among members of the B. cereus complex. Nevertheless, our application of MLST served to reinforce the notion that B. cereus and B. thuringiensis should be considered as the same species.

Multilocus sequence analysis of nectar pseudomonads reveals high genetic diversity and contrasting recombination patterns

The genetic and evolutionary relationships among floral nectar-dwelling Pseudomonas 'sensu stricto' isolates associated to South African and Mediterranean plants were investigated by multilocus sequence analysis (MLSA) of four core housekeeping genes (rrs, gyrB, rpoB and rpoD). A total of 35 different sequence types were found for the 38 nectar bacterial isolates characterised. Phylogenetic analyses resulted in the identification of three main clades [nectar groups (NGs) 1, 2 and 3] of nectar pseudomonads, which were closely related to five intrageneric groups: Pseudomonas oryzihabitans (NG 1); P. fluorescens, P. lutea and P. syringae (NG 2); and P. rhizosphaerae (NG 3). Linkage disequilibrium analysis pointed to a mostly clonal population structure, even when the analysis was restricted to isolates from the same floristic region or belonging to the same NG. Nevertheless, signatures of recombination were observed for NG 3, which exclusively included isolates retrieved from the floral nectar of insect-pollinated Mediterranean plants. In contrast, the other two NGs comprised both South African and Mediterranean isolates. Analyses relating diversification to floristic region and pollinator type revealed that there has been more unique evolution of the nectar pseudomonads within the Mediterranean region than would be expected by chance. This is the first work analysing the sequence of multiple loci to reveal geno- and ecotypes of nectar bacteria.

Differentiation of strains from the Bacillus cereus group by RFLP-PFGE genomic fingerprinting

Bacillus mycoides, Bacillus pseudomycoides, Bacillus weihenstephanensis, Bacillus anthracis, Bacillus thuringiensis, and Bacillus cereus belong to the B. cereus group. The last three species are characterized by different phenotype features and pathogenicity spectrum, but it has been shown that these species are genetically closely related. The macrorestriction analysis of the genomic DNA with the NotI enzyme was used to generate polymorphism of restriction profiles for 39 food-borne isolates (B. cereus, B. mycoides) and seven reference strains (B. mycoides, B. thuringiensis, B. weihenstephanensis, and B. cereus). The PFGE method was applied to differentiate the examined strains of the B. cereus group. On the basis of the unweighted pair group method with the arithmetic mean method and Dice coefficient, the strains were divided into five clusters (types A-E), and the most numerous group was group A (25 strains). A total of 21 distinct pulsotypes were observed. The RFLP-PFGE analysis was successfully used for the differentiation and characterization of B. cereus and B. mycoides strains isolated from different food products.