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

Cereulide production capacities and genetic properties of 31 emetic Bacillus cereus group strains

The highly potent toxin cereulide is a frequent cause of foodborne intoxications. This extremely resistant toxin is produced by Bacillus cereus group strains carrying the plasmid encoded cesHPTABCD gene cluster. It is known that the capacities to produce cereulide vary greatly between different strains but the genetic background of these variations is not clear. In this study, cereulide production capacities were associated with genetic characteristics. For this, cereulide levels in cultures of 31 strains were determined after incubation in tryptic soy broth for 24 h at 24 °C, 30 °C and 37 °C. Whole genome sequencing based data were used for an in-depth characterization of gene sequences related to cereulide production. The taxonomy, population structure and phylogenetic relationships of the strains were evaluated based on average nucleotide identity, multi-locus sequence typing (MLST), core genome MLST and single nucleotide polymorphism analyses. Despite a limited strain number, the approach of a genome wide association study (GWAS) was tested to link genetic variation with cereulide quantities. Our study confirms strain-dependent differences in cereulide production. For most strains, these differences were not explainable by sequence variations in the cesHPTABCD gene cluster or the regulatory genes abrB, spo0A, codY and pagRBc. Likewise, the population structure and phylogeny of the tested strains did not comprehensively reflect the cereulide production capacities. GWAS yielded first hints for associated proteins, while their possible effect on cereulide synthesis remains to be further investigated.

Analysis of Microbial Composition of Edible Insect Products Available for Human Consumption within the United States Using Traditional Microbiological Methods and Whole Genome Sequencing

Edible insects offer a promising protein source for humans, but their food safety risks have not been previously investigated within the United States. Therefore, the aim of this study was to investigate the microbial content of processed edible insect products. A total of eight different types of edible insect products, including diving beetles, silkworms, grasshoppers, Jamaican crickets, mealworms, mole crickets, whole roasted crickets, and 100% pure cricket powder, were purchased from a large online retailer for the analysis. All the products were purchased in August 2022 and examined between August 2022 and November 2022. Traditional microbiological methods were employed to determine microbial counts for each product type using three replicates (total number of samples = 24). This included assessing aerobic bacterial spore, lactic acid bacteria, Enterobacteriaceae, total viable counts, and the presence of Salmonella. Additionally, whole genome sequencing was employed to further characterize selected colonies (n = 96). Microbial counts data were statistically analyzed using one-way ANOVA, while sequence data were taxonomically classified using Sepia.Bacilluscereusgroup isolates underwent additional characterization with Btyper3. Product type significantly influenced total viable counts, bacterial spore counts, and lactic acid bacteria counts (P = 0.00391, P = 0.0065, and P < 0.001, respectively), with counts ranging from < 1.70 to 6.01 Log10 CFU/g, <1.70 to 5.25 Log10 CFU/g, and < 1.70 to 4.86 Log10 CFU/g, respectively. Enterobacteriaceae were only detected in mole crickets (<2.30 Log10 CFU/g) and house cricket powder (<2.15 Log10 CFU/g). All samples were negative for Salmonella. Whole genome sequencing revealed the presence of 12 different bacterial genera among the analyzed isolates, with a majority belonging to the Bacillus genus. Some of the isolates of Bacillus cereus group were identified as biovar Emeticus. Overall, although edible insects offer a promising food alternative, the presence of Bacillus cereus group in some products could raise concerns regarding food safety.

Microbial pesticides - challenges and future perspectives for testing and safety assessment with respect to human health

Plant protection measures are necessary to prevent pests and diseases from attacking and destroying crop plants and to meet consumer demands for agricultural produce. In the last decades the use of chemical pesticides has largely increased. Farmers are looking for alternatives. Biopesticides should be considered a sustainable solution. They may be less toxic than chemical pesticides, be very specific to the target pest, decompose quickly, and be less likely to cause resistance. On the other hand, lower efficacy and higher costs are two disadvantages of many biopesticides. Biopesticides include macroorganisms, natural compounds and microorganisms. Microbial pesticides are the most widely used and studied class of biopesticides. The greatest difference between microbial and chemical pesticides is the ability of the former to potentially multiply in the environment and on the crop plant after application. The data requirements for the European Union and the United States Environmental Protection Agency are highlighted, as these regulatory processes are the most followed in regions where local regulations for biopesticide products are not available or vague. New Approach Methods already proposed or harmonized for chemical pesticides are presented and discussed with respect to their use in evaluating microbial pesticide formulations. Evaluating the microbials themselves is not as simple as using the same validated New Approach Methods as for synthetic pesticides. Therefore, the authors suggest considering New Approach Method strategies specifically for microbials and global harmonization with acceptability with the advancements of such approaches. Further discussion is needed and greatly appreciated by the experts.

A Burkholderia cenocepacia-like environmental isolate strongly inhibits the plant fungal pathogen Zymoseptoria tritici

Fungal phytopathogens cause significant reductions in agricultural yields annually, and overusing chemical fungicides for their control leads to environmental pollution and the emergence of resistant pathogens. Exploring natural isolates with strong antagonistic effects against pathogens can improve our understanding of their ecology and develop new treatments for the future. We isolated and characterized a novel bacterial strain associated with the species Burkholderia cenocepacia, termed APO9, which strongly inhibits Zymoseptoria tritici, a commercially important pathogenic fungus causing Septoria tritici blotch in wheat. Additionally, this strain exhibits inhibitory activity against four other phytopathogens. We found that physical contact plays a crucial role for APO9's antagonistic capacity. Genome sequencing of APO9 and biosynthetic gene cluster (BGC) analysis identified nine classes of BGCs and three types of secretion systems (types II, III, and IV), which may be involved in the inhibition of Z. tritici and other pathogens. To identify genes driving APO9's inhibitory activity, we screened a library containing 1,602 transposon mutants and identified five genes whose inactivation reduced inhibition efficiency. One such gene encodes for a diaminopimelate decarboxylase located in a terpenoid biosynthesis gene cluster. Phylogenetic analysis revealed that while some of these genes are also found across the Burkholderia genus, as well as in other Betaproteobacteria, the combination of these genes is unique to the Burkholderia cepacia complex. These findings suggest that the inhibitory capacity of APO9 is complex and not limited to a single mechanism, and may play a role in the interaction between various Burkholderia species and various phytopathogens within diverse plant ecosystems.

IMPORTANCE

The detrimental effects of fungal pathogens on crop yields are substantial. The overuse of chemical fungicides contributes not only to environmental pollution but also to the emergence of resistant pathogens. Investigating natural isolates with strong antagonistic effects against pathogens can improve our understanding of their ecology and develop new treatments for the future. We discovered and examined a unique bacterial strain that demonstrates significant inhibitory activity against several phytopathogens. Our research demonstrates that this strain has a wide spectrum of inhibitory actions against plant pathogens, functioning through a complex mechanism. This plays a vital role in the interactions between plant microbiota and phytopathogens.

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.

Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

New Approach Methods to Assess the Enteropathogenic Potential of Strains of the Bacillus cereus Group, including Bacillus thuringiensis

Bacillus cereus (Bc) is a wide group of Gram-positive and spore-forming bacteria, known to be the etiological agents of various human infections, primarily food poisoning. The Bc group includes enteropathogenic strains able to germinate in the digestive tract and to produce enterotoxins such as Nhe, Hbl, and CytK. One species of the group, Bacillus thuringiensis (Bt), has the unique feature of producing insecticidal crystals during sporulation, making it an important alternative to chemical pesticides to protect crops from insect pest larvae. Nevertheless, several studies have suggested a link between the ingestion of pesticide strains and human cases of food poisoning, calling their safety into question. Consequently, reliable tools for virulence assessment are worth developing to aid decision making in pesticide regulation. Here, we propose complementary approaches based on two biological models, the human intestinal Caco-2 cell line and the insect Drosophila melanogaster, to assess and rank the enteric virulence potency of Bt strains in comparison with other Bc group members. Using a dataset of 48 Bacillus spp. strains, we showed that some Bc group strains, including Bt, were able to induce cytotoxicity in Caco-2 cells with concomitant release of IL-8 cytokine, a landmark of pro-inflammatory response. In the D. melanogaster model, we were able to sort a panel of 39 strains into four different classes of virulence, ranging from no virulence to strong virulence. Importantly, for the most virulent strains, mortality was associated with a loss of intestinal barrier integrity. Interestingly, although strains can share a common toxinotype, they display different degrees of virulence, suggesting the existence of specific mechanisms of virulence expression in vivo in the intestine.

(Pan)genomic analysis of two Rhodococcus isolates and their role in phenolic compound degradation

The genus Rhodococcus is recognized for its potential to degrade a large range of aromatic substances, including plant-derived phenolic compounds. We used comparative genomics in the context of the broader Rhodococcus pan-genome to study genomic traits of two newly described Rhodococcus strains (type-strain Rhodococcus pseudokoreensis R79T and Rhodococcus koreensis R85) isolated from apple rhizosphere. Of particular interest was their ability to degrade phenolic compounds as part of an integrated approach to treat apple replant disease (ARD) syndrome. The pan-genome of the genus Rhodococcus based on 109 high-quality genomes was open with a small core (1.3%) consisting of genes assigned to basic cell functioning. The range of genome sizes in Rhodococcus was high, from 3.7 to 10.9 Mbp. Genomes from host-associated strains were generally smaller compared to environmental isolates which were characterized by exceptionally large genome sizes. Due to large genomic differences, we propose the reclassification of distinct groups of rhodococci like the Rhodococcus equi cluster to new genera. Taxonomic species affiliation was the most important factor in predicting genetic content and clustering of the genomes. Additionally, we found genes that discriminated between the strains based on habitat. All members of the genus Rhodococcus had at least one gene involved in the pathway for the degradation of benzoate, while biphenyl degradation was mainly restricted to strains in close phylogenetic relationships with our isolates. The ~40% of genes still unclassified in larger Rhodococcus genomes, particularly those of environmental isolates, need more research to explore the metabolic potential of this genus.IMPORTANCERhodococcus is a diverse, metabolically powerful genus, with high potential to adapt to different habitats due to the linear plasmids and large genome sizes. The analysis of its pan-genome allowed us to separate host-associated from environmental strains, supporting taxonomic reclassification. It was shown which genes contribute to the differentiation of the genomes based on habitat, which can possibly be used for targeted isolation and screening for desired traits. With respect to apple replant disease (ARD), our isolates showed genome traits that suggest potential for application in reducing plant-derived phenolic substances in soil, which makes them good candidates for further testing against ARD.

The evolutionary genomics of adaptation to stress in wild rhizobium bacteria

Microbiota comprise the bulk of life's diversity, yet we know little about how populations of microbes accumulate adaptive diversity across natural landscapes. Adaptation to stressful soil conditions in plants provides seminal examples of adaptation in response to natural selection via allelic substitution. For microbes symbiotic with plants however, horizontal gene transfer allows for adaptation via gene gain and loss, which could generate fundamentally different evolutionary dynamics. We use comparative genomics and genetics to elucidate the evolutionary mechanisms of adaptation to physiologically stressful serpentine soils in rhizobial bacteria in western North American grasslands. In vitro experiments demonstrate that the presence of a locus of major effect, the nre operon, is necessary and sufficient to confer adaptation to nickel, a heavy metal enriched to toxic levels in serpentine soil, and a major axis of environmental soil chemistry variation. We find discordance between inferred evolutionary histories of the core genome and nreAXY genes, which often reside in putative genomic islands. This suggests that the evolutionary history of this adaptive variant is marked by frequent losses, and/or gains via horizontal acquisition across divergent rhizobium clades. However, different nre alleles confer distinct levels of nickel resistance, suggesting allelic substitution could also play a role in rhizobium adaptation to serpentine soil. These results illustrate that the interplay between evolution via gene gain and loss and evolution via allelic substitution may underlie adaptation in wild soil microbiota. Both processes are important to consider for understanding adaptive diversity in microbes and improving stress-adapted microbial inocula for human use.

Exploring the Genomic Landscape of Bacillus paranthracis PUMB_17 as a Proficient Phosphatidylcholine-Specific Phospholipase C Producer

Phospholipases find versatile applications across industries, including detergent production, food modification, pharmaceuticals (especially in drug delivery systems), and cell signaling research. In this study, we present a strain of Bacillus paranthracis for the first time, demonstrating significant potential in the production of phosphatidylcholine-specific phospholipase C (PC-PLC). The investigation thoroughly examines the B. paranthracis PUMB_17 strain, focusing on the activity of PC-PLC and its purification process. Notably, the PUMB_17 strain displays extracellular PC-PLC production with high specific activity during the late exponential growth phase. To unravel the genetic makeup of PUMB_17, we employed nanopore-based whole-genome sequencing and subsequently conducted a detailed genome annotation. The genome comprises a solitary circular chromosome spanning 5,250,970 bp, featuring a guanine-cytosine ratio of 35.49. Additionally, two plasmids of sizes 64,250 bp and 5845 bp were identified. The annotation analysis reveals the presence of 5328 genes, encompassing 5186 protein-coding sequences, and 142 RNA genes, including 39 rRNAs, 103 tRNAs, and 5 ncRNAs. The aim of this study was to make a comprehensive genomic exploration that promises to enhance our understanding of the previously understudied and recently documented capabilities of Bacillus paranthracis and to shed light on a potential use of the strain in the industrial production of PC-PLC.

Investigating environmental transmission to resolve a Bacillus cereus group outbreak in a neonatal intensive care unit using core genome multilocus sequence typing

BACKGROUND

We analyzed an outbreak of Bacillus cereus group (Bcg) at a single-center neonatal intensive care unit level IV by conducting comprehensive sampling of both patients and the environment.

METHODS

Between 06/2020 and 10/2021, all Bcg isolates identified by both regular colonization screening and additional sampling of the environment were subjected to whole-genome sequencing, followed by in vitro extraction of MLST ST, resistance genes and virulence factors. Using publicly available genome sequences, we defined an ad hoc core genome multilocus sequence typing (cgMLST) scheme comprising 2759 target genes for Bcg typing, which we applied to the detected isolates. We have compared the results with a stable cgMLST that was published in the meantime and completed the investigation with a SNP analysis.

RESULTS

We analyzed 28 Bcg isolates from patient and environmental samples using MLST and cgMLST. This revealed multiple sequence types, with ST127 being the most common (n = 13). Both cgMLST schemes grouped ten of the 13 ST127 isolates into a cluster, including two invasive isolates from two different patients and several environmental samples. SNP analysis postulated a screen from a ventilation machine as a possible reservoir.

CONCLUSION

In sensitive settings such as neonatal intensive care units, considering the environment in outbreak analyses is crucial, especially when investigating potential transmission routes through shared devices. When dealing with widespread bacteria such as Bcg, high-resolution typing techniques are necessary. In this study, we successfully resolved an outbreak of Bcg infections using a custom cgMLST scheme combined with a SNP analysis.

Bacillus thuringiensis Bt_UNVM-84, a Novel Strain Showing Insecticidal Activity against Anthonomus grandis Boheman (Coleoptera: Curculionidae)

Bacillus thuringiensis is a Gram-positive bacterium known for its insecticidal proteins effective against various insect pests. However, limited strains and proteins target coleopteran pests like Anthonomous grandis Boheman, causing substantial economic losses in the cotton industry. This study focuses on characterizing a Bacillus sp. strain, isolated from Oncativo (Argentina), which exhibits ovoid to amorphous parasporal crystals and was designated Bt_UNVM-84. Its genome encodes genes for the production of two pairs of binary Vpb1/Vpa2 proteins and three Cry-like proteins showing similarity with different Cry8 proteins. Interestingly, this gene content was found to be conserved in a previously characterized Argentine isolate of B. thuringiensis designated INTA Fr7-4. SDS-PAGE analysis revealed a major band of 130 kDa that is proteolytically processed to an approximately 66-kDa protein fragment by trypsin. Bioassays performed with spore-crystal mixtures demonstrated an interesting insecticidal activity against the cotton boll weevil A. grandis neonate larvae, resulting in 91% mortality. Strain Bt_UNVM-84 is, therefore, an interesting candidate for the efficient biological control of this species, causing significant economic losses in the cotton industry in the Americas.

Analysis of Sporulation in Bacillus cereus Biovar anthracis Which Contains an Insertion in the Gene for the Sporulation Factor σK

Bacillus cereus biovar anthracis (Bcbva) is an untypical pathogen causing a fatal anthrax-like disease in a variety of wildlife species in African rainforest areas. In contrast to Bacillus anthracis and most species of the B. cereus group, all strains of the Bcbva cluster contain a 22 kb insertion in the sigK gene which encodes the essential late sporulation sigma factor σK. This insertion is excised during sporulation in a site-specific recombination process resulting in an intact sigK gene and a circular molecule. The sporulation kinetics of two strains each of Bcbva and B. anthracis were compared by the expression analysis of eight sporulation-associated genes, including sigK, using reverse transcriptase quantitative real-time PCR. In addition, morphological sporulation stages were analyzed and quantified by electron microscopy. Our results indicated that the necessary excision of the insertion in Bcbva neither delayed nor inhibited its sporulation. In two spontaneous mutants of Bcbva, the excision of the sigK insertion and sporulation were impeded due to mutations in the spo0A and spoVG regulator genes, respectively. The spo0A frameshift mutation was overcome by intragenic suppression in a revertant which was able to sporulate normally, despite an M171S amino acid exchange in the global regulator Spo0A. A screening of the NCBI database identified further strains of the B. cereus group which possess unrelated insertions in the sigK gene, and two strains containing almost identical insertions at the same gene position. Some of the sigK insertions encode putative prophages, whereas the Bcbva insertion encoded a type I restriction-modification system. The function of these insertions and if they are possibly essential for sporulation remains to be assessed.

Get to Know Your Neighbors: Characterization of Close Bacillus anthracis Isolates and Toxin Profile Diversity in the Bacillus cereus Group

Unexpected atypical isolates of Bacillus cereus s.l. occasionally challenge conventional microbiology and even the most advanced techniques for anthrax detection. For anticipating and gaining trust, 65 isolates of Bacillus cereus s.l. of diverse origin were sequenced and characterized. The BTyper3 tool was used for assignation to genomospecies B. mosaicus (34), B. cereus s.s (29) and B. toyonensis (2), as well as virulence factors and toxin profiling. None of them carried any capsule or anthrax-toxin genes. All harbored the non-hemolytic toxin nheABC and sphygomyelinase spH genes, whereas 41 (63%), 30 (46%), 11 (17%) and 6 (9%) isolates harbored cytK-2, hblABCD, cesABCD and at least one insecticidal toxin gene, respectively. Matrix-assisted laser desorption ionization-time of flight mass spectrometry confirmed the production of cereulide (ces genes). Phylogeny inferred from single-nucleotide polymorphisms positioned isolates relative to the B. anthracis lineage. One isolate (BC38B) was of particular interest as it appeared to be the closest B. anthracis neighbor described so far. It harbored a large plasmid similar to other previously described B. cereus s.l. megaplasmids and at a lower extent to pXO1. Whereas bacterial collection is enriched, these high-quality public genetic data offer additional knowledge for better risk assessment using future NGS-based technologies of detection.

Characterization of the Toxigenic Potential of Bacillus cereus sensu lato Isolated from Raw Berries and Their Products

Bacillus cereus is estimated to be responsible for 1.4-12% of all food poisoning outbreaks worldwide. The objective of this study was to investigate the toxigenic potential of 181 isolates of B. cereus previously recovered from different types of berries and berry products (strawberries, raspberries, blackberries, and blueberries) by assessing the presence of enterotoxin genes (hblA, hblC, hblD, nheA, nheB, nheC, and cytK) and an emetic toxin cereulide synthetase gene (ces). The cytotoxic activity on Caco-2 cells was also evaluated for the two isolates containing the gene cytK. Twenty-three toxigenic profiles were found. The nheABC (91.7%) and hblACD (89.0%) complexes were the most prevalent among the isolates, while the cytK and ces genes were detected in low percentages, 1.1% and 3.3%, respectively. In addition, the nheABC/hblACD complex and ces genes were detected in isolates recovered throughout the production process of blackberries and strawberries. The cytotoxic activity on Caco-2 cells was also observed to be greater than 60% for isolates containing the cytK gene.

Identification, phylogenetic analysis, and genome mining of the tetracycline-resistant Bacillus thuringiensis strain m401 reveal its potential for biotechnological and biocontrol applications

Bacillus thuringiensis is an entomopathogen belonging to the Bacillus cereus clade. We isolated a tetracycline-resistant strain called m401, recovered it from honey, and identified it as Bacillus thuringiensis sv. kumamotoensis based on the average nucleotide identity calculations (ANIb) comparison and the analysis of the gyrB gene sequences of different B. thuringiensis serovars. Sequences with homology to virulence factors [cytK, nheA, nheB, nheC, hblA, hblB, hblC, hblD, entFM, and inhA] and tetracycline resistance genes [tet(45), tet(V), and tet(M)/tet(W)/tet(O)/tet(S) family] were identified in the bacterial chromosome. The prediction of plasmid-coding regions revealed homolog sequences to the MarR and TetR/AcrR family of transcriptional regulators, toxins, and lantipeptides. The genome mining analysis revealed 12 regions of biosynthetic gene clusters responsible for synthesizing secondary metabolites. We identified biosynthetic gene clusters coding for bacteriocins, siderophores, ribosomally synthesized post-translationally modified peptide products, and non-ribosomal peptide synthetase clusters that provide evidence for the possible use of Bt m401 as a biocontrol agent. Furthermore, Bt m401 showed high inhibition against all Paenibacillus larvae genotypes tested in vitro. In conclusion, Bt m401 owns various genes involved in different biological processes, such as transductional regulators associated with antibiotic resistance, toxins, and antimicrobial peptides with potential biotechnological and biocontrol applications.

Genus-wide genomic characterization of Macrococcus: insights into evolution, population structure, and functional potential

Introduction

Macrococcus species have been isolated from a range of mammals and mammal-derived food products. While they are largely considered to be animal commensals, Macrococcus spp. can be opportunistic pathogens in both veterinary and human clinical settings. This study aimed to provide insight into the evolution, population structure, and functional potential of the Macrococcus genus, with an emphasis on antimicrobial resistance (AMR) and virulence potential.

Methods

All high-quality, publicly available Macrococcus genomes (n = 104, accessed 27 August 2022), plus six South African genomes sequenced here (two strains from bovine clinical mastitis cases and four strains from beef products), underwent taxonomic assignment (using four different approaches), AMR determinant detection (via AMRFinderPlus), and virulence factor detection (using DIAMOND and the core Virulence Factor Database).

Results

Overall, the 110 Macrococcus genomes were of animal commensal, veterinary clinical, food-associated (including food spoilage), and environmental origins; five genomes (4.5%) originated from human clinical cases. Notably, none of the taxonomic assignment methods produced identical results, highlighting the potential for Macrococcus species misidentifications. The most common predicted antimicrobial classes associated with AMR determinants identified across Macrococcus included macrolides, beta-lactams, and aminoglycosides (n = 81, 61, and 44 of 110 genomes; 73.6, 55.5, and 40.0%, respectively). Genes showing homology to Staphylococcus aureus exoenzyme aureolysin were detected across multiple species (using 90% coverage, n = 40 and 77 genomes harboring aureolysin-like genes at 60 and 40% amino acid [AA] identity, respectively). S. aureus Panton-Valentine leucocidin toxin-associated lukF-PV and lukS-PV homologs were identified in eight M. canis genomes (≥40% AA identity, >85% coverage). Using a method that delineates populations using recent gene flow (PopCOGenT), two species (M. caseolyticus and M. armenti) were composed of multiple within-species populations. Notably, M. armenti was partitioned into two populations, which differed in functional potential (e.g., one harbored beta-lactamase family, type II toxin-antitoxin system, and stress response proteins, while the other possessed a Type VII secretion system; PopCOGenT p < 0.05).

Discussion

Overall, this study leverages all publicly available Macrococcus genomes in addition to newly sequenced genomes from South Africa to identify genomic elements associated with AMR or virulence potential, which can be queried in future experiments.

Development of an Online Genome Sequence Comparison Resource for Bacillus cereus sensu lato Strains Using the Efficient Composition Vector Method

An automated method was developed for differentiating closely related B. cereus sensu lato (s.l.) species, especially biopesticide Bacillus thuringiensis, from other human pathogens, B. anthracis and B. cereus sensu stricto (s.s.). In the current research, four typing methods were initially compared, including multi-locus sequence typing (MLST), single-copy core genes phylogenetic analysis (SCCGPA), dispensable genes content pattern analysis (DGCPA) and composition vector tree (CVTree), to analyze the genomic variability of 23 B. thuringiensis strains from aizawai, kurstaki, israelensis, thuringiensis and morrisoni serovars. The CVTree method was the best option to be used for typing B. thuringiensis strains since it proved to be the fastest method, whilst giving high-resolution data about the strains. In addition, CVTree agrees well with ANI-based method, revealing the relationship between B. thuringiensis and other B. cereus s.l. species. Based on these data, an online genome sequence comparison resource was built for Bacillus strains called the Bacillus Typing Bioinformatics Database to facilitate strain identification and characterization.

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.

Prevalence, attachment ability and strength of the biological control agent Bacillus thuringiensis on tomato

Bacillus thuringiensis (Bt) is commonly used as a biological control agent (BCA) to control insect pests in edible plant production and can as such be introduced into the food chain of fresh produce. Using standard food diagnostics Bt will be detected and reported as presumptive B. cereus. Tomato plants are often sprayed with Bt biopesticides for insect control, thus these Bt BCAs can also reach the tomato fruits and persist until consumption. In this study, vine tomatoes from the retail in Belgium (Flanders) were investigated for the occurrence and residual numbers of presumptive B. cereus and Bt. Of 109 tomato samples, 61 (56%) were tested positive for presumptive B. cereus. Of the presumptive B. cereus isolates (n = 213) recovered from these samples, 98% were identified as Bt by the production of parasporal crystals. Further quantitative real-time PCR assays on a subselection of Bt isolates (n = 61) showed that 95% of Bt isolates were indistinguishable from Bt biopesticide strains that are approved to be used on crops in the EU. Furthermore, the attachment strength of tested Bt biopesticide strains showed easier wash-off properties if using the commercial Bt granule formulation than the unformulated lab-cultured Bt or B. cereus spore suspensions.

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.

qPCR detection of viable Bacillus cereus group cells in cosmetic products

Reference methods for microbiological safety assessments of cosmetics rely on culture methods that reveal colonies of live microorganisms on growth media. Rapid molecular technologies, such as qPCR, detects the presence of target DNA in samples from dead and viable cells. DNA intercalating dyes, such as propidium monoazide (PMAxx), are capable of restricting PCR amplification to viable microbial cells. Here we developed singleplex and multiplex real time (qPCR) assays for the detection of Bacillus cereus (B. cereus) using 16S rRNA and phosphatidylcholine-specific phospholipase C (PLC) gene specific sequences coupled with PMAxx. The limit of detection was determined to be ~ 1 log CFU/ml for 16S rRNA and 3 log CFU/ml for PLC detection in pure culture using an eye shadow isolate, B. cereus 3A. We assessed the inclusivity and exclusivity of our qPCR assays using 212 strains, including 143 members of B. cereus, 38 non- B. cereus. and 31 non-Bacillus species; inclusivity was 100% for the 16S rRNA and 97.9% for the PLC targets; the exclusivity was 100% for 16S rRNA and 98.6% for PLC targets. These qPCR assays were then used to assess samples of commercial cosmetics: one set of liquid face toners (N = 3), artificially contaminated with B. cereus 3A, and one set of powdered cosmetics (N = 8), previously determined to be contaminated with B. cereus. For some samples, test portions were analyzed by qPCR in parallel, with and without PMAxx treatment. All test portions were simultaneously streaked on BACARA plates to confirm viable cells of B. cereus, according to the culture method. We found no difference in sensitivity between the singleplex and the multiplex qPCR assays (P > 0.05). Inoculated samples that did not recover B. cereus on plates still showed amplification of the DNA targets. However, that amplification was significantly delayed in PMAxx -treated samples (P < 0.0001) with C_T value differences of 7.82 for 16S rRNA and 7.22 for PLC. Likewise, amplification delay was significant (P < 0.0001) with inoculated samples that recovered B. cereus on plates with C_T value differences of 2.96 and 2.36 for 16S rRNA and PLC, respectively, demonstrating the presence of dead cells in the samples. All our qPCR results correlated with detection on BACARA plates (kappa, k = 0.99), independently of the presence of PMAxx in the PCR assays. Nevertheless, the amplification threshold with PMAxx dyes was significantly higher than the non-PMAxx dyes. Our findings confirm qPCR can be used for more rapid detection of microorganisms in cosmetics, including B. cereus, and selective detection of viable cells can be improved using PMAxx dyes.

In Search of Proximate Triggers of Anthrax Outbreaks in Wildlife: A Hypothetical Individual-Based Model of Plasmid Transfer within Bacillus Communities

Bacillus anthracis, the causative agent of anthrax in humans, livestock, and wildlife, exists in a community with hundreds of other species of bacteria in the environment. Work on the genetics of these communities has shown that B. anthracis shares a high percentage of chromosomal genes with both B. thuringiensis and B. cereus, and that phenotypic differences among these bacteria can result from extra-chromosomal DNA in the form of plasmids. We developed a simple hypothetical individual-based model to simulate the likelihood of detecting plasmids with genes encoding anthrax toxins within bacterial communities composed of B. anthracis, B. thuringiensis, and B. cereus, and the surrounding matrix of extra-cellular polymeric substances. Simulation results suggest the horizontal transfer of plasmids with genes encoding anthrax toxins among Bacillus species persisting outside the host could function as a proximate factor triggering anthrax outbreaks.

Genetic determinants of host tropism in Klebsiella phages

Bacteriophages play key roles in bacterial ecology and evolution and are potential antimicrobials. However, the determinants of phage-host specificity remain elusive. Here, we isolate 46 phages to challenge 138 representative clinical isolates of Klebsiella pneumoniae, a widespread opportunistic pathogen. Spot tests show a narrow host range for most phages, with <2% of 6,319 phage-host combinations tested yielding detectable interactions. Bacterial capsule diversity is the main factor restricting phage host range. Consequently, phage-encoded depolymerases are key determinants of host tropism, and depolymerase sequence types are associated with the ability to infect specific capsular types across phage families. However, all phages with a broader host range found do not encode canonical depolymerases, suggesting alternative modes of entry. These findings expand our knowledge of the complex interactions between bacteria and their viruses and point out the feasibility of predicting the first steps of phage infection using bacterial and phage genome sequences.

Soil Collected from a Single Great Smoky Mountains Trail Contains a Diversity of Listeria monocytogenes and Listeria spp

Listeria monocytogenes, a foodborne pathogen, and other Listeria spp. are present in natural environments. Isolating and characterizing strains from natural reservoirs can provide insight into the prevalence and diversity of Listeria spp. in these environments, elucidate their contribution to contamination of agricultural and food processing environments and food products, and lead to the discovery of novel species. In this study, we evaluated the diversity of Listeria spp. isolated from soil in a small region of the Great Smoky Mountains National Park, the most biodiverse national park in the U.S. National Park system. Of the 17 Listeria isolates recovered, whole-genome sequencing revealed that 14 were distinct strains. The strains represented a diversity of Listeria species (L. monocytogenes [n = 9], L. cossartiae subsp. cossartiae [n = 1], L. marthii [n = 1], L. booriae [n = 1], and a potentially novel Listeria sp. [n = 2]), as well as a diversity of sequence types based on multilocus sequence typing (MLST) and core genome MLST, including many novel designations. The isolates were not closely related (≥99.99% average nucleotide identity) to any isolates in public databases (NCBI, PATRIC), which also indicated novelty. The Listeria samples isolated in this study were collected from high-elevation sites near a creek that ultimately leads to the Mississippi River; thus, Listeria present in this natural environment could potentially travel downstream to a large region that includes portions of nine southeastern and midwestern U.S. states. This study provides insight into the diversity of Listeria spp. in the Great Smoky Mountains and indicates that this environment is a reservoir of novel Listeria spp. IMPORTANCE Listeria monocytogenes is a foodborne pathogen that can cause serious systemic illness that, although rare, usually results in hospitalization and has a relatively high mortality rate compared to other foodborne pathogens. Identification of novel and diverse Listeria spp. can provide insights into the genomic evolution, ecology, and evolution and variance of pathogenicity of this genus, especially in natural environments. Comparing L. monocytogenes and Listeria spp. isolates from natural environments, such as those recovered in this study, to contamination and/or outbreak strains may provide more information about the original natural sources of these strains and the pathways and mechanisms that lead to contamination of food products and agricultural or food processing environments.

Diversity of the genus Cryobacterium and proposal of 19 novel species isolated from glaciers

The bacterial genus Cryobacterium includes at present 14 species that live in cryospheric environments. In this study, we analyzed 101 genomes of Cryobacterium with pure cultures obtained from GenBank. They could be classified into 44 species based on average nucleotide identity (ANI) analysis, showing the diversity of Cryobacterium. Among these, 19 strains in our laboratory were isolated from the glacier samples in China. The pairwise ANI values of these 19 strains and known species were <95%, indicating that they represented 19 novel species. The comparative genomic analysis showed significant differences in gene content between the two groups with a maximum growth temperature (T _max) of ≤ 20°C and a T _max of >20°C. A comprehensive and robust phylogenetic tree, including 14 known species and 19 novel species, was constructed and showed five phylogenetic branches based on 265 concatenated single-copy gene sequences. The T _max parameter had a strong phylogenetic signal, indicating that the temperature adaptation of Cryobacterium was largely through vertical transfer rather than horizontal gene transfer and was affected by selection. Furthermore, using polyphasic taxonomy combined with phylogenomic analysis, we proposed 19 novel species of the genus Cryobacterium by the following 19 names: Cryobacterium serini sp. nov., Cryobacterium lactosi sp. nov., Cryobacterium gelidum sp. nov., Cryobacterium suzukii sp. nov., Cryobacterium fucosi sp. nov., Cryobacterium frigoriphilum sp. nov., Cryobacterium cryoconiti sp. nov., Cryobacterium lyxosi sp. nov., Cryobacterium sinapicolor sp. nov., Cryobacterium sandaracinum sp. nov., Cryobacterium cheniae sp. nov., Cryobacterium shii sp. nov., Cryobacterium glucosi sp. nov., Cryobacterium algoritolerans sp. nov., Cryobacterium mannosilyticum sp. nov., Cryobacterium adonitolivorans sp. nov., Cryobacterium algoricola sp. nov., Cryobacterium tagatosivorans sp. nov., and Cryobacterium glaciale sp. nov. Overall, the taxonomy and genomic analysis can improve our knowledge of phenotypic diversity, genetic diversity, and evolutionary characteristics of Cryobacterium.

Identification of Genetic Markers for the Detection of Bacillus thuringiensis Strains of Interest for Food Safety

Bacillus thuringiensis (Bt), belonging to the Bacillus cereus (Bc) group, is commonly used as a biopesticide worldwide due to its ability to produce insecticidal crystals during sporulation. The use of Bt, especially subspecies aizawai and kurstaki, to control pests such as Lepidoptera, generally involves spraying mixtures containing spores and crystals on crops intended for human consumption. Recent studies have suggested that the consumption of commercial Bt strains may be responsible for foodborne outbreaks (FBOs). However, its genetic proximity to Bc strains has hindered the development of routine tests to discriminate Bt from other Bc, especially Bacillus cereus sensu stricto (Bc ss), well known for its involvement in FBOs. Here, to develop tools for the detection and the discrimination of Bt in food, we carried out a genome-wide association study (GWAS) on 286 complete genomes of Bc group strains to identify and validate in silico new molecular markers specific to different Bt subtypes. The analyses led to the determination and the in silico validation of 128 molecular markers specific to Bt, its subspecies aizawai, kurstaki and four previously described proximity clusters associated with these subspecies. We developed a command line tool based on a 14-marker workflow, to carry out a computational search for Bt-related markers from a putative Bc genome, thereby facilitating the detection of Bt of interest for food safety, especially in the context of FBOs.

Identification and characterization of Bacillus thuringiensis and other Bacillus cereus group isolates from spinach by whole genome sequencing

Bacillus thuringiensis (Bt), used as a biological control agent (BCA), can persist on plants, and from there can be introduced into the final food product. In routine food safety diagnostics, these Bt residues cannot be distinguished from natural populations of Bacillus cereus present in plants and all are enumerated as "presumptive B. cereus." In this study, information on eventual use of Bt biopesticides, brand, application times and intervals provided by three food processing companies in Belgium, were integrated with quantitative data on presumptive B. cereus measured from fresh to frozen food products. This information together with data on genomic similarity obtained via whole genome sequencing (WGS) and cry gene profiling using a quantitative real-time PCR (qPCR) assay, confirmed that six out of 11 Bt isolates originated from the applied Bt biocontrol products. These identified Bt strains were shown to carry enterotoxin genes (nhe, hbl, cytK-2) and express Hbl enterotoxin in vitro. It was also noted that these Bt biopesticide strains showed no growth at standard refrigeration temperatures and a low or moderate biofilm-forming ability and cytotoxic activity. Our results also showed that the use of Bt as a BCA on spinach plants in the field led to higher residual counts of Bt in spinach (fresh or frozen) in the food supply chain, but the residual counts exceeding at present commonly assumed safety limit of 10⁵ CFU/g was only found in one fresh spinach sample. It is therefore recommended to establish a pre-harvest interval for Bt biopesticide application in the field to lower the likelihood of noncompliance to the generic B. cereus safety limit. Furthermore, WGS was found to be the best way to identify Bt biopesticide isolates at the strain level for foodborne outbreaks and clinical surveillance. The developed qPCR assay for screening on the presence of cry genes in presumptive B. cereus can be applied as a rapid routine test as an amendment to the already existing test on Bt crystal proteins determined via phase-contrast microscopy.

Gut Microbiome and Its Cofactors Are Linked to Lipoprotein Distribution Profiles

Increasing evidence indicates that the gut microbiome (GM) plays an important role in dyslipidemia. To date, however, no in-depth characterization of the associations between GM with lipoproteins distributions (LPD) among adult individuals with diverse BMI has been conducted. To determine such associations, we studied blood-plasma LPD, fecal short-chain fatty acids (SCFA) and GM of 262 Danes aged 19–89 years. Stratification of LPD segregated subjects into three clusters displaying recommended levels of lipoproteins and explained by age and body-mass-index. Higher levels of HDL2a and HDL2b were associated with a higher abundance of Ruminococcaceae and Christensenellaceae. Increasing levels of total cholesterol and LDL-1 and LDL-2 were positively associated with Lachnospiraceae and Coriobacteriaceae, and negatively with Bacteroidaceae and Bifidobacteriaceae. Metagenome-sequencing showed a higher abundance of biosynthesis of multiple B-vitamins and SCFA metabolism genes among healthier LPD profiles. Metagenomic-assembled genomes (MAGs) affiliated to Eggerthellaceae and Clostridiales were contributors of these genes and their relative abundance correlated positively with larger HDL subfractions. The study demonstrates that differences in composition and metabolic traits of the GM are associated with variations in LPD among the recruited subjects. These findings provide evidence for GM considerations in future research aiming to shed light on mechanisms of the GM–dyslipidemia axis.

Pathogenic Bacilli as an Emerging Biothreat?

Bacillus anthracis, present as a very durable endospore in soil, causes zoonotic illness which is mainly associated with herbivores and domestic animals. Human cases are scarce and often involve populations close to infected livestock. If anthrax is no longer of public health concern in developed countries, B. anthracis is one of the top-tier biological weapon agents. It is classified by the CDC as a category A agent. Since 1994, emerging strains of Bacillus cereus have been associated with anthrax-like disease in mammals. Some clinical strains of B. cereus harbor anthrax-like plasmid genes (pXO1 and pXO2) associated with non-human primate and human infections, with the same clinical presentation of inhalation anthrax and mortality rates. Although currently restricted to certain limited areas of circulation, the emergence of these new strains of B. cereus extends the list of potential agents possibly usable for bioterrorism or as a biological weapon. It is therefore important to improve our knowledge of the phylogeny within the B. cereus sensu lato group to better understand the origin of these strains. We can then more efficiently monitor the emergence of new strains to better control the risk of infection and limit potentially malicious uses.

Pan-Genome Analysis Reveals Functional Divergences in Gut-Restricted Gilliamella and Snodgrassella

Gilliamella and Snodgrassella, members of core gut microbiota in corbiculate bees, have high species diversity and adaptability to a wide range of hosts. In this study, we performed species taxonomy and phylogenetic analysis for Gilliamella and Snodgrassella strains that we isolated in our laboratory, in combination with published whole-genome. Functional effects of accessory and unique genes were investigated by KEGG category and pathway annotation in pan-genome analysis. Consequently, in Gilliamella, we inferred the importance of carbohydrate metabolism, amino acid metabolism, membrane transport, energy metabolism, and metabolism of cofactors and vitamins in accessory or unique genes. The pathway mentioned above, plus infectious disease, lipid metabolism, nucleotide metabolism as well as replication and repair exert a pivotal role in accessory or unique genes of Snodgrassella. Further analysis revealed the existence of functional differentiation of accessory and unique genes among Apis-derived genomes and Bombus-derived genomes. We also identified eight and four biosynthetic gene clusters in all Gilliamella and Snodgrassella genomes, respectively. Our study provides a good insight to better understand how host heterogeneity influences the bacterial speciation and affects the versatility of the genome of the gut bacteria.

Laboratory Misidentifications Resulting from Taxonomic Changes to Bacillus cereus Group Species, 2018-2022

Whole-genome sequencing (WGS) is being applied increasingly to Bacillus cereus group species; however, misinterpretation of WGS results may have severe consequences. We report 3 cases, 1 of which was an outbreak, in which misinterpretation of B. cereus group WGS results hindered communication within public health and industrial laboratories.

Bacillus cereus in Dairy Products and Production Plants

Spore-forming Bacillus cereus is a common contaminant of dairy products. As the microorganism is widespread in the environment, it can contaminate milk at the time of milking, but it can also reach the dairy products in each phase of production, storage and ripening. Milk pasteurization treatment is not effective in reducing contamination and can instead act as an activator of spore germination, and a potential associated risk still exists with the consumption of some processed foods. Prevalences and concentrations of B. cereus in milk and dairy products are extremely variable worldwide: in pasteurized milk, prevalences from 2% to 65.3% were reported, with concentrations of up to 3 × 10⁵ cfu/g, whereas prevalences in cheeses ranged from 0 to 95%, with concentrations of up to 4.2 × 10⁶ cfu/g. Bacillus cereus is also well known to produce biofilms, a serious concern for the dairy industry, with up to 90% of spores that are resistant to cleaning and are easily transferred. As the contamination of raw materials is not completely avoidable, and the application of decontamination treatments is only possible for some ingredients and is limited by both commercial and regulatory reasons, it is clear that the correct application of hygienic procedures is extremely important in order to avoid and manage the circulation of B. cereus along the dairy supply chain. Future developments in interventions must consider the synergic application of different mild technologies to prevent biofilm formation and to remove or inactivate the microorganism on the equipment.

Toxigenic Genes, Pathogenic Potential and Antimicrobial Resistance of Bacillus cereus Group Isolated from Ice Cream and Characterized by Whole Genome Sequencing

Bacillus cereus is isolated from a variety of foods where it may cause food spoilage and/or food poisoning due to its toxigenic and pathogenic nature. In this study, we identified members of B. cereus groups in 65% of the ice cream samples analyzed, which were characterized based on multi locus variable number tandem repeats analysis (MLVA) and whole genome sequencing (WGS). The MLVA revealed that 36 strains showed different allelic profiles. Analyses of WGS data enabled the identification of three members of the B. cereus group: B. cereus sensu stricto, B. mosaicus and B. thuringiensis. Based on the multi locus sequence typing (MLST) scheme, the strains were classified in 27 sequence types (STs), including ST26 that causes food poisoning. Toxin genes’ detection revealed the presence of the genes encoding nonhemolytic enterotoxin (NHE), hemolysin BL (HBL), cytotoxin K (cytK) and cereulide (ces) in 100%, 44%, 42% and 8% of the strains, respectively. The identification of eleven antimicrobial resistance (AMR) genes predicted the resistance to five different antimicrobials, and the resistance to beta-lactam antibiotics was confirmed with a phenotypic antimicrobial test. Taken together, the results showed that the B. cereus strains isolated from ice cream were a potential hazard for consumer safety.

Simultaneous Discrimination of Cereulide-Producing Bacillus cereus and Psychrotolerant B. cereus Group by Matrix-Assisted Laser Desorption Ionization-Time-of-Flight Mass Spectrometry

ABSTRACT

Cereulide-producing Bacillus cereus, which causes foodborne illnesses with vomiting, and psychrotolerant B. cereus group strains such as Bacillus mycoides, which can grow at ≥7°C and cause spoilage of refrigerated foods, are significant concerns for the food industry. Rapid and simple methods to discriminate the cereulide-producing B. cereus and psychrotolerant B. cereus group strains from other B. cereus group strains are needed. We developed a novel, rapid, and simple method with matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MALDI-TOF MS) analysis for simultaneous discrimination of these two groups from other B. cereus group strains. A potassium adduct of cereulide was used to detect cereulide-producing B. cereus, and three ribosomal subunit proteins (L30, S16, and S20) were used to detect psychrotolerant B. cereus group. A total of 51 B. cereus group strains were analyzed by MALDI-TOF MS. The biomarkers allowed successful discrimination of 16 cereulide-producing B. cereus and 15 psychrotolerant B. cereus group strains from other B. cereus group strains. The results showed that this MALDI-TOF MS analysis allows simultaneous discrimination of cereulide-producing B. cereus and psychrotolerant B. cereus group strains from other B. cereus group strains. This efficient method has the potential to be a valuable tool for ensuring food safety.

HIGHLIGHTS

Strains Associated with Two 2020 Welder Anthrax Cases in the United States Belong to Separate Lineages within Bacillus cereus sensu lato

Anthrax-causing members of Bacillus cereus sensu lato (s.l.) pose a serious threat to public health. While most anthrax-causing strains resemble B. anthracis phenotypically, rare cases of anthrax-like illness caused by strains resembling “B. cereus” have been reported. Here, whole-genome sequencing was used to characterize three B. cereus s.l. isolates associated with two 2020 welder anthrax cases in the United States, which resembled “B. cereus” phenotypically. Comparison of the three genomes sequenced here to all publicly available, high-quality B. cereus s.l. genomes (n = 2890 total genomes) demonstrated that genomes associated with each case effectively belonged to separate species at the conventional 95% average nucleotide identity prokaryotic species threshold. Two PubMLST sequence type 78 (ST78) genomes affiliated with a case in Louisiana were most closely related to B. tropicus and possessed genes encoding the Bps exopolysaccharide capsule, as well as hemolysin BL (Hbl) and cytotoxin K (CytK). Comparatively, a ST108 genome associated with a case in Texas was most closely related to B. anthracis; however, like other anthrax-causing strains most closely related to B. anthracis, this genome did not possess Bps-, Hbl-, or CytK-encoding genes. Overall, results presented here provide insights into the evolution of anthrax-causing B. cereus s.l.

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.

Genomic Sequencing of Bacillus cereus Sensu Lato Strains Isolated from Meat and Poultry Products in South Africa Enables Inter- and Intranational Surveillance and Source Tracking

Members of the Bacillus cereus sensu lato species complex, also known as the B. cereus group, vary in their ability to cause illness but are frequently isolated from foods, including meat products; however, food safety surveillance efforts that use whole-genome sequencing (WGS) often neglect these potential pathogens. Here, we evaluate the surveillance and source tracking potential of WGS as applied to B. cereus sensu lato by (i) using WGS to characterize B. cereus sensu lato strains isolated during routine surveillance of meat products across South Africa (n = 25) and (ii) comparing the genomes sequenced here to all publicly available, high-quality B. cereus sensu lato genomes (n = 2,887 total genomes). Strains sequenced here were collected from meat products obtained from (i) retail outlets, processing plants, and butcheries across six South African provinces (n = 23) and (ii) imports held at port of entry (n = 2). The 25 strains sequenced here were partitioned into 15 lineages via in silico seven-gene multilocus sequence typing (MLST). While none of the South African B. cereus sensu lato strains sequenced here were identical to publicly available genomes, six MLST lineages contained multiple strains sequenced in this study, which were identical or nearly identical at the whole-genome scale (≤3 core single nucleotide polymorphisms). Five MLST lineages contained (nearly) identical genomes collected from two or three South African provinces; one MLST lineage contained nearly identical genomes from two countries (South Africa and the Netherlands), indicating that B. cereus sensu lato can spread intra- and internationally via foodstuffs. IMPORTANCE Nationwide foodborne pathogen surveillance programs that use high-resolution genomic methods have been shown to provide vast public health and economic benefits. However, Bacillus cereus sensu lato is often overlooked during large-scale routine WGS efforts. Thus, to our knowledge, no studies to date have evaluated the potential utility of WGS for B. cereus sensu lato surveillance and source tracking in foodstuffs. In this preliminary proof-of-concept study, we applied WGS to B. cereus sensu lato strains collected via South Africa's national surveillance program of domestic and imported meat products, and we provide strong evidence that B. cereus sensu lato can be disseminated intra- and internationally via the agro-food supply chain. Our results showcase that WGS has the potential to be used for source tracking of B. cereus sensu lato in foods, although future WGS and metadata collection efforts are needed to ensure that B. cereus sensu lato surveillance initiatives are on par with those of other foodborne pathogens.

Establishment and Validation of a Two-Step LAMP Assay for Detection of Bacillus cereus-Group Isolates in Food and Their Possibility of Non-haemolytic Enterotoxin Production

The closely related members of the Bacillus cereus-group can mainly only be differentiated by whole genome sequencing. Among them, there are potentially toxin-producing bacteria. When consumed with food, these can cause vomiting or diarrhea and abdominal cramps. To date, although no EU-wide threshold exists, a bacterial count of 10⁵ CFU/g can be regarded as critical. Specific and rapid detection of the bacteria is difficult due to their close relationship, and no loop-mediated isothermal amplification (LAMP) assay has been developed so far to detect potentially toxin-producing members of the B. cereus-group. Aim of this study was to develop a LAMP method to detect critical cell counts specifically and rapidly of potentially non-haemolytic enterotoxin (NHE)-producing cells of this group. A two-step LAMP assay was developed. First, the target sequence groEL was used to determine the representatives of the B. cereus-group. Second, since bacteria in which nheB is present are basically capable of producing enterotoxins, this gene was chosen for detection. The specificity of the developed assay was 100% for B. cereus-group isolates and 93.7% for the detection of nheB. The analytical sensitivity was 0.1 pg DNA/μl. Using simplified DNA extraction by boiling, cell-based sensitivity was determined. Targeting groEL and nheB, 11.35-27.05 CFU/reaction and 11.35-270.5 CFU/reaction were detectable, respectively. Artificially contaminated samples were investigated to prove the application in foods. Direct detection of the critical value of B. cereus-group cells was possible in 83.3% of samples and detecting the toxin-gene 50% thereof. After a 6-h incubation period, the detection rate increased to 100 and 91.7%, respectively. Additionally, 100 natively contaminated food samples were tested, also quantitatively and culturally. Samples with relevant contamination levels were reliably detected using groEL-LAMP. After a 6-h incubation period, isolates bearing the toxin gene nheB could also be reliably detected. In addition, colony material was boiled and used as a LAMP template for simple detection. Specificity for the B. cereus-group was 100 and 93.22% detecting nheB. The study demonstrated that screening of food samples with the groEL/nheB-LAMP assay can be performed within 1 day, making it possible to detect critical levels of potentially NHE-toxin-producing cells of the B. cereus-group.

Phenotypic characterization and analysis of complete genomes of two distinct strains of the proposed species "L. swaminathanii"

Recently, a new Listeria species, “Listeria swaminathanii”, was proposed. Here, we phenotypically and genotypically characterize two additional strains that were previously obtained from soil samples and compare the results to the type strain. Complete genomes for both strains were assembled from hybrid Illumina and Nanopore sequencing reads and annotated. Further genomic analysis including average nucleotide identity (ANI) and detection of mobile genetic elements and genes of interest (e.g., virulence-associated) were conducted. The strains showed 98.7–98.8% ANI with the type strain. The UTK C1-0015 genome contained a partial monocin locus and a plasmid, while the UTK C1-0024 genome contained a full monocin locus and a prophage. Phenotypic characterization consistent with those performed on the proposed type strain was conducted to assess consistency of phenotypes across a greater diversity of the proposed species (n = 3 instead of n = 1). Only a few findings were notably different from those of the type strain, such as catalase activity, glycerol metabolism, starch metabolism, and growth at 41 °C. This study further expands our understanding of this newly proposed sensu stricto Listeria species.

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.

Bacillus cereus Invasive Infections in Preterm Neonates: an Up-to-Date Review of the Literature

Bacillus cereus group species are widespread, Gram-positive, spore-forming environmental bacteria. B. cereus sensu stricto is one of the major causes of food poisoning worldwide. In high-risk individuals, such as preterm neonates, B. cereus infections can cause fatal infections. It is important to note that the phenotypic identification methods commonly used in clinical microbiology laboratories make no distinction between B. cereus sensu stricto and the other members of the group (Bacillus anthracis excluded). As a result, all the invasive infections attributed to B. cereus are not necessarily due to B. cereus sensu stricto but likely to other closely related species of the B. cereus group. Next-generation sequencing (NGS) should be used to characterize the whole genome of the strains belonging to the B. cereus group. This could confirm whether the strains involved in previously reported B. cereus invasive infections preferentially belong to formerly known or emerging individual species. Moreover, infections related to B. cereus group species have probably been overlooked, since their isolation in human bacteriological samples has for a long time been regarded as an environmental contaminant of the cultures. Recent studies have questioned the emergence or reemergence of B. cereus invasive infections in preterm infants. This review reports our current understanding of B. cereus infections in neonates, including taxonomical updates, microbiological characteristics, bacterial identification, clinical features, host-pathogen interactions, environmental sources of contamination, and antimicrobial resistance.

In-depth genome analysis of Bacillus sp. BH32, a salt stress-tolerant endophyte obtained from a halophyte in a semiarid region

Endophytic strains belonging to the Bacillus cereus group were isolated from the halophytes Atriplex halimus L. (Amaranthaceae) and Tamarix aphylla L. (Tamaricaceae) from costal and continental regions in Algeria. Based on their salt tolerance (up to 5%), the strains were tested for their ability to alleviate salt stress in tomato and wheat. Bacillus sp. strain BH32 showed the highest potential to reduce salinity stress (up to + 50% and + 58% of dry weight improvement, in tomato and wheat, respectively, compared to the control). To determine putative mechanisms involved in salt tolerance and plant growth promotion, the whole genome of Bacillus sp. BH32 was sequenced, annotated, and used for comparative genomics against the genomes of closely related strains. The pangenome of Bacillus sp. BH32 and its closest relative was further analyzed. The phylogenomic analyses confirmed its taxonomic position, a member of the Bacillus cereus group, with intergenomic distances (GBDP analysis) pinpointing to a new taxon (digital DNA-DNA hybridization, dDDH < 70%). Genome mining unveiled several genes involved in stress tolerance, production of anti-oxidants and genes involved in plant growth promotion as well as in the production of secondary metabolites. KEY POINTS : • Bacillus sp. BH32 and other bacterial endophytes were isolated from halophytes, to be tested on tomato and wheat and to limit salt stress adverse effects. • The strain with the highest potential was then studied at the genomic level to highlight numerous genes linked to plant growth promotion and stress tolerance. • Pangenome approaches suggest that the strain belongs to a new taxon within the Bacillus cereus group.

The Complete Genome Sequence of Bacillus thuringiensis Strain ATCC 13367 Features a cry-Containing Chromosome

We report the high-quality closed genome sequence of Bacillus thuringiensis ATCC 13367, which is used as a bioinsecticide. The genome features a 328-kbp plasmid and a 5.8-Mbp chromosome. The genome is atypical for the species since it lacks a plasmid pesticidal crystal protein (cry) gene and instead features a chromosomal cry gene like that of B. thuringiensis HER1410.

The Bacillus cereus Strain EC9 Primes the Plant Immune System for Superior Biocontrol of Fusarium oxysporum

Antibiosis is a key feature widely exploited to develop biofungicides based on the ability of biological control agents (BCAs) to produce fungitoxic compounds. A less recognised attribute of plant-associated beneficial microorganisms is their ability to stimulate the plant immune system, which may provide long-term, systemic self-protection against different types of pathogens. By using conventional antifungal in vitro screening coupled with in planta assays, we found antifungal and non-antifungal Bacillus strains that protected the ornamental plant Kalanchoe against the soil-borne pathogen Fusarium oxysporum in experimental and commercial production settings. Further examination of one antifungal and one non-antifungal strain indicated that high protection efficacy in planta did not correlate with antifungal activity in vitro. Whole-genome sequencing showed that the non-antifungal strain EC9 lacked the biosynthetic gene clusters associated with typical antimicrobial compounds. Instead, this bacterium triggers the expression of marker genes for the jasmonic and salicylic acid defence pathways, but only after pathogen challenge, indicating that this strain may protect Kalanchoe plants by priming immunity. We suggest that the stimulation of the plant immune system is a promising mode of action of BCAs for the development of novel biological crop protection products.