Multiple-locus sequence typing and analysis of toxin genes in Bacillus cereus food-borne isolates

Genetic Characterization, Antibiotic Resistance, and Virulence Genes Profiling of Bacillus cereus Strains from Various Foods in Japan

Bacillus cereus sensu stricto is a foodborne pathogen that causes food poisoning. Their spore and biofilm-forming abilities persist in various environments and foods. This study investigated the prevalence, virulence, antibiotic resistance, and genetic diversity of B. cereus s. s. strains isolated from various food samples. Of 179 samples, 22.34% were positive for B. cereus s. s., with significantly high detection rates in milk products and raw chicken meat. Forty strains were isolated from positive samples. Matrix-assisted laser desorption ionization/time of flight mass spectrometry analysis revealed nine distinct clusters and multi-locus sequence typing revealed 34 sequence types including 23 novel sequences, demonstrating high genetic diversity among the isolates. PCR analysis revealed that all the strains contained at least one toxin gene, but none contained the cytK gene. Antibiotic resistance tests revealed that all isolates were classified as multidrug-resistant, with high resistance levels, particularly to β-lactam antibiotics and vancomycin, but were susceptible to gentamicin. All isolates showed variations in biofilm formation. This study highlights the significant public health risk due to B. cereus s. s. and underscores the need for stringent monitoring and control measures in food production to manage antimicrobial resistance and ensure food safety.

High Genetic Diversity and Virulence Potential in Bacillus cereus sensu lato Isolated from Milk and Cheeses in Apulia Region, Southern Italy

The Bacillus cereus group includes species that act as food-borne pathogens causing diarrheal and emetic symptoms. They are widely distributed and can be found in various foods. In this study, out of 550 samples of milk and cheeses, 139 (25.3%) were found to be contaminated by B. cereus sensu lato (s.l.). One isolate per positive sample was characterized by Multilocus Sequence Typing (MLST) and for the presence of ten virulence genes. Based on MLST, all isolates were classified into 73 different sequence types (STs), of which 12 isolates were assigned to new STs. Virulence genes detection revealed that 90% and 61% of the isolates harboured the nheABC and the hblCDA gene cluster, respectively. Ninety-four percent of the isolates harboured the enterotoxin genes entS and entFM; 8% of the isolates possessed the ces gene. Thirty-eight different genetic profiles were identified, suggesting a high genetic diversity. Our study clearly shows the widespread diffusion of potentially toxigenic isolates of B. cereus s.l. in milk and cheeses in the Apulia region highlighting the need to adopt GMP and HACCP procedures along every step of the milk and cheese production chain in order to reduce the public health risk linked to the consumption of foods contaminated by B. cereus s.l.

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.

Phylogenetic and protein prediction analysis reveals the taxonomically diverse distribution of virulence factors in Bacillus cereus strains

Bacillus cereus is a food contaminant with widely varying enterotoxic potential due to its virulence proteins. In this article, phylogenetic analysis of the amino acid sequences from the whole-genomes of 41 strains, evolutionary distance calculation of the amino acid sequences of the virulence genes, and functional and structural predictions of the virulence proteins were performed to reveal the taxonomically diverse distribution of virulence factors. The genome evolution of the strains showed a clustering trend based on the protein-coding virulence genes. The strains of B. cereus have evolved into non-toxic risk and toxic risk clusters with medium-high- and medium-low-risk subclusters. The evolutionary transfer distances of incomplete virulence genes relative to housekeeping genes were greater than those of complete virulence genes, and the distance values of HblACD were higher than those of nheABC and CytK among the complete virulence genes. Cytoplasmic localization was impossible for all the virulence proteins, and NheB, NheC, Hbl-B, and Hbl-L₁ were predicted to be extracellular. Nhe and Hbl proteins except CytK had similar spatial structures. The predicted structures of Nhe and Hbl mainly showed ‘head’ and ‘tail’ domains. The ‘head’ of NheA and Hbl-B, including two α-helices separated by β-tongue strands, might play a special role in the formation of Nhe trimers and Hbl trimers, respectively. The ‘cap’ of CytK, which includes two ‘latches’ with many β-sheets, formed a β-barrel structure with pores, and a ‘rim’ balanced the structure. The evolution of B. cereus strains showed a clustering tendency based on the protein-coding virulence genes, and the complete virulence-gene operon combination had higher relative genetic stability. The beta-tongue or latch associated with β-sheet folding might play an important role in the binding of virulence structures and pore-forming toxins in B. cereus.

Characterization of Vegetative Bacillus cereus and Bacillus subtilis Strains Isolated from Processed Cheese Products in an Italian Dairy Plant

Processed cheese is a commercial product characterized by high microbiological stability and extended shelf life obtained through the application of severe heat treatment. However, spore-forming bacteria can survive through thermal processes. Among them, microorganisms belonging to Bacillus genus have been reported. In this study, we examined the microbiological population of the first hours' production of processed cheeses in an Italian dairy plant during two seasons, between June and October 2020. Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) was used to identify bacteria colonies, allowing the isolation of Bacillus cereus and Bacillussubtilis strains. These results were further confirmed by amplification and sequencing of 16 rRNA bacterial region. A multi-locus sequence type (MLST) analysis was performed to assess the genetic similarity among a selection of isolates. The fourteen B. cereus strains showed two sequence types: ST-32 was observed in only one strain and the ST-371 in the remaining thirteen isolates. On the contrary, all twenty-one B. subtlis strains, included in the study, showed a new allelic profile for the pycA gene, resulting in a new sequence type: ST-249. For B. cereus strains, analysis of toxin genes was performed. All isolates were positive for nheABC, entFM, and cytK, while hblABCD, bceT, and ces were not detected. Moreover, the biofilm-forming ability of B. cereus and B. subtilis strains was assessed, and all selected isolates proved to be biofilm formers (most of them were stronger producers). Considering the genetical similarity between isolates, jointly with the capacity to produce biofilm, the presence of a recurring Bacillus population could be hypothesized.

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.

Keeping up with the Bacillus cereus group: taxonomy through the genomics era and beyond

The Bacillus cereus group, also known as B. cereus sensu lato (s.l.), is a species complex that contains numerous closely related lineages, which vary in their ability to cause illness in humans and animals. The classification of B. cereus s.l. isolates into species-level taxonomic units is thus essential for informing public health and food safety efforts. However, taxonomic classification of these organisms is challenging. Numerous-often conflicting-taxonomic changes to the group have been proposed over the past two decades, making it difficult to remain up to date. In this review, we discuss the major nomenclatural changes that have accumulated in the B. cereus s.l. taxonomic space prior to 2020, particularly in the genomic sequencing era, and outline the resulting problems. We discuss several contemporary taxonomic frameworks as applied to B. cereus s.l., including (i) phenotypic, (ii) genomic, and (iii) hybrid nomenclatural frameworks, and we discuss the advantages and disadvantages of each. We offer suggestions as to how readers can avoid B. cereus s.l. taxonomic ambiguities, regardless of the nomenclatural framework(s) they choose to employ. Finally, we discuss future directions and open problems in the B. cereus s.l. taxonomic realm, including those that cannot be solved by genomic approaches alone.

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.

Multi-locus sequence typing and virulence profile in Bacillus cereus sensu lato strains isolated from dairy products

Members of Bacillus cereus group are important food contaminants and they are of relevant interest in food safety and public heath due to their ability to cause two distinct forms of food poisoning, emetic and diarrhoeal syndrome. In the present study, 90 strains of B. cereus isolated from dairy products, have been typed using Multilocus Sequence Typing (MLST) analysis and investigated for the occurrence of 10 enterotoxigenic genes (hblA, hblC, hblD, nheA, nheB, nheC, cytK, entFM, entS and bceT) and one emetogenic gene (ces), to determine their genetic diversity. A total of 58 sequence types were identified and among these 17 were signalled as new profiles. Among the virulence genes, the majority of our strains carried the entS (92%), entFM (86%), nhe (82%) and cytK (72%) genes. All remaining genes were identified in at least one strain with different prevalence, stressing the genetic diversity, how even the different grade of pathogenicity of B. cereus isolated from dairy products.

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.

No Assembly Required: Using BTyper3 to Assess the Congruency of a Proposed Taxonomic Framework for the Bacillus cereus Group With Historical Typing Methods

The Bacillus cereus group, also known as B. cereus sensu lato (s.l.), is a species complex comprising numerous closely related lineages, which vary in their ability to cause illness in humans and animals. The classification of B. cereus s.l. isolates into species-level taxonomic units is essential for facilitating communication between and among microbiologists, clinicians, public health officials, and industry professionals, but is not always straightforward. A recently proposed genomospecies-subspecies-biovar taxonomic framework aims to provide a standardized nomenclature for this species complex but relies heavily on whole-genome sequencing (WGS). It thus is unclear whether popular, low-cost typing methods (e.g., single- and multi-locus sequence typing) remain congruent with the proposed taxonomy. Here, we characterize 2,231 B. cereus s.l. genomes using a combination of in silico (i) average-nucleotide identity (ANI)-based genomospecies assignment, (ii) ANI-based subspecies assignment, (iii) seven-gene multi-locus sequence typing (MLST), (iv) single-locus panC group assignment, (v) rpoB allelic typing, and (vi) virulence factor detection. We show that sequence types (STs) assigned using MLST can be used for genomospecies assignment, and we provide a comprehensive list of ST/genomospecies associations. For panC group assignment, we show that an adjusted, eight-group framework is largely, albeit not perfectly, congruent with the proposed eight-genomospecies taxonomy, as panC alone may not distinguish (i) B. luti from Group II B. mosaicus and (ii) B. paramycoides from Group VI B. mycoides. We additionally provide a list of loci that capture the topology of the whole-genome B. cereus s.l. phylogeny that may be used in future sequence typing efforts. For researchers with access to WGS, MLST, and/or panC data, we showcase how our recently released software, BTyper3 (https://github.com/lmc297/BTyper3), can be used to assign B. cereus s.l. isolates to taxonomic units within this proposed framework with little-to-no user intervention or domain-specific knowledge of B. cereus s.l. taxonomy. We additionally outline a novel method for assigning B. cereus s.l. genomes to pseudo-gene flow units within proposed genomospecies. The results presented here highlight the backward-compatibility and accessibility of the recently proposed genomospecies-subspecies-biovar taxonomic framework and illustrate that WGS is not a necessity for microbiologists who want to use the proposed nomenclature effectively.

A Study on Prevalence and Characterization of Bacillus cereus in Ready-to-Eat Foods in China

Bacillus cereus is widely distributed in different food products and can cause a variety of symptoms associated with food poisoning. Since ready-to-eat (RTE) foods are not commonly sterilized by heat treatment before consumption, B. cereus contamination may cause severe food safety problems. In this study, we investigated the prevalence of B. cereus in RTE food samples from different regions of China and evaluated the levels of bacterial contamination, antibiotic resistance, virulence gene distribution, and genetic polymorphisms of these isolates. Of the tested retail RTE foods, 35% were positive for B. cereus, with 39 and 83% of the isolated strains harboring the enterotoxin-encoding hblACD and nheABC gene clusters, respectively. The entFM gene was detected in all B. cereus strains. The cytK gene was present in 68% of isolates, but only 7% harbored the emetic toxin-encoding gene cesB. Antimicrobial susceptibility testing revealed that the majority of the isolates were resistant not only to most β-lactam antibiotics, but also to rifamycin. Multilocus sequence typing (MLST) revealed that the 368 isolates belonged to 192 different sequence types (STs) including 93 new STs, the most prevalent of which was ST26. Collectively, our study indicates the prevalence, bacterial contamination levels, and biological characteristics of B. cereus isolated from RTE foods in China and demonstrates the potential hazards of B. cereus in RTE foods.

Bacillus cereus-Attributable Primary Cutaneous Anthrax-Like Infection in Newborn Infants, India

During March 13-June 23, 2018, anthrax-like cutaneous lesions attributed to the Bacillus cereus group of organisms developed in 12 newborns in India. We traced the source of infection to the healthcare kits used for newborn care. We used multilocus sequence typing to characterize the 19 selected strains from various sources in hospital settings, including the healthcare kits. This analysis revealed the existence of a genetically diverse population comprising mostly new sequence types. Phylogenetic analysis clustered most strains into the previously defined clade I, composed primarily of pathogenic bacilli. We suggest that the synergistic interaction of nonhemolytic enterotoxin and sphingomyelinase might have a role in the development of cutaneous lesions. The infection was controlled by removing the healthcare kits and by implementing an ideal housekeeping program. All the newborns recovered after treatment with ciprofloxacin and amikacin.

Comparative genomic survey of Bacillus cereus sensu stricto isolates from the dairy production chain in Brazil

The genomes of 262 Bacillus cereus isolates were analyzed including 69 isolates sampled from equipment, raw milk and dairy products from Brazil. The population structure of isolates showed strains belonging to known phylogenetic groups II, III, IV, V and VI. Almost all the isolates obtained from dairy products belonged to group III. Investigation of specific alleles revealed high numbers of isolates carrying toxin-associated genes including cytK (53.62%), hblA (59.42%), hblC (44.93%), hblD (53.62%), nheA (84.06%), nheB (89.86%) and nheC (84.06%) with isolates belonging to groups IV and V having significant higher prevalence of hblACD and group IV of CytK genes. Strains from dairy products had significantly lower prevalence of CytK and hblACD genes compared to isolates from equipment and raw milk/bulk tanks. Genes related to sucrose metabolism were detected at higher frequency in isolates obtained from raw milk compared to strains from equipment and utensils. The population genomic analysis demonstrated the diversity of strains and variability of putative function among B. cereus group isolates in Brazilian dairy production, with large numbers of strains potentially able to cause foodborne illness. This detailed information will contribute to targeted interventions to reduce milk contamination and spoilage associated with B. cereus in Brazil.

Bacillus cereus Isolated From Vegetables in China: Incidence, Genetic Diversity, Virulence Genes, and Antimicrobial Resistance

Bacillus cereus is a food-borne opportunistic pathogen that can induce diarrheal and emetic symptoms. It is widely distributed in different environments and can be found in various foods, including fresh vegetables. As their popularity grows worldwide, the risk of bacterial contamination in fresh vegetables should be fully evaluated, particularly in vegetables that are consumed raw or processed minimally, which are not commonly sterilized by enough heat treatment. Thereby, it is necessary to perform potential risk evaluation of B. cereus in vegetables. In this study, 294 B. cereus strains were isolated from vegetables in different cities in China to analyze incidence, genetic polymorphism, presence of virulence genes, and antimicrobial resistance. B. cereus was detected in 50% of all the samples, and 21/211 (9.95%) of all the samples had contamination levels of more than 1,100 MPN/g. Virulence gene detection revealed that 95 and 82% of the isolates harbored nheABC and hblACD gene clusters, respectively. Additionally, 87% of the isolates harbored cytK gene, and 3% of the isolates possessed cesB. Most strains were resistant to rifampicin and β-lactam antimicrobials but were sensitive to imipenem, gentamicin, ciprofloxacin, kanamycin, telithromycin, ciprofloxacin, and chloramphenicol. In addition, more than 95.6% of the isolates displayed resistance to three kinds of antibiotics. Based on multilocus sequence typing, all strains were classified into 210 different sequence types (STs), of which 145 isolates were assigned to 137 new STs. The most prevalent ST was ST770, but it included only eight isolates. Taken together, our research provides the first reference for the incidence and characteristics of B. cereus in vegetables collected throughout China, indicating a potential hazard of B. cereus when consuming vegetables without proper handling.

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.

RNA-Seq Analysis of Antibiotic-Producing Bacillus subtilis SC-8 Reveals a Role for Small Peptides in Controlling PapR Signaling

Bacillus subtilis SC-8 (BSSC8) shows a narrow antimicrobial activity against the Bacillus cereus group. Previously, B. cereus-derived PapR as a signal peptide to stimulate PlcR, which plays a significant role in regulating the transcription of virulence factors, was assumed to stimulate antibiotic production in BSSC8. To better understand the functional role of PapR in the antibiotic production of BSSC8 and the interspecies interaction, the global transcriptomic profiling of BSSC8 was investigated using RNA-Seq in this study. Small peptides derived from B. cereus wild type (WTBC) and a papR-deleted mutant strain (MTBC) were individually supplied to BSSC8 cultures, and changes in global transcription levels were compared by RNA-Seq. In the presence of WTBC small peptides, more genes (80.9%) were significantly upregulated than in cells exposed to MTBC small peptides. Specifically, 48.8 and 83.4% of genes involved in glycolysis and the TCA cycle, respectively, showed changes in transcription levels in response to small peptides from both strains. Of the genes showing the alterations, 35.0% (glycolysis) and 60.0% (TCA cycle) of transcripts were significantly regulated only in response to WTBC-derived small peptides. Furthermore, the expression of biosynthetic genes encoding several known antibiotics in BSSC8 was further decreased in response to WTBC small peptides.

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.

Comparative Genomics of Bacillus thuringiensis Reveals a Path to Specialized Exploitation of Multiple Invertebrate Hosts

Understanding the genetic basis of host shifts is a key genomic question for pathogen and parasite biology. The Bacillus cereus group, which encompasses Bacillus thuringiensis and Bacillus anthracis, contains pathogens that can infect insects, nematodes, and vertebrates. Since the target range of the essential virulence factors (Cry toxins) and many isolates is well known, this group presents a powerful system for investigating how pathogens can diversify and adapt to phylogenetically distant hosts. Specialization to exploit insects occurs at the level of the major clade and is associated with substantial changes in the core genome, and host switching between insect orders has occurred repeatedly within subclades. The transfer of plasmids with linked cry genes may account for much of the adaptation to particular insect orders, and network analysis implies that host specialization has produced strong associations between key toxin genes with similar targets. Analysis of the distribution of plasmid minireplicons shows that plasmids with orf156 and orf157, which carry genes encoding toxins against Lepidoptera or Diptera, were contained only by B. thuringiensis in the specialized insect clade (clade 2), indicating that tight genome/plasmid associations have been important in adaptation to invertebrate hosts. Moreover, the accumulation of multiple virulence factors on transposable elements suggests that cotransfer of diverse virulence factors is advantageous in terms of expanding the insecticidal spectrum, overcoming insect resistance, or through gains in pathogenicity via synergistic interactions between toxins.IMPORTANCE Population genomics have provided many new insights into the formation, evolution, and dynamics of bacterial pathogens of humans and other higher animals, but these pathogens usually have very narrow host ranges. As a pathogen of insects and nematodes, Bacillus thuringiensis, which produces toxins showing toxicity to many orders of insects and other invertebrates, can be used as a model to study the evolution of pathogens with wide host ranges. Phylogenomic analysis revealed that host specialization and switching occur at the level of the major clade and subclade, respectively. A toxin gene co-occurrence network indicates that multiple toxins with similar targets were accumulated by the same cell in the whole species. This accumulation may be one of the strategies that B. thuringiensis has used to fight against host resistance. This kind of formation and evolution of pathogens represents a different path used against multiple invertebrate hosts from that used against higher animals.

Surviving Between Hosts: Sporulation and Transmission

To survive adverse conditions, some bacterial species are capable of developing into a cell type, the "spore," which exhibits minimal metabolic activity and remains viable in the presence of multiple environmental challenges. For some pathogenic bacteria, this developmental state serves as a means of survival during transmission from one host to another. Spores are the highly infectious form of these bacteria. Upon entrance into a host, specific signals facilitate germination into metabolically active replicating organisms, resulting in disease pathogenesis. In this article, we will review spore structure and function in well-studied pathogens of two genera, Bacillus and Clostridium, focusing on Bacillus anthracis and Clostridium difficile, and explore current data regarding the lifestyles of these bacteria outside the host and transmission from one host to another.

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.

Prevalence and Characterization of Campylobacter jejuni Isolated from Retail Chicken in Tianjin, China

Campylobacter jejuni is an important foodborne pathogen worldwide; however, there is a lack of information on the prevalence and antibiotic-resistant profile of C. jejuni in the People's Republic of China. We determined the prevalence and characteristics of C. jejuni on the retail level in Tianjin, one of the five national central cities in China. A total of 227 samples of chicken wings, legs, and breasts were collected from supermarkets and wet markets; 42 of these samples were confirmed to be positive for Campylobacter contamination. The contamination rates of C. jejuni and other Campylobacter species were 13.7% (31 of 227 samples) and 5.7% (13 of 227 samples), respectively. A group of 31 C. jejuni isolates was subjected to antimicrobial susceptibility testing. All (100%) the selected isolates were resistant to ciprofloxacin and nalidixic acid; 77.4% were resistant to tetracycline, 67.7% to doxycycline, 35.5% to gentamicin, 25.8% to clindamycin and florfenicol, 19.4% to chloramphenicol, and 12.9% to erythromycin and azithromycin. A remarkably high proportion (41.9%) of multidrug-resistant isolates was identified. Multilocus sequence typing was conducted to study the population structure of the C. jejuni strains and their relationship to human isolates. The correlation between antimicrobial resistance traits and certain sequence types (STs) or clonal complexes was determined as well. A great genetic diversity of poultry isolates was identified, with 11 STs belonging to 6 clonal complexes and 11 singleton STs. The novel STs accounted for 40.9% (n = 9) of the 22 STs. ST-21, ST-353, ST-354, ST-443, ST-607, and ST-828 complexes had been previously identified from human isolates. This study revealed an extensive level of antimicrobial resistance and genetic diversity in C. jejuni isolated from chicken products in Tianjin, highlighting the necessity of performing enforced interventions to reduce Campylobacter prevalence in China.

Genetic diversity and population structure of the Bacillus cereus group bacteria from diverse marine environments

The phylogenetic diversity of marine bacteria belonged to the Bacillus cereus group has not been well investigated. Here, we present the genetic diversity and population structure of 71 bacteria from diverse marine environments, using a multilocus sequence typing (MLST) approach and the analyses of digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) based on some representative genomic sequences. The MLST analysis demonstrated that these isolates were highly diverse and a wide distribution in marine environments and some of them showed niche specificity to some extent. They were assigned to 27 sequence types (STs) with 23 novel STs. Phylogenetic analysis of 82 bacteria containing 11 type strains based on MLST discriminated them as 20 clusters including 10 new ones. Both the dDDH and ANI results supported the proposition that each of 20 clusters represented one independent species, including 10 putative novel species. Values of 98.3% of MLST similarity and 96.2% of ANI were proposed as the standard for the species definition of this group. In summary, the first insight into the phylogenetic diversity of the group bacteria from marine environments will contribute to better understanding of their ecological role and evolution in contrast with terrestrial environments.

Production of hemolysin BL by Bacillus cereus group isolates of dairy origin is associated with whole-genome phylogenetic clade

BACKGROUND

Bacillus cereus group isolates that produce diarrheal or emetic toxins are frequently isolated from raw milk and, in spore form, can survive pasteurization. Several species within the B. cereus group are closely related and cannot be reliably differentiated by established taxonomical criteria. While B. cereus is traditionally recognized as the principal causative agent of foodborne disease in this group, there is a need to better understand the distribution and expression of different toxin and virulence genes among B. cereus group food isolates to facilitate reliable characterization that allows for assessment of the likelihood of a given isolate to cause a foodborne disease.

RESULTS

We performed whole genome sequencing of 22 B. cereus group dairy isolates, which represented considerable genetic diversity not covered by other isolates characterized to date. Maximum likelihood analysis of these genomes along with 47 reference genomes representing eight validly published species revealed nine phylogenetic clades. Three of these clades were represented by a single species (B. toyonensis -clade V, B. weihenstephanensis - clade VI, B. cytotoxicus - VII), one by two dairy-associated isolates (clade II; representing a putative new species), one by two species (B. mycoides, B. pseudomycoides - clade I) and four by three species (B. cereus, B. thuringiensis, B. anthracis - clades III-a, b, c and IV). Homologues of genes encoding a principal diarrheal enterotoxin (hemolysin BL) were distributed across all, except the B. cytotoxicus clade. Using a lateral flow immunoassay, hemolysin BL was detected in 13 out of 18 isolates that carried hblACD genes. Isolates from clade III-c (which included B. cereus and B. thuringiensis) consistently did not carry hblACD and did not produce hemolysin BL. Isolates from clade IV (B. cereus, B. thuringiensis) consistently carried hblACD and produced hemolysin BL. Compared to others, clade IV was significantly (p = 0.0001) more likely to produce this toxin. Isolates from clade VI (B. weihenstephanensis) carried hblACD homologues, but did not produce hemolysin BL, possibly due to amino acid substitutions in different toxin-encoding genes.

CONCLUSIONS

Our results demonstrate that production of diarrheal enterotoxin hemolysin BL is neither inclusive nor exclusive to B. cereus sensu stricto, and that phylogenetic classification of isolates may be better than taxonomic identification for assessment of B. cereus group isolates risk for causing a diarrheal foodborne disease.

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.

Microbial dynamics during shelf-life of industrial Ricotta cheese and identification of a Bacillus strain as a cause of a pink discolouration

Dairy products are perishable and have to be preserved from spoilage during the food chain to achieve the desired shelf-life. Ricotta is a typical Italian soft dairy food produced by heat coagulation of whey proteins and is considered to be a light and healthy product. The shelf-life of Ricotta could be extended, as required by the international food trade market; however, heat resistant microflora causes spoilage and poses issues regarding the safety of the product. Next-generation sequencing (NGS) applied to the Ricotta samples defined the composition of the microbial community in-depth during the shelf-life. The analysis demonstrated the predominance of spore-forming bacteria throughout the shelf-life, mostly belonging to Bacillus, Paenibacillus and Clostridium genera. A strain involved in spoilage and causing a pink discolouration of Ricotta was isolated and characterised as Bacillus mycoides/weihenstephanensis. This is the first report of a food discolouration caused by a toxigenic strain belonging to the Bacillus cereus group that resulted the predominant strain in the community of the defective ricotta. These results suggest that the processing of raw materials to eliminate spores and residual microflora could be essential for improving the quality and the safety of the product and to extend the shelf-life of industrial Ricotta.

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.

A Genomic View of Lactobacilli and Pediococci Demonstrates that Phylogeny Matches Ecology and Physiology

Lactobacilli are used widely in food, feed, and health applications. The taxonomy of the genus Lactobacillus, however, is confounded by the apparent lack of physiological markers for phylogenetic groups of lactobacilli and the unclear relationships between the diverse phylogenetic groups. This study used the core and pan-genomes of 174 type strains of Lactobacillus and Pediococcus to establish phylogenetic relationships and to identify metabolic properties differentiating phylogenetic groups. The core genome phylogenetic tree separated homofermentative lactobacilli and pediococci from heterofermentative lactobacilli. Aldolase and phosphofructokinase were generally present in homofermentative but not in heterofermentative lactobacilli; a two-domain alcohol dehydrogenase and mannitol dehydrogenase were present in most heterofermentative lactobacilli but absent in most homofermentative organisms. Other genes were predominantly present in homofermentative lactobacilli (pyruvate formate lyase) or heterofermentative lactobacilli (lactaldehyde dehydrogenase and glycerol dehydratase). Cluster analysis of the phylogenomic tree and the average nucleotide identity grouped the genus Lactobacillus sensu lato into 24 phylogenetic groups, including pediococci, with stable intra- and intergroup relationships. Individual groups may be differentiated by characteristic metabolic properties. The link between phylogeny and physiology that is proposed in this study facilitates future studies on the ecology, physiology, and industrial applications of lactobacilli.

Characterization of the spore-forming Bacillus cereus sensu lato group and Clostridium perfringens bacteria isolated from the Australian dairy farm environment

Background

The Bacillus cereus sensu lato group and Clostridium perfringens are spore-forming bacteria often associated with food spoilage and which can cause emetic and diarrheal syndromes in humans and ruminants. This study characterised the phenotypes and genotypes of 50 Bacillus cereus s. l. isolates and 26 Clostridium perfringens isolates from dairy farms environments in Victoria, Australia.

Results

Five of the seven B. cereus s. l. species were isolated, and analysis of the population diversity using Pulsed-Field Gel Electrophoresis (PFGE) suggested that the populations are largely distinct to each farm. Enterotoxin production by representative isolates of each B. cereus s. l. species identified was typically found to be reduced in milk, compared with broth. Among the C. perfringens isolates, only two different toxin types were identified, type A and D. Bovine and ovine farms harbored only type A whereas both type A and D were found on two of the three caprine farms.

Conclusions

This study showed that the B. cereus s. l. populations on the sampled farms exhibit a broad diversity in both species and genotypes. The risk of toxin-induced diarrheal illness through consumption of contaminated milk may be limited, in comparison with other food matrices. Type A strains of C. perfringens were the most abundant on dairy farms in Victoria, however type D may be of concern on caprine farms as it can cause enterotoxemia in goats.

Electronic supplementary material

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

Bacillus "next generation" diagnostics: moving from detection toward subtyping and risk-related strain profiling

The highly heterogeneous genus Bacillus comprises the largest species group of endospore forming bacteria. Because of their ubiquitous nature, Bacillus spores can enter food production at several stages resulting in significant economic losses and posing a potential risk to consumers due the capacity of certain Bacillus strains for toxin production. In the past, food microbiological diagnostics was focused on the determination of species using conventional culture-based methods, which are still widely used. However, due to the extreme intra-species diversity found in the genus Bacillus, DNA-based identification and typing methods are gaining increasing importance in routine diagnostics. Several studies showed that certain characteristics are rather strain-dependent than species-specific. Therefore, the challenge for current and future Bacillus diagnostics is not only the efficient and accurate identification on species level but also the development of rapid methods to identify strains with specific characteristics (such as stress resistance or spoilage potential), trace contamination sources, and last but not least discriminate potential hazardous strains from non-toxic strains.

Divergence of the SigB regulon and pathogenesis of the Bacillus cereus sensu lato group

Background

The Bacillus cereus sensu lato group currently includes seven species (B. cereus, B. anthracis, B. mycoides, B. pseudomycoides, B. thuringiensis, B. weihenstephanensis and B. cytotoxicus) that recent phylogenetic and phylogenomic analyses suggest are likely a single species, despite their varied phenotypes. Although horizontal gene transfer and insertion-deletion events are clearly important for promoting divergence among these genomes, recent studies have demonstrated that a major basis for phenotypic diversity in these organisms may be differential regulation of the highly similar gene content shared by these organisms. To explore this hypothesis, we used an in silico approach to evaluate the relationship of pathogenic potential and the divergence of the SigB-dependent general stress response within the B. cereus sensu lato group, since SigB has been demonstrated to support pathogenesis in Bacillus, Listeria and Staphylococcus species.

Results

During the divergence of these organisms from a common “SigB-less” ancestor, the placement of SigB promoters at varied locations in the B. cereus sensu lato genomes predict alternative structures for the SigB regulon in different organisms. Predicted promoter changes suggesting differential transcriptional control of a common gene pool predominate over evidence of indels or horizontal gene transfer for explaining SigB regulon divergence.

Conclusions

Four lineages of the SigB regulon have arisen that encompass different gene contents and suggest different strategies for supporting pathogenesis. This is consistent with the hypothesis that divergence within the B. cereus sensu lato group rests in part on alternative strategies for regulation of a common gene pool.

Molecular characterization of the probiotic strain Bacillus cereus var. toyoi NCIMB 40112 and differentiation from food poisoning strains

Bacillus cereus var. toyoi strain NCIMB 40112 (Toyocerin), a probiotic authorized in the European Union as feed additive for swine, bovines, poultry, and rabbits, was characterized by DNA fingerprinting applying pulsed-field gel electrophoresis and multilocus sequence typing and was compared with reference strains (of clinical and environmental origins). The probiotic strain was clearly characterized by pulsed-field gel electrophoresis using the restriction enzymes Apa I and Sma I resulting in unique DNA patterns. The comparison to the clinical reference strain B. cereus DSM 4312 was done with the same restriction enzymes, and again a clear differentiation of the two strains was possible by the resulting DNA patterns. The use of the restriction enzymes Apa I and Sma I is recommended for further studies. Furthermore, multilocus sequence typing analysis revealed a sequence type (ST 111) that was different from all known STs of B. cereus strains from food poisoning incidents. Thus, a strain characterization and differentiation from food poisoning strains for the probiotic strain was possible.

HyperCAT: an extension of the SuperCAT database for global multi-scheme and multi-datatype phylogenetic analysis of the Bacillus cereus group population

The Bacillus cereus group of bacteria includes species that are of significant medical and economic importance. We previously developed the SuperCAT database, which integrates data from all five multilocus sequence typing (MLST) schemes available to infer the genetic relatedness within this group. Since large numbers of isolates have been typed by other techniques, these should be incorporated in order to provide the most comprehensive and truly global view of the B. cereus group population. The SuperCAT system has been extended into a new database, HyperCAT, with two main additions. First, an extended supertree approach was applied to combine the phylogenetic information available from MLST, amplified fragment length polymorphism and multilocus enzyme electrophoresis. Secondly, a tree-independent clustering algorithm was designed to build superclusters of genetically closely related isolates sharing identical genotyping data. The superclusters were then mapped onto the supertree to generate an integrative genetic and phylogenetic snapshot of the B. cereus group population currently incorporating 2143 isolates. HyperCAT is freely accessible at the University of Oslo's typing website, which has also been upgraded with TNT software, allowing improved and ultra-fast supertree reconstructions. In addition, novel and advanced tools have been included for interactive viewing and navigation of trees, clusters and networks. Database URL: http://mlstoslo.uio.no/

Extended and global phylogenetic view of the Bacillus cereus group population by combination of MLST, AFLP, and MLEE genotyping data

The Bacillus cereus group of bacteria includes species that can cause food-poisoning or spoilage, such as B. cereus, as well as Bacillus anthracis, the cause of anthrax. In the present report we have conducted a multi-datatype analysis using tools from the HyperCAT database (http://mlstoslo.uio.no/) that we recently developed, combining data from multilocus sequence typing (Tourasse et al., 2010), amplified fragment length polymorphism, and multilocus enzyme electrophoresis typing techniques. We provide a comprehensive snapshot of the B. cereus group population, incorporating 2213 isolates including 450 from food and dairy products, in the form of both phylogenetic supertrees and superclusters of genetically closely related isolates. Our main findings include the detection of phylogenetically separated groups of isolates possibly representing novel evolutionary lineages within the B. cereus group, a putative new branch of B. anthracis, as well as new groups of related strains containing both environmental and clinical isolates. In addition, the multi-datatype analysis revealed to a larger extent than previously recognized that food-borne isolates can share identical genotyping profiles with strains from various other origins. Altogether, the global analysis confirms and extends the results underlining the opportunistic nature of B. cereus group organisms, and the fact that isolates responsible for disease outbreaks and contamination of foodstuffs can originate from various genetic backgrounds.

Production, secretion and biological activity of Bacillus cereus enterotoxins

Bacillus cereus behaves as an opportunistic pathogen frequently causing gastrointestinal diseases, and it is increasingly recognized to be responsible for severe local or systemic infections. Pathogenicity of B. cereus mainly relies on the secretion of a wide array of toxins and enzymes and also on the ability to undergo swarming differentiation in response to surface-sensing. In this report, the pathogenicity exerted by B. cereus toxins is described with particular attention to the regulatory mechanisms of production and secretion of HBL, Nhe and CytK enterotoxins.