Genes Associated With Psychrotolerant Bacillus cereus Group 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.

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

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

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

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

The 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.

Paenibacillus odorifer, the Predominant Paenibacillus Species Isolated from Milk in the United States, Demonstrates Genetic and Phenotypic Conservation of Psychrotolerance but Clade-Associated Differences in Nitrogen Metabolic Pathways

Although Paenibacillus species isolates are frequently isolated from pasteurized fluid milk, the link between the genetic diversity and phenotypic characteristics of these isolates was not well understood, especially as some Bacillales isolated from milk are unable to grow at refrigeration temperatures. Our data demonstrate that Paenibacillus spp. isolated from fluid milk represent tremendous interspecies diversity, with P. odorifer being the predominant Paenibacillus sp. isolated. Furthermore, genetic and phenotypic data support that P. odorifer is well suited to transition from a soil-dwelling environment, where nitrogen fixation (and other nitrate/nitrite reduction pathways present only in clade A) may facilitate growth, to fluid milk, where its multiple cold shock-associated adaptations enable it to grow at refrigeration temperatures throughout the storage of milk. Therefore, efforts to reduce bacterial contamination of milk will require a systematic approach to reduce P. odorifer contamination of raw milk.

Paenibacillus is a spore-forming bacterial genus that is frequently isolated from fluid milk and is proposed to play a role in spoilage. To characterize the genetic and phenotypic diversity of Paenibacillus spp., we first used rpoB allelic typing data for a preexisting collection of 1,228 Paenibacillus species isolates collected from raw and processed milk, milk products, and dairy environmental sources. Whole-genome sequencing (WGS) and average nucleotide identity by BLAST (ANIb) analyses performed for a subset of 58 isolates representing unique and overrepresented rpoB allelic types in the collection revealed that these isolates represent 21 different Paenibacillus spp., with P. odorifer being the predominant species. Further genomic characterization of P. odorifer isolates identified two distinct phylogenetic clades, clades A and B, which showed significant overrepresentation of 172 and 164 ortholog clusters and 94 and 52 gene ontology (GO) terms, respectively. While nitrogen fixation genes were found in both clades, multiple genes associated with nitrate and nitrite reduction were overrepresented in clade A isolates; additional phenotypic testing demonstrated that nitrate reduction is specific to isolates in clade A. Hidden Markov models detected 9 to 10 different classes of cold shock-associated genetic elements in all P. odorifer isolates. Phenotypic testing revealed that all isolates tested here can grow in skim milk broth at 6°C, suggesting that psychrotolerance is conserved in P. odorifer. Overall, our data suggest that Paenibacillus spp. isolated from milk in the United States represent broad genetic diversity, which may provide challenges for targeted-control strategies aimed at reducing fluid milk spoilage.

IMPORTANCE Although Paenibacillus species isolates are frequently isolated from pasteurized fluid milk, the link between the genetic diversity and phenotypic characteristics of these isolates was not well understood, especially as some Bacillales isolated from milk are unable to grow at refrigeration temperatures. Our data demonstrate that Paenibacillus spp. isolated from fluid milk represent tremendous interspecies diversity, with P. odorifer being the predominant Paenibacillus sp. isolated. Furthermore, genetic and phenotypic data support that P. odorifer is well suited to transition from a soil-dwelling environment, where nitrogen fixation (and other nitrate/nitrite reduction pathways present only in clade A) may facilitate growth, to fluid milk, where its multiple cold shock-associated adaptations enable it to grow at refrigeration temperatures throughout the storage of milk. Therefore, efforts to reduce bacterial contamination of milk will require a systematic approach to reduce P. odorifer contamination of raw milk.