Intraclade Variability in Toxin Production and Cytotoxicity of Bacillus cereus Group Type Strains and Dairy-Associated Isolates

Assessment of the exposure to cytotoxic Bacillus cereus group genotypes through HTST milk consumption

This study addresses the limited tools available for assessing food safety risks from cytotoxic Bacillus cereus group strains in contaminated food. We quantified the growth, in skim milk broth, of 17 cytotoxic B. cereus strains across 6 phylogenetic groups with various virulence gene profiles. The strains did not grow in HTST milk at 4 or 6°C. At 10°C, 15 strains exhibited growth; at 8°C, one strain grew; and all strains grew at temperatures ≥ 14°C. Using growth data from 16 strains, we developed linear secondary growth models and an exposure assessment model. This model, simulating a 5-stage HTST milk supply chain and up to 35 d of consumer storage with an initial contamination of 100 cfu/mL, estimated that 2.81 ± 0.66% and 4.13 ± 2.53% of milk containers would surpass 105 cfu/mL of B. cereus by d 21 and 35, respectively. A sensitivity analysis identified the initial physiological state of cells (Q0) as the most influential variable affecting predictions for specific isolates. What-if scenarios indicated that increases in mean and variability of consumer storage temperatures significantly affected the predicted B. cereus concentrations in milk. This model serves as an initial tool for risk-based food safety decision making regarding low-level B. cereus contamination.

Molecular characterization and toxigenic profiles of Bacillus cereus isolates from foodstuff and food poisoning outbreaks in Brazil

Bacillus cereus sensu stricto (s.s.) is a well-known foodborne pathogen that produces a range of enterotoxins and is able to cause two different types of foodborne illnesses-the emetic and the diarrheal syndromes. In this study, 54 B. cereus s.s. strains isolated from foodstuff and foods involved in food poisoning outbreaks were characterized according to the presence of toxin-encoding genes, virulence-encoding genes, and panC typing. Most isolates were assigned to panC groups IV (61.1%) and III (25.9%), but members of groups II and V could also be found. Investigation of specific alleles revealed high numbers of isolates carrying toxin and other virulence genes including nheA (100%), nheB (100%), hblA (79.6%), hblC (79.6%), hblD (74.1%), cytK-2 (61.1%), clo (100%), pc-plc (75.9%), sph (68.5%), pi-plc (66.6%), hlyIII (62.9%), and hlyII (24.1%). All isolates were negative for ces and cytK-1. In summary, we detected various enterotoxin and other virulence factor genes associated with diarrheal syndrome in strains analyzed, implicated or not with food poisoning. Furthermore, the most isolates analyzed belong to high-risk phylogenetic groups' panC types III and IV. Our study provides a convenient molecular scheme for characterization of B. cereus s.s. strains responsible for food poisoning outbreaks in order to improve the monitoring and investigation and assess emerging clusters and diversity of strains.

Comparison of the performance of multiple whole-genome sequence-based tools for the identification of Bacillus cereus sensu stricto biovar Thuringiensis

The Bacillus cereus sensu stricto (s.s.) species comprises strains of biovar Thuringiensis (Bt) known for their bioinsecticidal activity, as well as strains with foodborne pathogenic potential. Bt strains are identified (i) based on the production of insecticidal crystal proteins, also known as Bt toxins, or (ii) based on the presence of cry, cyt, and vip genes, which encode Bt toxins. Multiple bioinformatics tools have been developed for the detection of crystal protein-encoding genes based on whole-genome sequencing (WGS) data. However, the performance of these tools is yet to be evaluated using phenotypic data. Thus, the goal of this study was to assess the performance of four bioinformatics tools for the detection of crystal protein-encoding genes. The accuracy of sequence-based identification of Bt was determined in reference to phenotypic microscope-based screening for the production of crystal proteins. A total of 58 diverse B. cereus sensu lato strains isolated from clinical, food, environmental, and commercial biopesticide products underwent WGS. Isolates were examined for crystal protein production using phase contrast microscopy. Crystal protein-encoding genes were detected using BtToxin_Digger, BTyper3, IDOPS (identification of pesticidal sequences), and Cry_processor. Out of 58 isolates, the phenotypic production of crystal proteins was confirmed for 18 isolates. Specificity and sensitivity of Bt identification based on sequences were 0.85 and 0.94 for BtToxin_Digger, 0.97 and 0.89 for BTyper3, 0.95 and 0.94 for IDOPS, and 0.88 and 1.00 for Cry_processor, respectively. Cry_processor predicted crystal protein production with the highest specificity, and BtToxin_Digger and IDOPS predicted crystal protein production with the highest sensitivity. Three out of four tested bioinformatics tools performed well overall, with IDOPS achieving high sensitivity and specificity (>0.90).IMPORTANCEStrains of Bacillus cereus sensu stricto (s.s.) biovar Thuringiensis (Bt) are used as organic biopesticides. Bt is differentiated from the foodborne pathogen Bacillus cereus s.s. by the production of insecticidal crystal proteins. Thus, reliable genomic identification of biovar Thuringiensis is necessary to ensure food safety and facilitate risk assessment. This study assessed the accuracy of whole-genome sequencing (WGS)-based identification of Bt compared to phenotypic microscopy-based screening for crystal protein production. Multiple bioinformatics tools were compared to assess their performance in predicting crystal protein production. Among them, identification of pesticidal sequences performed best overall at WGS-based Bt identification.

Whole Genome Sequencing Reveals Antimicrobial Resistance and Virulence Genes of Both Pathogenic and Non-Pathogenic B. cereus Group Isolates from Foodstuffs in Thailand

Members of the Bacillus cereus group are spore-forming Gram-positive bacilli that are commonly associated with diarrheal or emetic food poisoning. They are widespread in nature and frequently present in both raw and processed food products. Here, we genetically characterized 24 B. cereus group isolates from foodstuffs. Whole-genome sequencing (WGS) revealed that most of the isolates were closely related to B. cereus sensu stricto (12 isolates), followed by B. pacificus (5 isolates), B. paranthracis (5 isolates), B. tropicus (1 isolate), and "B. bingmayongensis" (1 isolate). The most detected virulence genes were BAS_RS06430, followed by bacillibactin biosynthesis genes (dhbA, dhbB, dhbC, dhbE, and dhbF), genes encoding the three-component non-hemolytic enterotoxin (nheA, nheB, and nheC), a gene encoding an iron-regulated leucine-rich surface protein (ilsA), and a gene encoding a metalloprotease (inhA). Various biofilm-associated genes were found, with high prevalences of tasA and sipW genes (matrix protein-encoding genes); purA, purC, and purL genes (eDNA synthesis genes); lytR and ugd genes (matrix polysaccharide synthesis genes); and abrB, codY, nprR, plcR, sinR, and spo0A genes (biofilm transcription regulator genes). Genes related to fosfomycin and beta-lactam resistance were identified in most of the isolates. We therefore demonstrated that WGS analysis represents a useful tool for rapidly identifying and characterizing B. cereus group strains. Determining the genetic epidemiology, the presence of virulence and antimicrobial resistance genes, and the pathogenic potential of each strain is crucial for improving the risk assessment of foodborne B. cereus group strains.

Pathogenicity and Genomic Characterization of a Novel Genospecies, Bacillus shihchuchen, of the Bacillus cereus Group Isolated from Chinese Softshell Turtle (Pelodiscus sinensis)

The Chinese softshell turtle (CST; Pelodiscus sinensis) is a freshwater aquaculture species of substantial economic importance that is commercially farmed across Asia, particularly in Taiwan. Although diseases caused by the Bacillus cereus group (Bcg) pose a major threat to commercial CST farming systems, information regarding its pathogenicity and genome remains limited. Here, we investigated the pathogenicity of Bcg strains isolated in a previous study and performed whole-genome sequencing. Pathogenicity analysis indicated that QF108-045 isolated from CSTs caused the highest mortality rate, and whole-genome sequencing revealed that it was an independent group distinct from other known Bcg genospecies. The average nucleotide identity compared to other known Bcg genospecies was below 95%, suggesting that QF108-045 belongs to a new genospecies, which we named Bacillus shihchuchen. Furthermore, genes annotation revealed the presence of anthrax toxins, such as edema factor and protective antigen, in QF108-045. Therefore, the biovar anthracis was assigned, and the full name of QF108-045 was Bacillus shihchuchen biovar anthracis. In addition to possessing multiple drug-resistant genes, QF108-045 demonstrated resistance to various types of antibiotics, including penicillins (amoxicillin and ampicillin), cephalosporins (ceftifour, cephalexin, and cephazolin), and polypeptides, such as vancomycin.

Modelling growth of Bacillus cereus in paneer by one-step parameter estimation

Bacillus cereus phylogenetic group III and IV strains are commonly associated with food products and cause toxin mediated foodborne diseases. These pathogenic strains have been identified from milk and dairy products, such as reconstituted infant formula and several cheeses. Paneer is a fresh, soft cheese originating from India that is prone to foodborne pathogen contamination, such as by Bacillus cereus. However, there are no reported studies of B. cereus toxin formation in paneer or predictive models quantifying growth of the pathogen in paneer under different environmental conditions. This study assessed enterotoxin-producing potential of B. cereus group III and IV strains, isolated from dairy farm environments, in fresh paneer. Growth of a four-strain cocktail of toxin-producing B. cereus strains was measured in freshly prepared paneer incubated at 5-55 °C and modelled using a one-step parameter estimation combined with bootstrap re-sampling to generate confidence intervals for model parameters. The pathogen grew in paneer between 10 and 50 °C and the developed model fit the observed data well (R² = 0.972, RMSE = 0.321 log₁₀ CFU/g). The cardinal parameters for B. cereus growth in paneer along with the 95% confidence intervals were: μ_opt 0.812 log₁₀ CFU/g/h (0.742, 0.917); T_opt is 44.177 °C (43.16, 45.49); T_min is 4.405 °C (3.973, 4.829); T_max is 50.676 °C (50.367, 51.144). The model developed can be used in food safety management plans and risk assessments to improve safety of paneer while also adding to limited information on B. cereus growth kinetics in dairy products.

Lactobacillus rhamnosus GR-1 attenuates foodborne Bacillus cereus-induced NLRP3 inflammasome activity in bovine mammary epithelial cells by protecting intercellular tight junctions

Background

Bacillus cereus is an important pathogen that causes human food poisoning, specifically diarrhea and vomiting. B. cereus can also induce mastitis in dairy cows and has a strong survival ability in milk, as it cannot be inactivated by high-temperature short-time pasteurization. Therefore, B. cereus can enter the market through pasteurized milk and other dairy products, imposing enormous hidden dangers on food safety and human health.

Results

In this study, B. cereus 2101 (BC) was isolated from milk samples of cows with mastitis. BC grew rapidly with strong hemolysis, making it difficult to prevent mastitis and ensure food security. MAC-T cells were treated with BC and/or Lactobacillus rhamnosus GR-1 (LGR-1). Pretreatment with LGR-1 protected the integrity of tight junctions and the expression of zonula occludens-1 (ZO-1) and occludin destroyed by BC. Furthermore, LGR-1 pretreatment reduced the expression of NOD-like receptor family member pyrin domain-containing protein 3 (NLRP3), caspase recruitment and activation domain (ASC), Caspase-1 p20, gasdermin D (GSDMD) p30, inflammatory factors (interleukin (IL)-1β and IL-18), and cell death induced by BC. Moreover, LGR-1 pretreatment reduced NLRP3 inflammasome activity and increased expressions of ZO-1 and occludin induced by lipopolysaccharides (LPS) + ATP stimulation. MAC-T cells were transfected with NLRP3 siRNA or MCC950 and/or treated with BC and/or LGR-1. NLRP3-siRNA transfection and MCC950 attenuated BC-induced NLRP3 inflammasome activity. Expression of inflammatory cytokines and cell death suggested that the inflammatory pathway might play an important role in the induction of the NLRP3 inflammasome by BC and the protection of LGR-1.

Conclusions

These results suggest that LGR-1 might be a probiotic alternative to antibiotics and could be administered to prevent mastitis in dairy cows, thus ensuring food security.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00752-w.

Prevalence, distribution, enterotoxin profiles, antimicrobial resistance, and genetic diversity of Bacillus cereus group isolates from lettuce farms in Korea

Lettuce wraps are popular in Korean cuisine for their high nutritional value and versatility as healthy additions to multiple dishes. Microbial contamination of lettuce is a major concern, as lettuce is consumed fresh without cooking. Among foodborne pathogens, the spore-forming, facultative anaerobic bacterium, Bacillus cereus is one of the frequently detected pathogen in lettuce in Korea. In this study, we investigated the prevalence and distribution of Bacillus cereus strains in lettuce production farms and further evaluated the enterotoxin gene profiles, antibiotic susceptibility, multidrug resistance pattern, and genetic differences among the B. cereus group isolates. Of the 140 samples isolated from 10 lettuce production farms, 30 samples (21.42%) were positive for B. cereus in which 19 (31.6%) and 10 (23.25%) were from soil and lettuce, respectively. The enterotoxin patterns A (hblCDA, nheABC, entFM, and cytK genes) and B (hblCDA, nheABC, and entFM genes) accounted for 50% and 20% of all the isolates, whereas the emetic gene cesB was not detected in any of the B. cereus group isolates. Antibiotic susceptibility testing of the B. cereus group isolates revealed that all the strains were predominantly resistant to β-lactam antibiotics except imipenem and generally susceptible to most of the non β-lactam antibiotics, including gentamycin, streptomycin, chloramphenicol, and tetracycline. ERIC-PCR and MLST analysis revealed high genetic diversity among the 30 B. cereus group isolates, which belonged to 26 different sequence types (STs) and seven new STs. Moreover, isolates with identical STs exhibited similar patterns of antibiotic resistance and enterotoxin profiles. Results of this study indicate a high prevalence of B. cereus group isolates in lettuce production farms in the Republic of Korea.

Regulation of Enterotoxins Associated with Bacillus cereus Sensu Lato Toxicoinfection

Bacillus cereus sensu lato (s.l.) includes foodborne pathogens, as well as beneficial microorganisms, such as bioinsecticides. Some of the beneficial and commercially used B. cereus s.l. strains have been shown to carry enterotoxin genes, the products of which can cause toxicoinfection in humans. Furthermore, recent epidemiological reports indicated that some bioinsecticidal strains have been linked with foodborne illness outbreaks. This demonstrates the need for improved surveillance of B. cereus s.l., which includes characterization of isolates' virulence capacity. However, the prediction of virulence capacity of B. cereus s.l. strains is challenging. Genetic screening for enterotoxin gene presence has proven to be insufficient for accurate discrimination between virulent and avirulent strains, given that nearly all B. cereus s.l. strains carry at least one enterotoxin gene. Furthermore, complex regulatory networks governing the expression of enterotoxins, and potential synergistic interactions between enterotoxins and other virulence factors make the prediction of toxicoinfection based on isolates' genome sequences challenging. In this review, we summarize and synthesize the current understanding of the regulation of enterotoxins associated with the B. cereus s.l. toxicoinfection and identify gaps in the knowledge that need to be addressed to facilitate identification of genetic markers predictive of cytotoxicity and toxicoinfection.

Bacillus cereus cytotoxin K triggers gasdermin D-dependent pyroptosis

Bacillus cereus is well known as a causative agent of foodborne gastrointestinal diseases and systemic non-gastrointestinal diseases. We have recently identified a pathogenic B. cereus (named H2) from a deep-sea cold-seep. H2 possesses the pyroptosis-inducing capacity and contains a number of enterotoxins including cytotoxin K (CytK). In the present work, we examined the cytotoxicity of the CytK of H2 to human macrophages. CytK bound macrophages by interaction with the plasma membrane and caused cellular structure damage. CytK−cell interaction triggered rapid pyroptosis mediated by caspase 1-activated gasdermin D (GSDMD). CytK-induced pyroptosis required NLRP3 inflammasome activation, K⁺ efflux, and intracellular Ca²⁺ accumulation. CytK exhibited apparent binding to several cytomembrane lipids, in particular phosphatidic acid, which proved to be essential to CytK-elicited cell death. Together, these results add new insights into the cytotoxic mechanism of CytK.

Comparative analysis of Bacillus cereus group isolates' resistance using disk diffusion and broth microdilution and the correlation between antimicrobial resistance phenotypes and genotypes

Bacillus cereus group isolates (n = 85) were screened for phenotypic resistance to 18 antibiotics using broth microdilution and CLSI M45 Bacillus spp. breakpoints. The susceptibility to 9 out of 18 antibiotics was tested also using disk diffusion method and M100 Staphylococcus spp. breakpoints. Overall, a high prevalence of susceptibility to clinically relevant antibiotics was identified using broth microdilution. For most tested antibiotics, a poor correlation was found between zones of inhibition and minimum inhibitory concentrations. Using the broth microdilution results as a reference for comparison, we identified high error rates and low categorical agreement between results produced using disk diffusion and broth microdilution for the seven tested antibiotics with defined breakpoints. This suggests that disk diffusion should be avoided for AST of B. cereus group isolates. Further, we detected antimicrobial resistance genes with ARIBA and ABRIcate to calculate the sensitivity and specificity for predicting phenotypic resistance determined using broth microdilution based on the presence of detected antimicrobial resistance genes (ARGs). ARGs with poor sensitivity and high specificity included rph (rifampicin, 0%, 93%), mph (erythromycin, 0%, 99%), bla1 (penicillin, 29%, 100), and blaZ (penicillin, 56%, 100%). Compared to penicillin, bla1 and blaZ had lower specificity for the prediction of ampicillin resistance. Overall, none of the ARGs had both high sensitivity and specificity, suggesting the need for further study of the mechanisms underlying phenotypic antimicrobial resistance in the B. cereus group. IMPORTANCE Bacillus cereus group includes human pathogens that can cause severe infections requiring antibiotic treatment. Screening of environmental and food isolates for antimicrobial resistance can provide insight into what antibiotics may be more effective therapeutic options based on the lower prevalence of resistance. Currently, the comparison of antimicrobial susceptibility testing results using the disk diffusion method is complicated by the fact that many previous studies have used Staphylococcus spp. breakpoints to interpret their results. In this study, we compared the results of disk diffusion interpreted using the Staphylococcus spp. breakpoints against the results of broth microdilution interpreted using Bacillus spp. breakpoints. We demonstrated that the disk diffusion method does not produce reliable results for B. cereus group isolates and should therefore be avoided. This study also provides new insight into poor associations between the presence of antimicrobial resistance genes and resistance phenotypes for the B. cereus group.

Performance Testing of Bacillus cereus Chromogenic Agar Media for Improved Detection in Milk and Other Food Samples

In this study, the performance of four alternative selective chromogenic B. cereus agar was compared to the reference mannitol-yolk polymyxin (MYP) agar (ISO 7932) using inclusion and exclusion test strains (n = 110) and by analyzing naturally contaminated milk and other food samples (n = 64). Subsequently, the panC group affiliation and toxin gene profile of Bacillus cereus senso lato (s.l.) isolates were determined. Our results corroborate that the overall best performing media CHROMagar™ B. cereus (93.6% inclusivity; 82.7% exclusivity) and BACARA^® (98.2% inclusivity, 62.7% exclusivity) are more sensitive and specific compared to Brilliance™ B. cereus, MYP and ChromoSelect Bacillus Agar. Both media allow unequivocal detection of B. cereus with low risks of misidentification. Media containing ß-D-glucosidase for the detection of presumptive B. cereus may form atypical colony morphologies resulting in a false negative evaluation of the sample. Naturally contaminated samples presented high numbers of background flora, while numbers of presumptive B. cereus were below the detection limit (<10 CFU g⁻¹ or mL⁻¹). Recovery after freezing resulted in the highest detection of B. cereus s.l. on BACARA^® (57.8%), CHROMagar™ B. cereus (56.3%) and MYP agar (54.7%). The panC/toxin profile combination IV/A was the most abundant (33.0%), followed by III/F (21.7%) and VI/C (10.4%). More panC and toxin combinations were present in 15.6% of samples when reanalyzed after freezing. In order to improve detection and confirmation of B. cereus s.l. in food samples, we recommend the parallel use of two complementary selective media followed by molecular characterization (e.g., panC typing combined with toxin gene profiling). When determining psychrotolerant or thermophilic members of the B. cereus group, the selective agar media should additionally be incubated at appropriate temperatures (5 °C, ≥45 °C). If high-risk toxin genes (e.g., ces or cytK-1) are detected, the strain-specific ability to produce toxin should be examined to decisively assess risk.

The Possible Transmission and Potential Enterotoxicity of Bacillus cereus on Lettuce Farms in Five Chinese Provinces

Bacillus cereus is a well-characterized human pathogen that produces toxins associated with diarrheal and emetic foodborne diseases. To investigate the possible transmission of B. cereus on lettuce farms in China and determine its enterotoxicity, (I) a total of 524 samples (lettuce: 332, soil: 69, water: 57, manure: 57, pesticide: 9) were collected from 46 lettuce farms in five Chinese provinces, (II) multilocus sequence typing (MLST) was used to classify B. cereus isolates and for trace analysis, and (III) the presence of toxin genes and enterotoxins (Hbl and Nhe) was detected in 68 strains. The results showed that one hundred and sixty-one lettuce samples (48.5%) tested positive for B. cereus at levels ranging from 10 to 5.3 × 10⁴ CFU/g. Among the environmental sample categories surveyed, the highest positive rate was that of the pesticide samples at 55.6%, followed by soil samples at 52.2% and manure samples at 12.3%. Moreover, one hundred isolates of B. cereus yielded 68 different sequence types (STs) and were classified into five phylogenetic clades. Furthermore, Nhe toxin genes (nheA, nheB, nheC) were broadly distributed and identified in all 68 strains (100%), while Hbl toxin genes (hblA, hblC, hblD) were present in 61 strains (89.7%), entFM was detected in 62 strains (91.2%), and cytK was found in 29 strains (42.6%). All strains were negative for ces. As for the enterotoxin, Nhe was observed in all 68 isolates carrying nheB, while Hbl was present in 76.5% (52/68) of the strains harboring hblC. This study is the first report of possible B. cereus transmission and of its potential enterotoxicity on lettuce farms in China. The results showed that soil and pesticides are the main sources of B. cereus on lettuce farms in China, and the possible transmission routes are as follows: soil-lettuce, manure-lettuce, pesticide-lettuce, manure-soil-lettuce, and water-manure-soil-lettuce. Furthermore, the B. cereus isolates, whether from lettuce or the environment, pose a potential risk to health.

Bacillus cytotoxicus-A potentially virulent food-associated microbe

Bacillus cytotoxicus is a member of the Bacillus cereus group with the ability to grow at high temperatures (up to 52℃) and to synthesize cytotoxin K-1, a diarrhoeagenic cytotoxin, which appears to be unique to this species and more cytotoxic than the cytotoxin K-2 produced by other members of this group. Only a few isolates of this species have been characterized with regard to their cytotoxic effects, and the role of cytotoxin K-1 as a causative agent of food poisoning remains largely unclear. Bacillus cytotoxicus was initially isolated from a food-borne outbreak, which led to three deaths, and the organism has since been linked to other outbreaks all involving plant-based food matrices. Other studies, as well as food-borne incidents reported to the UK Food Standards Agency, detected B. cytotoxicus in insect-related products and in dried food products. With insect-related food becoming increasingly popular, the association with this pathogen is concerning, requiring further investigation and evidence to protect public health. This review summarizes the current knowledge around B. cytotoxicus and highlights gaps in the literature from a food safety perspective.

Bacillus cereus food intoxication and toxicoinfection

Bacillus cereus is one of the leading etiological agents of toxin-induced foodborne diseases. Its omnipresence in different environments, spore formation, and its ability to adapt to varying conditions and produce harmful toxins make this pathogen a health hazard that should not be underestimated. Food poisoning by B. cereus can manifest itself as an emetic or diarrheal syndrome. The former is caused by the release of the potent peptide toxin cereulide, whereas the latter is the result of proteinaceous enterotoxins (e.g., hemolysin BL, nonhemolytic enterotoxin, and cytotoxin K). The final harmful effect is not only toxin and strain dependent, but is also affected by the stress responses, accessory virulence factors, and phenotypic properties under extrinsic, intrinsic, and explicit food conditions and host-related environment. Infamous portrait of B. cereus as a foodborne pathogen, as well as a causative agent of nongastrointestinal infections and even nosocomial complications, has inspired vast volumes of multidisciplinary research in food and clinical domains. As a result, extensive original data became available asking for a new, both broad and deep, multifaceted look into the current state-of-the art regarding the role of B. cereus in food safety. In this review, we first provide an overview of the latest knowledge on B. cereus toxins and accessory virulence factors. Second, we describe the novel taxonomy and some of the most pertinent phenotypic characteristics of B. cereus related to food safety. We link these aspects to toxin production, overall pathogenesis, and interactions with its human host. Then we reflect on the prevalence of different toxinotypes in foods opening the scene for epidemiological aspects of B. cereus foodborne diseases and methods available to prevent food poisoning including overview of the different available methods to detect B. cereus and its toxins.

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.

Characterization of a Deep Sea Bacillus toyonensis Isolate: Genomic and Pathogenic Features

Bacillus toyonensis is a group of Gram-positive bacteria belonging to the Bacillus cereus group and used in some cases as probiotics or biocontrol agents. To our knowledge, B. toyonensis from the deep sea (depth >1,000 m) has not been documented. Here, we report the isolation and characterization of a B. toyonensis strain, P18, from a deep sea hydrothermal field. P18 is aerobic, motile, and able to grow at low temperatures (4°C) and high concentrations of NaCl (8%). P18 possesses a circular chromosome of 5,250,895 bp and a plasmid of 536,892 bp, which encode 5,380 and 523 genes, respectively. Of these genes, 2,229 encode hypothetical proteins that could not be annotated based on the COG database. Comparative genomic analysis showed that P18 is most closely related to the type strain of B. toyonensis, BCT-7112^T. Compared to BCT-7112^T, P18 contains 1,401 unique genes, 441 of which were classified into 20 COG functional categories, and the remaining 960 genes could not be annotated. A total of 319 putative virulence genes were identified in P18, including toxin-related genes, and 24 of these genes are absent in BCT-7112^T. P18 exerted strong cytopathic effects on fish and mammalian cells that led to rapid cell death. When inoculated via injection into fish and mice, P18 rapidly disseminated in host tissues and induced acute infection and mortality. Histopathology revealed varying degrees of tissue lesions in the infected animals. Furthermore, P18 could survive in fish and mouse sera and possessed hemolytic activity. Taken together, these results provide the first evidence that virulent B. toyonensis exists in deep sea environments.

The genomic characterisation and comparison of Bacillus cereus strains isolated from indoor air

Background

Bacillus cereus is ubiquitous in nature, found in environments such as soil, plants, air, and part of the insect and human gut microbiome. The ability to produce endospores and biofilms contribute to their pathogenicity, classified in two types of food poisoning: diarrheal and emetic syndromes. Here we report gap-free, whole-genome sequences of two B. cereus strains isolated from air samples and analyse their emetic and diarrheal potential.

Results

Genome assemblies of the B. cereus strains consist of one chromosome and seven plasmids each. The genome size of strain SGAir0260 is 6.30-Mb with 6590 predicted coding sequences (CDS) and strain SGAir0263 is 6.47-Mb with 6811 predicted CDS. Macrosynteny analysis showed 99% collinearity between the strains isolated from air and 90.2% with the reference genome. Comparative genomics with 57 complete B. cereus genomes suggests these strains from air are closely associated with strains isolated from foodborne illnesses outbreaks. Due to virulence potential of B. cereus and its reported involvement in nosocomial infections, antibiotic resistance analyses were performed and confirmed resistance to ampicillin and fosfomycin, with susceptibility to ciprofloxacin, tetracycline and vancomycin in both strains.

Conclusion

Phylogenetic analysis combined with detection of haemolytic (hblA, hblC, and hblD) and non-haemolytic (nheA, nheB, and nheC) enterotoxin genes in both air-isolated strains point to the diarrheic potential of the air isolates, though not emetic. Characterization of these airborne strains and investigation of their potential disease-causing genes could facilitate identification of environmental sources of contamination leading to foodborne illnesses and nosocomial infections transported by air.

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.

The assessment of leading traits in the taxonomy of the Bacillus cereus group

Bacillus cereus sensu lato strains (B. cereus group) are widely distributed in nature and have received interest for decades due to their importance in insect pest management, food production and their positive and negative repercussions in human health. Consideration of practical uses such as virulence, physiology, morphology, or ill-defined features have been applied to describe and classify species of the group. However, current comparative studies have exposed inconsistencies between evolutionary relatedness and biological significance among genomospecies of the B. cereus group. Here, the combined analyses of core-based phylogeny and all versus all Average Nucleotide Identity values based on 2116 strains were conducted to update the genomospecies circumscriptions within B. cereus group. These analyses suggested the existence of 57 genomospecies, 37 of which are novel, thus indicating that the taxonomic identities of more than 39% of the analyzed strains should be revised or updated. In addition, we found that whole-genome in silico analyses were suitable to differentiate genomospecies such as B. anthracis, B. cereus and B. thuringiensis. The prevalence of toxin and virulence factors coding genes in each of the genomospecies of the B. cereus group was also examined, using phylogeny-aware methods at wide-genome scale. Remarkably, Cry and emetic toxins, commonly assumed to be associated with B. thuringiensis and emetic B. paranthracis, respectively, did not show a positive correlation with those genomospecies. On the other hand, anthrax-like toxin and capsule-biosynthesis coding genes were positively correlated with B. anthracis genomospecies, despite not being present in all strains, and with presumably non-pathogenic genomospecies. Hence, despite these features have been so far considered relevant for industrial or medical classification of related species of the B. cereus group, they were inappropriate for their circumscription. In this study, genomospecies of the group were accurately affiliated and representative strains defined, generating a rational framework that will allow comparative analysis in epidemiological or ecological studies. Based on this classification the role of specific markers such as Type VII secretion system, cytolysin, bacillolysin, and siderophores such as petrobactin were pointed out for further analysis.

Enteropathogenic Potential of Bacillus thuringiensis Isolates from Soil, Animals, Food and Biopesticides

Despite its benefits as biological insecticide, Bacillus thuringiensis bears enterotoxins, which can be responsible for a diarrhoeal type of food poisoning. Thus, all 24 isolates from foodstuffs, animals, soil and commercially used biopesticides tested in this study showed the genetic prerequisites necessary to provoke the disease. Moreover, though highly strain-specific, various isolates were able to germinate and also to actively move, which are further requirements for the onset of the disease. Most importantly, all isolates could grow under simulated intestinal conditions and produce significant amounts of enterotoxins. Cytotoxicity assays classified 14 isolates as highly, eight as medium and only two as low toxic. Additionally, growth inhibition by essential oils (EOs) was investigated as preventive measure against putatively enteropathogenic B. thuringiensis. Cinnamon Chinese cassia showed the highest antimicrobial activity, followed by citral, oregano and winter savory. In all tests, high strain-specific variations appeared and must be taken into account when evaluating the hazardous potential of B. thuringiensis and using EOs as antimicrobials. Altogether, the present study shows a non-negligible pathogenic potential of B. thuringiensis, independently from the origin of isolation. Generally, biopesticide strains were indistinguishable from other isolates. Thus, the use of these pesticides might indeed increase the risk for consumers’ health. Until complete information about the safety of the applied strains and formulations is available, consumers or manufacturers might benefit from the antimicrobial activity of EOs to reduce the level of contamination.

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.

Cereulide Synthetase Acquisition and Loss Events within the Evolutionary History of Group III Bacillus cereus Sensu Lato Facilitate the Transition between Emetic and Diarrheal Foodborne Pathogens

Cereulide-producing members of Bacillus cereus sensu lato group III (also known as emetic B. cereus) possess cereulide synthetase, a plasmid-encoded, nonribosomal peptide synthetase encoded by the ces gene cluster. Despite the documented risks that cereulide-producing strains pose to public health, the level of genomic diversity encompassed by emetic B. cereus has never been evaluated at a whole-genome scale. Here, we employ a phylogenomic approach to characterize group III B. cereus sensu lato genomes which possess ces (ces positive) alongside their closely related, ces-negative counterparts (i) to assess the genomic diversity encompassed by emetic B. cereus and (ii) to identify potential ces loss and/or gain events within the evolutionary history of the high-risk and medically relevant sequence type (ST) 26 lineage often associated with emetic foodborne illness. Using all publicly available ces-positive group III B. cereus sensu lato genomes and the ces-negative genomes interspersed among them (n = 159), we show that emetic B. cereus is not clonal; rather, multiple lineages within group III harbor cereulide-producing strains, all of which share an ancestor incapable of producing cereulide (posterior probability = 0.86 to 0.89). Members of ST 26 share an ancestor that existed circa 1748 (95% highest posterior density [HPD] interval = 1246.89 to 1915.64) and first acquired the ability to produce cereulide before 1876 (95% HPD = 1641.43 to 1946.70). Within ST 26 alone, two subsequent ces gain events were observed, as well as three ces loss events, including among isolates responsible for B. cereus sensu lato toxicoinfection (i.e., "diarrheal" illness).IMPORTANCE B. cereus is responsible for thousands of cases of foodborne disease each year worldwide, causing two distinct forms of illness: (i) intoxication via cereulide (i.e., emetic syndrome) or (ii) toxicoinfection via multiple enterotoxins (i.e., diarrheal syndrome). Here, we show that emetic B. cereus is not a clonal, homogenous unit that resulted from a single cereulide synthetase gain event followed by subsequent proliferation; rather, cereulide synthetase acquisition and loss is a dynamic, ongoing process that occurs across lineages, allowing some group III B. cereus sensu lato populations to oscillate between diarrheal and emetic foodborne pathogens over the course of their evolutionary histories. We also highlight the care that must be taken when selecting a reference genome for whole-genome sequencing-based investigation of emetic B. cereus sensu lato outbreaks, since some reference genome selections can lead to a confounding loss of resolution and potentially hinder epidemiological investigations.

Unique inducible filamentous motility identified in pathogenic Bacillus cereus group species

Active migration across semi-solid surfaces is important for bacterial success by facilitating colonization of unoccupied niches and is often associated with altered virulence and antibiotic resistance profiles. We isolated an atmospheric contaminant, subsequently identified as a new strain of Bacillus mobilis, which showed a unique, robust, rapid, and inducible filamentous surface motility. This flagella-independent migration was characterized by formation of elongated cells at the expanding edge and was induced when cells were inoculated onto lawns of metabolically inactive Campylobacter jejuni cells, autoclaved bacterial biomass, adsorbed milk, and adsorbed blood atop hard agar plates. Phosphatidylcholine (PC), bacterial membrane components, and sterile human fecal extracts were also sufficient to induce filamentous expansion. Screening of eight other Bacillus spp. showed that filamentous motility was conserved amongst B. cereus group species to varying degrees. RNA-Seq of elongated expanding cells collected from adsorbed milk and PC lawns versus control rod-shaped cells revealed dysregulation of genes involved in metabolism and membrane transport, sporulation, quorum sensing, antibiotic synthesis, and virulence (e.g., hblA/B/C/D and plcR). These findings characterize the robustness and ecological significance of filamentous surface motility in B. cereus group species and lay the foundation for understanding the biological role it may play during environment and host colonization.

Potential Enterotoxicity of Phylogenetically Diverse Bacillus cereus Sensu Lato Soil Isolates from Different Geographical Locations

Bacillus cereus sensu lato comprises Gram-positive spore-forming bacteria producing toxins associated with foodborne diseases. Three pore-forming enterotoxins, nonhemolytic enterotoxin (Nhe), hemolysin BL (Hbl), and cytotoxin K (CytK), are considered the primary factors in B. cereus sensu lato diarrhea. The aim of this study was to determine the potential risk of enterotoxicity among soil B. cereus sensu lato isolates representing diverse phylogroups and originated from different geographic locations with various climates (Burkina Faso, Kenya, Argentina, Kazakhstan, and Poland). While nheA- and hblA-positive isolates were present among all B. cereus sensu lato populations and distributed across all phylogenetic groups, cytK-2-positive strains predominated in geographic regions with an arid hot climate (Africa) and clustered together on a phylogenetic tree mainly within mesophilic groups III and IV. The highest in vitro cytotoxicity to Caco-2 and HeLa cells was demonstrated by the strains clustered within phylogroups II and IV. Overall, our results suggest that B. cereus sensu lato pathogenicity is a comprehensive process conditioned by many intracellular factors and diverse environmental conditions.IMPORTANCE This research offers a new route for a wider understanding of the dependency between pathogenicity and phylogeny of a natural bacterial population, specifically within Bacillus cereus sensu lato, that is widely distributed around the world and easily transferred into food products. Our study indicates differences in the phylogenetic and geographical distributions of potential enterotoxigenic B. cereus sensu lato strains. Hence, these bacilli possess a risk for human health, and rapid testing methods for their identification are greatly needed. In particular, the detection of the CytK enterotoxin should be a supporting strategy for the identification of pathogenic B. cereus sensu lato.

Proposal of a Taxonomic Nomenclature for the Bacillus cereus Group Which Reconciles Genomic Definitions of Bacterial Species with Clinical and Industrial Phenotypes

The Bacillus cereus group comprises numerous closely related species, including bioterrorism agent B. anthracis, foodborne pathogen B. cereus, and biopesticide B. thuringiensis Differentiating organisms capable of causing illness or death from those used in industry is essential for risk assessment and outbreak preparedness. However, current species definitions facilitate species-phenotype incongruences, particularly when horizontally acquired genes are responsible for a phenotype. Using all publicly available B. cereus group genomes (n = 2,231), we show that current species definitions lead to overlapping genomospecies clusters, in which 66.2% of genomes belong to multiple genomospecies at a conventional 95 average nucleotide identity (ANI) genomospecies threshold. A genomospecies threshold of ≈92.5 ANI is shown to reflect a natural gap in genome similarity for the B. cereus group, and medoid genomes identified at this threshold are shown to yield resolvable genomospecies clusters with minimal overlap (six of 2,231 genomes assigned to multiple genomospecies; 0.269%). We thus propose a nomenclatural framework for the B. cereus group which accounts for (i) genomospecies using resolvable genomospecies clusters obtained at ≈92.5 ANI, (ii) established lineages of medical importance using a formal collection of subspecies names, and (iii) heterogeneity of clinically and industrially important phenotypes using a formalized and extended collection of biovar terms. We anticipate that the proposed nomenclature will remain interpretable to clinicians, without sacrificing genomic species definitions, which can in turn aid in pathogen surveillance; early detection of emerging, high-risk genotypes; and outbreak preparedness.IMPORTANCE Historical species definitions for many prokaryotes, including pathogens, have relied on phenotypic characteristics that are inconsistent with genome evolution. This scenario forces microbiologists and clinicians to face a tradeoff between taxonomic rigor and clinical interpretability. Using the Bacillus cereus group as a model, a conceptual framework for the taxonomic delineation of prokaryotes which reconciles genomic definitions of species with clinically and industrially relevant phenotypes is presented. The nomenclatural framework outlined here serves as a model for genomics-based bacterial taxonomy that moves beyond arbitrarily set genomospecies thresholds while maintaining congruence with phenotypes and historically important species names.

Enterotoxin Genes, Antibiotic Susceptibility, and Biofilm Formation of Low-Temperature-Tolerant Bacillus cereus Isolated from Green Leaf Lettuce in the Cold Chain

The prevalence and characteristics of low-temperature-tolerant Bacillus cereus (psychrotolerant B. cereus) in green leaf lettuce collected during cold chain were investigated. Among the 101 isolated B. cereus samples, only 18 were capable of growth at 7 °C, and these isolates shared potential health hazard characteristics with mesophilic isolates. Most psychrotolerant B. cereus isolates contained various combinations of nheA, nheB, nheC, hblA, hblA, hblC, hblD, cytK, and entFM. Most isolates of psychrotolerant B. cereus possessed at least two enterotoxin genes and 28% of isolates harbored tested nine enterotoxin genes. Additionally, the psychrotolerant B. cereus isolates showed resistance to tetracycline and rifampin and intermediate levels of resistance to clindamycin. A total of 23% of isolates among psychrotolerant B. cereus displayed a high level of biofilm formation at 7 °C than at 10 °C or 30 °C. The results of this study indicate that cold distribution and storage for green leaf lettuce may fail to maintain food safety due to the presence of enterotoxigenic, antibiotic-resistant, and strong biofilm forming psychrotolerant B. cereus isolates, which therefore poses a potential health risk to the consumer. Our findings provide the first account of the prevalence and characteristics of psychrotolerant B. cereus isolated from green leaf lettuce during cold storage, suggesting a potential hazard of psychrotolerant B. cereus isolates to public health and the food industry.

Food Sensing: Detection of Bacillus cereus Spores in Dairy Products

Milk is a source of essential nutrients for infants and adults, and its production has increased worldwide over the past years. Despite developments in the dairy industry, premature spoilage of milk due to the contamination by Bacillus cereus continues to be a problem and causes considerable economic losses. B. cereus is ubiquitously present in nature and can contaminate milk through a variety of means from the farm to the processing plant, during transport or distribution. There is a need to detect and quantify spores directly in food samples, because B. cereus might be present in food only in the sporulated form. Traditional microbiological detection methods used in dairy industries to detect spores show limits of time (they are time consuming), efficiency and sensitivity. The low level of B. cereus spores in milk implies that highly sensitive detection methods should be applied for dairy products screening for spore contamination. This review describes the advantages and disadvantages of classical microbiological methods used to detect B. cereus spores in milk and milk products, related to novel methods based on molecular biology, biosensors and nanotechnology.

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.

Risk presented to minimally processed chilled foods by psychrotrophic Bacillus cereus

BACKGROUND

Spores of psychrotrophic Bacillus cereus may survive the mild heat treatments given to minimally processed chilled foods. Subsequent germination and cell multiplication during refrigerated storage may lead to bacterial concentrations that are hazardous to health.

SCOPE AND APPROACH

This review is concerned with the characterisation of factors that prevent psychrotrophic B. cereus reaching hazardous concentrations in minimally processed chilled foods and associated foodborne illness. A risk assessment framework is used to quantify the risk associated with B. cereus and minimally processed chilled foods.

KEY FINDINGS AND CONCLUSIONS

Bacillus cereus is responsible for two types of food poisoning, diarrhoeal (an infection) and emetic (an intoxication); however, no reported outbreaks of food poisoning have been associated with B. cereus and correctly stored commercially-produced minimally processed chilled foods. In the UK alone, more than 1010 packs of these foods have been sold in recent years without reported illness, thus the risk presented is very low. Further quantification of the risk is merited, and this requires additional data. The lack of association between diarrhoeal food poisoning and correctly stored commercially-produced minimally processed chilled foods indicates that an infectious dose has not been reached. This may reflect low pathogenicity of psychrotrophic strains. The lack of reported association of psychrotrophic B. cereus with emetic illness and correctly stored commercially-produced minimally processed chilled foods indicates that a toxic dose of the emetic toxin has not been formed. Laboratory studies show that strains form very small quantities of emetic toxin at chilled temperatures.

Toxigenic potential and antimicrobial susceptibility of Bacillus cereus group bacteria isolated from Tunisian foodstuffs

BACKGROUND

Despite the importance of the B. cereus group as major foodborne pathogens that may cause diarrheal and/or emetic syndrome(s), no study in Tunisia has been conducted in order to characterize the pathogenic potential of the B. cereus group. The aim of this study was to assess the sanitary potential risks of 174 B. cereus group strains isolated from different foodstuffs by detecting and profiling virulence genes (hblA, hblB, hblC, hblD, nheA, nheB, nheC, cytK, bceT and ces), testing the isolates cytotoxic activity on Caco-2 cells and antimicrobial susceptibility towards 11 antibiotics.

RESULTS

The entertoxin genes detected among B. cereus isolates were, in decreasing order, nheA (98.9%), nheC (97.7%) and nheB (86.8%) versus hblC (54.6%), hblD (54.6%), hblA (29.9%) and hblB (14.9%), respectively encoding for Non-hemolytic enterotoxin (NHE) and Hemolysin BL (HBL). The isolates are multi-toxigenic, harbouring at least one gene of each NHE and HBL complexes associated or not to bceT, cytK-2 and ces genes. Based on the incidence of virulence genes, the strains were separated into 12 toxigenic groups. Isolates positive for cytK (37,9%) harbored the cytK-2 variant. The detection rates of bceT and ces genes were 50.6 and 4%, respectively. When bacteria were incubated in BHI-YE at 30 °C for 18 h and for 5 d, 70.7 and 35% of the strains were shown to be cytotoxic to Caco-2 cells, respectively. The cytotoxicity of B. cereus strains depended on the food source of isolation. The presence of virulence factors is not always consistent with cytotoxicity. However, different combinations of enterotoxin genetic determinants are significantly associated to the cytotoxic potential of the bacteria. All strains were fully sensitive to rifampicin, chloramphenicol, ciprofloxacin, and gentamycin. The majority of the isolates were susceptible to streptomycin, kanamycin, erythromycin, vancomycin and tetracycline but showed resistance to ampicillin and novobiocin.

CONCLUSION

Our results contribute data that are primary to facilitate risk assessments in order to prevent food poisoning due to B. cereus group.

Comparative Genomic and Phylogenomic Analyses Clarify Relationships Within and Between Bacillus cereus and Bacillus thuringiensis: Proposal for the Recognition of Two Bacillus thuringiensis Genomovars

The present study was designed to clarify the taxonomic status of two species classified as Bacillus cereus sensu lato, namely B. cereus sensu stricto and Bacillus thuringiensis. To this end, nearly 900 whole genome sequences of strains assigned to these taxa were the subject of comparative genomic and phylogenomic analyses. A phylogenomic tree based on core gene sequences showed that the type strains of B. cereus and B. thuringiensis formed a well-supported monophyletic clade that was clearly separated from corresponding clades composed of the remaining validly published species classified as B. cereus sensu lato. However, since average nucleotide identity and digital DNA-DNA hybridization similarities between the two types of Bacillus were slightly higher than the thresholds used to distinguish between closely related species we conclude that B. cereus and B. thuringiensis should continue to be recognized as validly published species. The B. thuringiensis strains were assigned to two genomically distinct groups, we propose that these taxa be recognized as genomovars, that is, as B. thuringiensis gv. thuringiensis and B. thuringiensis gv. cytolyticus. The extensive comparative genomic data clearly show that the distribution of pesticidal genes is irregular as strains identified as B. thuringiensis were assigned to several polyphyletic groups/subclades within the B. cereus-B. thuringiensis clade. Consequently, we recommend that genomic or equivalent molecular systematic features should be used to identify B. thuringiensis strains as the presence of pesticidal genes cannot be used as a diagnostic marker for this species. Comparative taxonomic studies are needed to find phenotypic properties that can be used to distinguish between the B. thuringiensis genomovars and between them and B. cereus.

Precision Food Safety: a Paradigm Shift in Detection and Control of Foodborne Pathogens

The implementation of whole-genome sequencing in food safety has revolutionized foodborne pathogen tracking and outbreak investigations. The vast amounts of genomic data that are being produced through ongoing surveillance efforts continue advancing our understanding of pathogen diversity and genome biology. Produced genomic data are also supporting the use of metagenomics and metatranscriptomics for detection and functional characterization of microbiological hazards in foods and food processing environments. In addition to that, many studies have shown that metabolic and pathogenic potential, antimicrobial resistance, and other phenotypes relevant to food safety can be predicted from whole-genome sequences, omitting the need for multiple laboratory tests. Nevertheless, further work in the area of functional inference is necessary to enable accurate interpretation of functional information inferred from genomic and metagenomic data, as well as real-time detection and tracking of high-risk pathogen subtypes and microbiomes.

Genes Associated With Psychrotolerant Bacillus cereus Group Isolates

The Bacillus cereus group comprises 18 different species, including human pathogens as well as psychrotolerant strains that are an important cause of fluid milk spoilage. To enhance our understanding of the genetic markers associated with psychrotolerance (defined here as > 1 log₁₀ increase in cfu/mL after 21 days incubation at 6°C) among dairy-associated B. cereus group isolates, we used genetic (whole genome sequencing) and phenotypic methods [growth in Skim Milk Broth (SMB) and Brain Heart Infusion (BHI) broth] to characterize 23 genetically-distinct representative isolates from a collection of 503 dairy-associated isolates. Quality threshold clustering identified three categories of psychrotolerance: (i) 14 isolates that were not psychrotolerant in BHI or SMB, (ii) 6 isolates that were psychrotolerant in BHI but not in SMB, and (iii) 2 isolates that were psychrotolerant in BHI and SMB. One isolate, which was psychrotolerant in BHI broth but was just below the cut-off of >1 log₁₀ cfu/mL increase in SMB was not assigned to a cluster. A maximum likelihood phylogeny constructed with core genome single nucleotide polymorphisms classified all psychrotolerant isolates (i.e., psychrotolerant in BHI) into clade VI (representing B. mycoides/weihenstephanensis). Analysis of correlations between gene ortholog presence or absence patterns and psychrotolerance identified 206 orthologous gene clusters that were significantly overrepresented among psychrotolerant strains, including two clusters of cold shock proteins, which were identified in 8/9 and 7/9 psychrotolerant isolates. Gene ontology analyses revealed 36 gene ontology terms that were overrepresented in psychrotolerant isolates, including putrescine catabolic processes and putrescine transmembrane transporter activity. Lastly, Hidden Markov Model searches identified three protein family motifs, including cold shock domain proteins and fatty acid hydroxylases that were significantly associated with psychrotolerance in BHI broth. Analyses of CspA sequences revealed a positive association between psychrotolerant strains and a previously identified “psychrotolerant” CspA sequence. Overall, our data highlight genetic and phenotypic differences in psychrotolerance among B. cereus group dairy-associated isolates and show that psychrotolerance is dependent on the growth medium. We also identified a number of gene targets that could be used for specific detection or control of psychrotolerant B. cereus group isolates.

Characterization of Emetic and Diarrheal Bacillus cereus Strains From a 2016 Foodborne Outbreak Using Whole-Genome Sequencing: Addressing the Microbiological, Epidemiological, and Bioinformatic Challenges

The Bacillus cereus group comprises multiple species capable of causing emetic or diarrheal foodborne illness. Despite being responsible for tens of thousands of illnesses each year in the U.S. alone, whole-genome sequencing (WGS) is not yet routinely employed to characterize B. cereus group isolates from foodborne outbreaks. Here, we describe the first WGS-based characterization of isolates linked to an outbreak caused by members of the B. cereus group. In conjunction with a 2016 outbreak traced to a supplier of refried beans served by a fast food restaurant chain in upstate New York, a total of 33 B. cereus group isolates were obtained from human cases (n = 7) and food samples (n = 26). Emetic (n = 30) and diarrheal (n = 3) isolates were most closely related to B. paranthracis (group III) and B. cereus sensu stricto (group IV), respectively. WGS indicated that the 30 emetic isolates (24 and 6 from food and humans, respectively) were closely related and formed a well-supported clade distinct from publicly available emetic group III genomes with an identical sequence type (ST 26). The 30 emetic group III isolates from this outbreak differed from each other by a mean of 8.3 to 11.9 core single nucleotide polymorphisms (SNPs), while differing from publicly available emetic group III ST 26 B. cereus group genomes by a mean of 301.7-528.0 core SNPs, depending on the SNP calling methodology used. Using a WST-1 cell proliferation assay, the strains isolated from this outbreak had only mild detrimental effects on HeLa cell metabolic activity compared to reference diarrheal strain B. cereus ATCC 14579. We hypothesize that the outbreak was a single source outbreak caused by emetic group III B. cereus belonging to the B. paranthracis species, although food samples were not tested for presence of the emetic toxin cereulide. In addition to showcasing how WGS can be used to characterize B. cereus group strains linked to a foodborne outbreak, we also discuss potential microbiological and epidemiological challenges presented by B. cereus group outbreaks, and we offer recommendations for analyzing WGS data from the isolates associated with them.

Identification of cereulide producing Bacillus cereus by MALDI-TOF MS

The Bacillus (B.) cereus group consists of nine recognized species which are present worldwide. B. cereus play an important role in food-borne diseases by producing different toxins. Yet, only a small percentage of B. cereus strains are able to produce the heat stable cereulide, the causative agent of emetic food poisoning. To minimize the entry of emetic B. cereus into the food chain, food business operators are dependent on efficient and reliable methods enabling the differentiation between emetic and non-emetic strains. Currently, only time-consuming cell bioassays, molecular methods and tandem mass spectrometry are available for this purpose. Thus, the aim of the present study was to establish a fast and reliable method for the differentiation between emetic/non-emetic strains by MALDI-TOF MS. Selected strains/isolates of the B. cereus group as well as other Bacillus spp. (total n = 121) were cultured on sheep blood agar for 48 h before analysis. Subsequently, the cultures were directly analyzed by MALDI-TOF MS without prior extraction steps. The samples were measured in the mass range of m/z 800-1800 Da. Using ClinProTools 3.0 statistical software and Flex analysis software (Bruker Daltonics GmbH, Bremen, Germany), a differentiation between emetic/non-emetic isolates was possible with a rate of correct identification of 99.1% by means of the evaluation of two specific biomarkers (m/z 1171 and 1187 Da).

Enterotoxin Production of Bacillus thuringiensis Isolates From Biopesticides, Foods, and Outbreaks

While the relevance of Bacillus (B.) cereus as a major cause of gastroenteritis is undisputed, the role of the closely related B. thuringiensis in foodborne disease is unclear. B. thuringiensis strains frequently harbor enterotoxin genes. However, the organism has only very rarely been associated with foodborne outbreaks, possibly due to the fact that during outbreak investigations, B. cereus is routinely not differentiated from B. thuringiensis. A recent EFSA scientific opinion stresses the urgent need for further data allowing for improved risk assessment, in particular as B. thuringiensis is a commonly used biopesticide. Therefore, the aim of this study was to gain further insights into the hazardous potential of B. thuringiensis. To this end, 39 B. thuringiensis isolates obtained from commercially used biopesticides, various food sources, as well as from foodborne outbreaks were characterized by panC typing, panC-based SplitsTree analysis, toxin gene profiling, FTIR spectroscopic analysis, a cytotoxicity assay screening for enterotoxic activity, and a sphingomyelinase assay. The majority of the tested B. thuringiensis isolates exhibited low (23%, n = 9) or mid level enterotoxicity (74%, n = 29), and produced either no (59%, n = 23) or low levels (33%, n = 13) of sphingomyelinase, which is reported to act synergistically with enterotoxins Nhe and Hbl. One strain isolated from rosemary was however classified as highly enterotoxic surpassing the cytotoxic activity of the high-level reference strain by a factor of 1.5. This strain also produced vast amounts of sphingomyelinase. Combining all results obtained in this study into a fingerprint pattern, several enterotoxic biopesticide strains were indistinguishable from those of isolates from foods or collected in association with outbreaks. Our study shows that many B. thuringiensis biopesticide strains exhibit mid-level cytotoxicity in a Vero cell assay and that some of these strains cannot be differentiated from isolates collected from foods or in association with outbreaks. Thus, we demonstrate that the use of B. thuringiensis strains as biopesticides can represent a food safety risk, underpinning the importance of assessing the hazardous potential of each strain and formulation used.

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.