Genomic insights into the taxonomic status of the Bacillus cereus group

Prevalence, Characteristics, and Selection of Bacillus cereus Subgroups from Dairy Products for Challenge Testing and Predictive Model Development

Prevalence, toxin gene profiles, lactose fermentation, and growth responses of B. cereus sensu lato subgroups in various dairy and dairy alternative products and ingredients were studied to identify relevant isolates for challenge testing and model development to predict and manage growth responses. Out of 71 examined products or ingredients, 51 B. cereus s.l. isolates were obtained from 35 positive samples (49% prevalence). These 51 isolates along with 18 additional dairy isolates and 12 B. cereus s.l. reference strains were identified using MALDI-TOF. The 81 isolates were further characterized by panC sequencing, testing for cold shock and toxin genes (cspA; hbl, nhe, CytK and ces), lactose fermentation, and study of growth rates (µmax) under various conditions (45 °C, 10 °C, 6% NaCl, pH 5.1), resulting in 298 µmax-values. These conditions were selected to differentiate mesophilic and psychrotolerant strains and to identify tolerant isolates. Dairy powders (83%), pasteurized upconcentrated cheese whey (43%), and cheeses (42%) had the highest prevalences of B. cereus s.l. and the highest concentrations in positive samples (5-100 CFU/g or ml). The panC groups II, III, IV, VI, and VIII were detected among the dairy isolates, with 97% harboring one or more toxin genes. Lactose fermentation was observed in 42% of isolates, with lactose-fermenting B. cereus s.l. isolates of panC groups III and IV dominant in dairy powders. Growth rates of B. cereus s.l. varied considerably, among and within panC groups of the dairy isolates. Based on the highest growth rates at 45 °C, 10 °C, 6% NaCl, pH 5.1, panC group membership, toxin genes profiles, and lactose fermentation ability, two cocktails of "mesophilic" or "psychrotolerant" isolates were selected. These strain cocktails can be used in future challenge testing and predictive food microbiology studies to evaluate and manage the growth of B. cereus s.l. in dairy products and ingredients.

First insights into the prevalence, genetic characteristics, and pathogenicity of Bacillus cereus from generations worldwide

Bacillus cereus, a global threat, is one of the major causes of toxin-induced foodborne diseases. However, a comprehensive assessment of the prevalence and characteristics of B. cereus worldwide is still lacking. Here, we applied whole-genome sequence analysis to 191 B. cereus collected in Africa, America, Asia, Europe, and Oceania from the 1900s to 2022, finding that CC142 dominated the global B. cereus clonal complex. The results provided direct evidence that B. cereus could spread through the food chain and intercontinentally. B. cereus from different generations worldwide showed coherence in the antibiotic-resistant gene and virulence and biofilm gene profiles, although with high genomic heterogeneity. The BCI-BCII-vanZF-fosB profiles and virulence and biofilm genes were detected at high rates, and we emphasized that B. cereus would pose a serious challenge to global public health and clinical treatment.IMPORTANCEThis study first emphasized the prevalence, genetic characteristics, and pathogenicity of Bacillus cereus worldwide from the 1900s to 2022 using whole-genome sequence analysis. The CC142 dominated the global Bacillus cereus clonal complex. Moreover, we revealed a close evolutionary relationship between the isolates from different sources. B. cereus isolates from different generations worldwide showed coherence in potential pathogenicity, although with high genomic heterogeneity. The BCI-BCII-vanZF-fosB profiles and virulence and biofilm genes were detected at high rates, and we emphasized that B. cereus would pose a serious challenge to global public health and clinical treatment.

Characterization of the Bacillus cereus Group Isolated from Ready-to-Eat Foods in Poland by Whole-Genome Sequencing

Bacillus cereus sensu lato can contaminate food and cause food poisoning by producing toxins such as cereulide, toxin BL, and cytotoxin K. In this study, we retrospectively analyzed B. cereus sensu lato from retail food products and food poisoning cases using PCR methods to determine their virulence profiles. A new toxin profile, encoding all four toxins (hbl, nhe, cytK, ces), was found in 0.4% of isolates. The toxin profiles, classified into A-J, revealed that 91.8% harbored nhe genes, while hbl, cytK, and ces were detected in 43.8%, 46.9%, and 4.2% of isolates, respectively. Whole-genome sequencing (WGS) identified four distinct species within the B. cereus group, with 21 isolates closely related to B. cereus sensu stricte, 25 to B. mosaicus, 2 to B. toyonensis, and 1 to B. mycoides. Three novel sequence types (STs 3297, 3298, 3299) were discovered. Antibiotic resistance genes were common, with 100% of isolates carrying beta-lactam resistance genes. Fosfomycin (80%), vancomycin (8%), streptothricin (6%), tetracycline (4%), and macrolide resistance (2%) genes were also detected. These results highlight the genetic diversity and antibiotic resistance potential of B. cereus sensu lato strains in Polish food products.

Investigating Anthrax-Associated Virulence Genes among Archival and Contemporary Bacillus cereus Group Genomes

Bacillus anthracis causes anthrax through virulence factors encoded on two plasmids. However, non-B. anthracis organisms within the closely related, environmentally ubiquitous Bacillus cereus group (BCG) may cause an anthrax-like disease in humans through the partial adoption of anthrax-associated virulence genes, challenging the definition of anthrax disease. To elucidate these phenomena and their evolutionary past, we performed whole-genome sequencing on non-anthracis BCG isolates, including 93 archival (1967-2003) and 5 contemporary isolates (2019-2023). We produced annotated genomic assemblies and performed a pan-genome analysis to identify evidence of virulence gene homology and virulence gene acquisition by linear inheritance or horizontal gene transfer. At least one anthrax-associated virulence gene was annotated in ten isolates. Most homologous sequences in archival isolates showed evidence of pseudogenization and subsequent gene loss. The presence or absence of accessory genes, including anthrax-associated virulence genes, aligned with the phylogenetic structure of the BCG core genome. These findings support the hypothesis that anthrax-associated virulence genes were inherited from a common ancestor in the BCG and were retained or lost across different lineages, and contribute to a growing body of work informing public health strategies related to anthrax surveillance and identification.

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

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

Characterization, whole-genome sequence analysis, and protease production of a new thermophilic Bacillus licheniformis strain isolated from Debagh hot spring, Algeria

A new thermophilic strain, designated as Bacillus sp. LMB3902, was isolated from Hammam Debagh, the hottest spring in Algeria (up to 98 °C). This isolate showed high protease production in skim milk media at 55 °C and exhibited significant specific protease activity by using azocasein as a substrate (157.50 U/mg). Through conventional methods, chemotaxonomic characteristics, 16S rRNA gene sequencing, and comparative genomic analysis with the closely related strain Bacillus licheniformis DSM 13 (ATCC 14580 T), the isolate Bacillus sp. LMB3902 was identified as a potentially new strain of Bacillus licheniformis. In addition, the gene functions of Bacillus sp. LMB3902 strain were predicted using the Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Clusters of Orthologous Groups, Non-Redundant Protein Sequence Database, Swiss-Prot, and Pfam databases. The results showed that the genome size of Bacillus sp. LMB3902 was 4.279.557 bp, with an average GC content of 46%. The genome contained 4.760 predicted genes, including 8 rRNAs, 78 tRNAs, and 24 sRNAs. A total of 235 protease genes were annotated including 50 proteases with transmembrane helix structures and eight secreted proteases with signal peptides. Additionally, the majority of secondary metabolites found by antiSMASH platform showed low similarity to identified natural products, such as fengicin (53%), lichenysin (57%), and surfactin (34%), suggesting that this strain may encode for novel uncharacterized natural products which can be useful for biotechnological applications. This study is the first report that describes the complete genome sequence, taxono-genomics, and gene annotation as well as protease production of the Bacillus genus in this hydrothermal vent.

Biofilm-forming Ability of Bacillus thuringiensis Strains from Biopesticides on Polystyrene and their Attachment on Spinach

Bacillus thuringiensis-based commercial products as a biopesticide have been used for more than 60 years in agriculture. However, as one of the species in B. cereus group, B. thuringiensis has been considered as an emerging hazard with the potential to cause food toxico-infections. The present study aimed to evaluate the biofilm-forming ability of B. thuringiensis biopesticide strains and their attachment on spinach, compared to foodborne B. cereus strains. Biofilm formations of tested strains were found to be strain-specific and affected by the nutrient conditions more than the incubation time. Nutrient starvation conditions generally reduced the biofilm formation of tested B. thuringiensis and B. cereus strains, particularly B. thuringiensis ABTS-1857 strain was found as the nonbiofilm former in starvation conditions. It is worth mentioning that B. thuringiensis SA-11 strain showed stronger biofilm-forming ability with more air-liquid interface biofilm than the other two B. thuringiensis biopesticide strains, but no such higher attachment of B. thuringiensis SA-11 to spinach was observed. These results indicate that B. thuringiensis SA-11 strain can enter the food processing lines by the attachment on spinach leaves, and it has the potential to form biofilms throughout the processing lines or the production environment when sufficient nutrients are available. However, more biofilm tests of B. thuringiensis biopesticide strains in the vegetable production chain should be performed. The dry formulation of commercial B. thuringiensis biopesticides enhanced their adhesion on spinach leaves, whereas the strength of adhesion was not improved by the formulation. In addition, 1-2 log reductions of spores after the intensive washing of spinach leaves in the lab were detected. However, the log reduction due to the actual washing done by the food processing companies in large-volume washing baths or by consumers at home would be limited and less than this lab simulation.

Evaluating the Safety of Bacillus cereus GW-01 Obtained from Sheep Rumen Chyme

Bacillus cereus is responsible for 1.4-12% food poisoning outbreaks worldwide. The safety concerns associated with the applications of B. cereus in health and medicine have been controversial due to its dual role as a pathogen for foodborne diseases and a probiotic in humans and animals. In this study, the pathogenicity of B. cereus GW-01 was assessed by comparative genomic, and transcriptome analysis. Phylogenetic analysis based on a single-copy gene showed clustering of the strain GW-01, and 54 B. cereus strains from the NCBI were classified into six major groups (I-VI), which were then associated with the source region and sequence types (STs). Transcriptome results indicated that the expression of most genes related with toxins secretion in GW-01 was downregulated compared to that in the lag phase. Overall, these findings suggest that GW-01 is not directly associated with pathogenic Bacillus cereus and highlight an insightful strategy for assessing the safety of novel B. cereus strains.

Antimicrobial Resistance Profile, Whole-Genome Sequencing and Core Genome Multilocus Sequence Typing of B. anthracis Isolates in Croatia from 2001 to 2022

Bacillus anthracis, the causative agent of anthrax disease, is a worldwide threat to livestock, wildlife and public health. It is also considered one of the most important pathogens of bioterrorism. Rapid and reliable diagnosis and administration of antimicrobials are essential for effective anthrax treatment. In this study, we determined the in vitro susceptibilities of 40 isolates of B. anthracis isolated in Croatia over the recent two decades to 18 antimicrobials. Whole-genome sequencing was performed, and bioinformatics tools were used to determine virulence factors and antimicrobial resistance genes. Core genome-based multilocus sequence typing was used for isolate comparison and phylogenetic analysis. All isolates were susceptible to all antimicrobials recommended for post-exposure prophylaxis or anthrax therapy. Susceptibility was found to all other tested antimicrobials that are an alternative for primary therapy. We found two beta-lactamase genes, but their expression is not sufficient to confer resistance. In all isolates used in this study, we found 21 virulence genes, 8 of which are responsible for toxin and capsule production. As far as phylogenetic analysis is concerned, the B. anthracis isolates from Croatia are categorised into two clades. The first is clade A, subclade Trans Eurasia, and the other is clade B, subclade B2.

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

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

Microbiological evaluation of the indigenous fermented condiment okpeye available at various retail markets in the south-eastern region of Nigeria

In Africa, indigenous fermented condiments contribute to food security as a low-cost source of protein. Okpeye is an indigenous fermented condiment produced from Prosopis africana seeds. The reliance on spontaneous fermentation processes and unhygienic practices during production often results in the contamination of the final product with microbial hazards. A microbiological evaluation of 18 commercial samples of okpeye purchased from six markets in two cities in southeastern Nigeria was conducted. Fifty-nine (59) bacteria were isolated and identified at the species level by phenotyping and sequencing the 16S rRNA, gyrB and rpoB genes. Bacillus (47.4 %) and Staphylococcus (42.3 %) were the predominant bacterial genera in okpeye. Overall, B. amyloliquefaciens and S. simulans were the most frequently occurring bacteria and were present in all samples. In addition, B. cereus was isolated in samples obtained from all markets. Other bacterial species included B. velezensis, Oceanobacillus caeni, S. cohnii, Escherichia fergusonni and Vagacoccus lutrae. The B. cereus isolates (10) were screened for the presence of 8 enterotoxin genes (hblA, hblC, hblD, nheA, nheB, nheC, cytK, entFM) and one emetic gene (cesB). The non-haemolytic enterotoxin (nheABC) and haemolytic enterotoxin (hblABD) complexes were present in 70 % and 50 % of B. cereus respectively. The positive rate of cytK and entFM genes was 70 %, while the cesB gene was 30 %. Antibiotic susceptibility assessment showed that most of the isolates were susceptible to gentamicin, tetracycline, streptomycin, and erythromycin but resistant to ciprofloxacin and vancomycin. These findings highlight the need for further controls to reduce contamination with potential pathogenic bacteria in indigenous fermented condiments such as okpeye. There is also a need to educate producers regarding hygienic practices to safeguard public health and food security.

Clostridium strain FAM25158, a unique endospore-forming bacterium related to Clostridium tyrobutyricum and isolated from Emmental cheese shows low tolerance to salt

The genus Clostridium is a large and diverse group of species that can cause food spoilage, including late blowing defect (LBD) in cheese. In this study, we investigated the taxonomic status of strain FAM25158 isolated from Emmental cheese with LBD using a polyphasic taxonomic and comparative genomic approach. A 16S rRNA gene sequence phylogeny suggested affiliation to the Clostridium sensu stricto cluster, with Clostridium tyrobutyricum DSM 2637T being the closest related type strain (99.16% sequence similarity). Average Nucleotide Identity (ANI) analysis revealed that strain FAM25158 is at the species threshold with C. tyrobutyricum, with ANI values ranging from 94.70 to 95.26%, while the digital DNA-DNA hybridization values were below the recommended threshold, suggesting that FAM25158 is significantly different from C. tyrobutyricum at the genomic level. Moreover, comparative genomic analysis between FAM25158 and its four closest C. tyrobutyricum relatives revealed a diversity of metabolic pathways, with FAM25158 differing from other C. tyrobutyricum strains by the presence of genes such as scrA, srcB, and scrK, responsible for sucrose utilization, and the absence of many important functional genes associated with cold and osmolality adaptation, which was further supported by phenotypic analyses. Surprisingly, strain FAM25158 exhibited unique physiologic traits, such as an optimal growth temperature of 30°C, in contrast to its closest relatives, C. tyrobutyricum species with an optimal growth temperature of 37°C. Additionally, the growth of FAM25158 was inhibited at NaCl concentrations higher than 0.5%, a remarkable observation considering its origin from cheese. While the results of this study provide novel information on the genetic content of strain FAM25158, the relationship between its genetic content and the observed phenotype remains a topic requiring further investigation.

Organization, conservation, and diversity of biosynthetic gene clusters in Bacillus sp. BH32 and its closest relatives in the Bacillus cereus group

This study explores the organization, conservation, and diversity of biosynthetic gene clusters (BGCs) among Bacillus sp. strain BH32, a plant-beneficial bacterial endophyte, and its closest nontype Bacillus cereus group strains. BGC profiles were predicted for each of the 17 selected strains using antiSMASH, resulting in the detection of a total of 198 BGCs. We quantitatively compared the BGCs and analysed their conservation, distribution, and evolutionary relationships. The study identified both conserved and singleton BGCs across the studied Bacillus strains, with minimal variation, and discovered two major BGC synteny blocks composed of homologous BGCs conserved within the B. cereus group. The identified BGC synteny blocks provide insight into the evolutionary relationships and diversity of BGCs within this complex group.

De novo assembly and comparative genome analysis for polyhydroxyalkanoates-producing Bacillus sp. BNPI-92 strain

BACKGROUND

Certain Bacillus species play a vital role in polyhydroxyalkanoate (PHA) production. However, most of these isolates did not properly identify to species level when scientifically had been reported.

RESULTS

From NGS analysis, 5719 genes were predicted in the de novo genome assembly. Based on genome annotation using RAST server, 5,527,513 bp sequences were predicted with 5679 bp number of protein-coding sequence. Its genome sequence contains 35.1% and 156 GC content and contigs, respectively. In RAST server analysis, subsystem (43%) and non-subsystem coverage (57%) were generated. Ortho Venn comparative genome analysis indicated that Bacillus sp. BNPI-92 shared 2930 gene cluster (core gene) with B. cereus ATCC 14579 T (AE016877), B. paranthracis Mn5T (MACE01000012), B. thuringiensis ATCC 10792 T (ACNF01000156), and B. antrics Amen T (AE016879) strains. For our strain, the maximum gene cluster (190) was shared with B. cereus ATCC 14579 T (AE016877). For Ortho Venn pair wise analysis, the maximum overlapping gene clusters thresholds have been detected between Bacillus s p.BNPI-92 and Ba. cereus ATCC 14579 T (5414). Average nucleotide identity (ANI) such as OriginalANI and OrthoANI, in silicon digital DND-DNA hybridization (isDDH), Type (Strain) Genome Server (TYGS), and Genome-Genome Distance Calculator (GGDC) were more essentially related Bacillus sp. BNPI-92 with B. cereus ATCC 14579 T strain. Therefore, based on the combination of RAST annotation, OrthoVenn server, ANI and isDDH result Bacillus sp.BNPI-92 strain was strongly confirmed to be a B. cereus type strain. It was designated as B. cereus BNPI-92 strain. In B. cereus BNPI-92 strain whole genome sequence, PHA biosynthesis encoding genes such as phaP, phaQ, phaR (PHA synthesis repressor phaR gene sequence), phaB/phbB, and phaC were predicted on the same operon. These gene clusters were designated as phaPQRBC. However, phaA was located on other operons.

CONCLUSIONS

This newly obtained isolate was found to be new a strain based on comparative genomic analysis and it was also observed as a potential candidate for PHA biosynthesis.

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.

Plant-Associated Representatives of the Bacillus cereus Group Are a Rich Source of Antimicrobial Compounds

Seventeen bacterial strains able to suppress plant pathogens have been isolated from healthy Vietnamese crop plants and taxonomically assigned as members of the Bacillus cereus group. In order to prove their potential as biocontrol agents, we perform a comprehensive analysis that included the whole-genome sequencing of selected strains and the mining for genes and gene clusters involved in the synthesis of endo- and exotoxins and secondary metabolites, such as antimicrobial peptides (AMPs). Kurstakin, thumolycin, and other AMPs were detected and characterized by different mass spectrometric methods, such as MALDI-TOF-MS and LIFT-MALDI-TOF/TOF fragment analysis. Based on their whole-genome sequences, the plant-associated isolates were assigned to the following species and subspecies: B. cereus subsp. cereus (6), B. cereus subsp. bombysepticus (5), Bacillus tropicus (2), and Bacillus pacificus. These three isolates represent novel genomospecies. Genes encoding entomopathogenic crystal and vegetative proteins were detected in B. cereus subsp. bombysepticus TK1. The in vitro assays revealed that many plant-associated isolates enhanced plant growth and suppressed plant pathogens. Our findings indicate that the plant-associated representatives of the B. cereus group are a rich source of putative antimicrobial compounds with potential in sustainable agriculture. However, the presence of virulence genes might restrict their application as biologicals in agriculture.

Bacillus dicomae sp. nov., a new member of the Bacillus cereus group isolated from medicinal plant Dicoma anomala

A Gram-stain-positive, endospore-forming endophytic bacterial strain designated MHSD28T was isolated from surface-sterilized leaves of Dicoma anomala collected from Eisleben, Botlokwa, Limpopo Province, South Africa. The phenotypic and phylogenetic characteristics of strain MHSD28T were consistent with those of members within the Bacillus cereus group. Comparative analysis between this strain and its relatives confirmed that it belongs to this group and forms a monophyletic branch. The digital DNA-DNA hybridization values between strain MHSD28T and its relatives were lower than the 70 % threshold for species delineation. To further determine its phylogenetic position, multi-locus sequence analysis (MLSA) based on five concatenated housekeeping gene (gyrB, atpD, DnaK, rpoB and rpoD) sequences, phenotypic analysis, matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) biotyper identification, fatty acid and polar lipid profile analyses were carried out. Phenotypic characterization, MLSA, whole genome sequence based analyses and MALDI-TOF results placed strain MHSD28T within the B. cereus group. The major fatty acids were iso-C15 : 0 and summed feature 3 and the main polar lipids were diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine. The respiratory quinone was menaquinone-7. The cell-wall peptidoglycan structure included meso-diaminopimelic acid. Considering the above results, strain MHSD28T represents a novel species of the B. cereus group, for which the name Bacillus dicomae sp. nov. is proposed. The type strain is MHSD28T (=BD 2262T=LMG 32287T=CECT 30671T).

Elucidation of the Bovine Intramammary Bacteriome and Resistome from healthy cows of Swiss dairy farms in the Canton Tessin

Healthy, untreated cows of nine dairy herds from the Swiss Canton Tessin were analyzed three times within one year to identify the most abundant species of the intramammary bacteriome. Aseptically collected milk samples were cultured and bacteria identified using MALDI-TOF. Of 256 cows analyzed, 96% were bacteriologically positive and 80% of the 1,024 quarters were positive for at least one bacterial species. 84.5% of the quarters were healthy with somatic cell counts (SCC) < 200,000 cells/mL, whereas 15.5% of the quarters showed a subclinical mastitis (SCC ≥ 200,000 cells/mL). We could assign 1,288 isolates to 104 different bacterial species including 23 predominant species. Non-aureus staphylococci and mammaliicocci (NASM) were most prevalent (14 different species; 73.5% quarters). Staphylococcus xylosus and Mammaliicoccus sciuri accounted for 74.7% of all NASM isolates. To describe the intramammary resistome, 350 isolates of the predominant species were selected and subjected to short-read whole genome sequencing (WGS) and phenotypic antibiotic resistance profiling. While complete genomes of eight type strains were available, the remaining 15 were de novo assembled with long reads as a resource for the community. The 23 complete genomes served for reference-based assembly of the Illumina WGS data. Both chromosomes and mobile genetic elements were examined for antibiotic resistance genes (ARGs) using in-house and online software tools. ARGs were then correlated with phenotypic antibiotic resistance data from minimum inhibitory concentration (MIC). Phenotypic and genomic antimicrobial resistance was isolate-specific. Resistance to clindamycin and oxacillin was most frequently observed (65 and 30%) in Staphylococcus xylosus but could not be linked to chromosomal or plasmid-borne ARGs. However, in several cases, the observed antimicrobial resistance could be explained by the presence of mobile genetic elements like tetK carried on small plasmids. This represents a possible mechanism of transfer between non-pathogenic bacteria and pathogens of the mammary gland within and between herds. The-to our knowledge-most extensive bacteriome reported and the first attempt to link it with the resistome promise to profoundly affect veterinary bacteriology in the future and are highly relevant in a One Health context, in particular for mastitis, the treatment of which still heavily relies on antibiotics.

Retrospective surveillance of viable Bacillus cereus group contaminations in commercial food and feed vitamin B2 products sold on the Belgian market using whole-genome sequencing

Bacillus cereus is a spore-forming bacterium that occurs as a contaminant in food and feed, occasionally resulting in food poisoning through the production of various toxins. In this study, we retrospectively characterized viable B. cereus sensu lato (s.l.) isolates originating from commercial vitamin B₂ feed and food additives collected between 2016 and 2022 by the Belgian Federal Agency for the Safety of the Food Chain from products sold on the Belgian market. In total, 75 collected product samples were cultured on a general medium and, in case of bacterial growth, two isolates per product sample were collected and characterized using whole-genome sequencing (WGS) and subsequently characterized in terms of sequence type (ST), virulence gene profile, antimicrobial resistance (AMR) gene profile, plasmid content, and phylogenomic relationships. Viable B. cereus was identified in 18 of the 75 (24%) tested products, resulting in 36 WGS datasets, which were classified into eleven different STs, with ST165 (n = 10) and ST32 (n = 8) being the most common. All isolates carried multiple genes encoding virulence factors, including cytotoxin K-2 (52.78%) and cereulide (22.22%). Most isolates were predicted to be resistant to beta-lactam antibiotics (100%) and fosfomycin (88.89%), and a subset was predicted to be resistant to streptothricin (30.56%). Phylogenomic analysis revealed that some isolates obtained from different products were closely related or even identical indicating a likely common origin, whereas for some products the two isolates obtained did not show any close relationship to each other or other isolates found in other products. This study reveals that potentially pathogenic and drug-resistant B. cereus s.l. can be present in food and feed vitamin B₂ additives that are commercially available, and that more research is warranted to assess whether their presence in these types of products poses a threat to consumers.

Retrospective Screening of Anthrax-like Disease Induced by Bacillus tropicus str. JMT from Chinese Soft-Shell Turtles in Taiwan

Bacillus cereus is ubiquitous in the environment and a well-known causative agent of foodborne disease. Surprisingly, more and more emerging strains of atypical B. cereus have been identified and related to severe disease in humans and mammals such as chimpanzees, apes, and bovine. Recently, the atypical B. cereus isolates, which mainly derive from North America and Africa, have drawn great attention due to the potential risk of zoonosis. The cluster of B. cereus carries several anthrax-like virulent genes that are implicated in lethal disease. However, in non-mammals, the distribution of atypical B. cereus is still unknown. In this study, we conducted a retrospective screening of the 32 isolates of Bacillus spp. from diseased Chinese soft-shelled turtles from 2016 to 2020. To recognize the causative agent, we used various methods, such as sequencing analysis using PCR-amplification of the 16S rRNA gene, multiplex PCR for discriminating, and colony morphology by following previous studies. Furthermore, the digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values were calculated, respectively, below the 70 and 96% cutoff to define species boundaries. According to the summarized results, the pathogen is taxonomically classified as Bacillus tropicus str. JMT (previous atypical Bacillus cereus). Subsequently, analyses such as targeting the unique genes using PCR and visual observation of the bacteria under various staining techniques were implemented in our study. Our findings show that all (32/32, 100%) isolates in this retrospective screening share similar phenotypical properties and carry the protective antigen (PA), edema factor (EF), hyaluronic acid (HA), and exopolysaccharide (Bps) genes on their plasmids. In this study, the results indicate that the geographic distribution and host range of B. tropicus were previously underestimated.

DNA Nanomachine (DNM) Biplex Assay for Differentiating Bacillus cereus Species

Conventional methods for the detection and differentiation of Bacillus cereus group species have drawbacks mostly due to the complexity of genetic discrimination between the Bacillus cereus species. Here, we describe a simple and straightforward assay based on the detected unamplified bacterial 16S rRNA by DNA nanomachine (DNM). The assay uses a universal fluorescent reporter and four all-DNA binding fragments, three of which are responsible for "opening up" the folded rRNA while the fourth stand is responsible for detecting single nucleotide variation (SNV) with high selectivity. Binding of the DNM to 16S rRNA results in the formation of the 10-23 deoxyribozyme catalytic core that cleaves the fluorescent reporter and produces a signal, which is amplified over time due to catalytic turnover. This developed biplex assay enables the detection of B. thuringiensis 16S rRNA at fluorescein and B. mycoides at Cy5 channels with a limit of detection of 30 × 10³ and 35 × 10³ CFU/mL, respectively, after 1.5 h with a hands-on time of ~10 min. The new assay may simplify the analysis of biological RNA samples and might be useful for environmental monitoring as a simple and inexpensive alternative to amplification-based nucleic acid analysis. The DNM proposed here may become an advantageous tool for detecting SNV in clinically significant DNA or RNA samples and can easily differentiate SNV under broadly variable experimental conditions and without prior amplification.

Elucidating the genomic history of commercially used Bacillus thuringiensis subsp. tenebrionis strain NB176

Bacillus thuringiensis subsp. tenebrionis (Btt) produces a coleopteran-specific crystal protoxin protein (Cry3Aa δ-endotoxin). After its discovery in 1982, the strain NB125 (DSM 5526) was eventually registered in 1990 to control the Colorado potato beetle (Leptinotarsa decemlineata). Gamma-irradiation of NB125 resulted in strain NB176-1 (DSM 5480) that exhibited higher cry3Aa production and became the active ingredient of the plant protection product Novodor^® FC. Here, we report a comparative genome analysis of the parental strain NB125, its derivative NB176-1 and the current commercial production strain NB176. The entire genome sequences of the parental and derivative strains were deciphered by a hybrid de novo approach using short (Illumina) and long (Nanopore) read sequencing techniques. Genome assembly revealed a chromosome of 5.4 to 5.6 Mbp and six plasmids with a size range from 14.9 to 250.5 kbp for each strain. The major differences among the original NB125 and the derivative strains NB176-1 and NB176 were an additional copy of the cry3Aa gene, which translocated to another plasmid as well as a chromosomal deletion (~ 178 kbp) in NB176. The assembled genome sequences were further analyzed in silico for the presence of virulence and antimicrobial resistance (AMR) genes.

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.

Characterization of the pathogenicity of a Bacillus cereus isolate from the Mariana Trench

Bacillus cereus is an important opportunistic pathogen widely distributed in the environment. In this study, we reported the isolation and characterization of a B. cereus isolate, MB1, from the Challenger Deep of the Mariana Trench. MB1 is aerobic, motile, and able to form endospores. It possesses 5966 genes distributed on a circular chromosome and two plasmids. The MB1 genome contains 14 sets of 23S, 5S, and 16S ribosomal RNA operons, 106 tRNA genes, 4 sRNA genes, 12 genomic islands, 13 prophages, and 302 putative virulence genes, including enterotoxins and cytolysins. Infection studies showed that MB1 was able to cause acute and lethal infection in fish and mice, and was highly toxic to mammalian cells. MB1 induced, in a dose-dependent manner, pyroptotic cell death, characterized by activation of caspase-1, cleavage of gasdermin D, and release of IL-1β and IL-18. MB1 spores exhibited swimming and haemolytic capacity, but were severely attenuated in pathogenicity, which, however, was regained to the full extent when the spores germinated under suitable conditions. Taken together, these results provide new insights into the biological and pathogenic mechanism of deep sea B. cereus.

High-Quality Genome Sequences of Six Actinobacterial Strains Isolated from Granite, Granodiorite, and Tourmaline Rock Surfaces Sampled from Tamil Nadu, India, and New England, United States

Here, we announce four contiguous and two high-quality draft genome sequences of six actinobacterial strains (Blastococcus, Georgenia, Nocardioides, Allobranchiibius, Yimella, and Williamsia) that were isolated from rock samples obtained from Indian historical ruins and colonial building stones in New England, United States. These new sequences expand the genome datasets recovered from stone-dwelling microbes and will allow the prediction of their potential role in the stone microbiome.

Draft Genome Sequences for Bacteria Associated with Root Nodules of Alnus incana in New England

Nine bacterial strains isolated from the root nodules of Alnus incana were sequenced to determine their potential roles in plant health. The selected bacterial isolates belonged to the genera Bacillus, Herbaspirillum , Pantoea , Paenibacillus , and Rothia . Here, we report the draft genome sequences.

Conjugation across Bacillus cereus and kin: A review

Horizontal gene transfer (HGT) is a major driving force in shaping bacterial communities. Key elements responsible for HGT are conjugation-like events and transmissible plasmids. Conjugative plasmids can promote their own transfer as well as that of co-resident plasmids. Bacillus cereus and relatives harbor a plethora of plasmids, including conjugative plasmids, which are at the heart of the group species differentiation and specification. Since the first report of a conjugation-like event between strains of B. cereus sensu lato (s.l.) 40 years ago, many have studied the potential of plasmid transfer across the group, especially for plasmids encoding major toxins. Over the years, more than 20 plasmids from B. cereus isolates have been reported as conjugative. However, with the increasing number of genomic data available, in silico analyses indicate that more plasmids from B. cereus s.l. genomes present self-transfer potential. B. cereus s.l. bacteria occupy diverse environmental niches, which were mimicked in laboratory conditions to study conjugation-related mechanisms. Laboratory mating conditions remain nonetheless simplistic compared to the complex interactions occurring in natural environments. Given the health, economic and ecological importance of strains of B. cereus s.l., it is of prime importance to consider the impact of conjugation within this bacterial group.

Ecophysiological and genomic analyses of a representative isolate of highly abundant Bacillus cereus strains in contaminated subsurface sediments

Bacillus cereus strain CPT56D-587-MTF (CPTF) was isolated from the highly contaminated Oak Ridge Reservation (ORR) subsurface. This site is contaminated with high levels of nitric acid and multiple heavy metals. Amplicon sequencing of the 16S rRNA genes (V4 region) in sediment from this area revealed an amplicon sequence variant (ASV) with 100% identity to the CPTF 16S rRNA sequence. Notably, this CPTF-matching ASV had the highest relative abundance in this community survey, with a median relative abundance of 3.77% and comprised 20%-40% of reads in some samples. Pangenomic analysis revealed that strain CPTF has expanded genomic content compared to other B. cereus species-largely due to plasmid acquisition and expansion of transposable elements. This suggests that these features are important for rapid adaptation to native environmental stressors. We connected genotype to phenotype in the context of the unique geochemistry of the site. These analyses revealed that certain genes (e.g. nitrate reductase, heavy metal efflux pumps) that allow this strain to successfully occupy the geochemically heterogenous microniches of its native site are characteristic of the B. cereus species while others such as acid tolerance are mobile genetic element associated and are generally unique to strain CPTF.

Identification of Universally Applicable and Species-Specific Marker Peptides for Bacillus anthracis

Anthrax is a zoonotic infection caused by the bacterium Bacillus anthracis (BA). Specific identification of this pathogen often relies on targeting genes located on two extrachromosomal plasmids, which represent the major pathogenicity factors of BA. However, more recent findings show that these plasmids have also been found in other closely related Bacillus species. In this study, we investigated the possibility of identifying species-specific and universally applicable marker peptides for BA. For this purpose, we applied a high-resolution mass spectrometry-based approach for 42 BA isolates. Along with the genomic sequencing data and by developing a bioinformatics data evaluation pipeline, which uses a database containing most of the publicly available protein sequences worldwide (UniParc), we were able to identify eleven universal marker peptides unique to BA. These markers are located on the chromosome and therefore, might overcome known problems, such as observable loss of plasmids in environmental species, plasmid loss during cultivation in the lab, and the fact that the virulence plasmids are not necessarily a unique feature of BA. The identified chromosomally encoded markers in this study could extend the small panel of already existing chromosomal targets and along with targets for the virulence plasmids, may pave the way to an even more reliable identification of BA using genomics- as well as proteomics-based techniques.

Draft Genomes Sequences of 11 Geodermatophilaceae Strains Isolated from Building Stones from New England and Indian Stone Ruins found at historic sites in Tamil Nadu, India

Metagenomic analysis of stone microbiome from samples collected in New England, USA and Tamil Nadu, India identified numerous Actinobacteria including Geodermatphilaceae. A culture-dependent approach was performed as a companion study with this culture-independent metagenomic analysis of these stone samples and resulted in the isolation of eleven Geodermatphilaceae strains (2 Geodermatophilus and 9 Blastococcus strains). The genomes of the 11 Geodermatphilaceae strains were sequenced and analyzed. The genomes for the two Geodermatophilus isolates, DF1-2 and TF2-6, were 4.45 and 4.75 Mb, respectively, while the Blastococcus genomes ranged in size from 3.98 to 5.48 Mb. Phylogenetic analysis, digital DNA:DNA hybridization (dDDH), and comparisons of the average nucleotide identities (ANI) suggest the isolates represent novel Geodermatophilus and Blastococcus species. Functional analysis of the Geodermatphilaceae genomes provides insight on the stone microbiome niche.

The SGNH hydrolase family: a template for carbohydrate diversity

The substitution and de-substitution of carbohydrate materials are important steps in the biosynthesis and/or breakdown of a wide variety of biologically important polymers. The SGNH hydrolase superfamily is a group of related and well-studied proteins with a highly conserved catalytic fold and mechanism composed of 16 member families. SGNH hydrolases can be found in vertebrates, plants, fungi, bacteria, and archaea, and play a variety of important biological roles related to biomass conversion, pathogenesis, and cell signaling. The SGNH hydrolase superfamily is chiefly composed of a diverse range of carbohydrate-modifying enzymes, including but not limited to the carbohydrate esterase families 2, 3, 6, 12 and 17 under the carbohydrate-active enzyme classification system and database (CAZy.org). In this review, we summarize the structural and functional features that delineate these subfamilies of SGNH hydrolases, and which generate the wide variety of substrate preferences and enzymatic activities observed of these proteins to date.

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

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

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.

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

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

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.

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.

Genomic insights into the diversity of non-coding RNAs in Bacillus cereus sensu lato

Bacillus cereus sensu lato is a group of bacteria of medical and agricultural importance in different ecological niches and with controversial taxonomic relationships. Studying the composition of non-coding RNAs (ncRNAs) in several bacterial groups has been an important tool for identifying genetic information and better understanding genetic regulation towards environment adaptation. However, to date, no comparative genomics study of ncRNA has been performed in this group. Thus, this study aimed to identify and characterize the set of ncRNAs from 132 strains of Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis to obtain an overview of the diversity and distribution of these genetic elements in these species. We observed that the number of ncRNAs differs in the chromosomes of the three species, but not in the plasmids, when species or phylogenetic clusters were compared. The prevailing functional/structural category was Cis-reg and the most frequent class was Riboswitch. However, in plasmids, the class Group II intron was the most frequent. Also, nine ncRNAs were selected for validation in the strain B. thuringiensis 407 by RT-PCR, which allowed to identify the expression of the ncRNAs. The wide distribution and diversity of ncRNAs in the B. cereus group, and more intensely in B. thuringiensis, may help improve the abilities of these species to adapt to various environmental changes. Further studies should address the expression of these genetic elements in different conditions.

The Fate of Bacteria of the Bacillus cereus Group in the Amoeba Environment

The Bacillus cereus sensu lato group consists of several closely related species, including B. anthracis, B. cereus sensu stricto, and B. thuringiensis. Spores of these pathogenic bacteria are commonly found in the soil but evidence suggests that they are unable to grow in such a natural environment in the absence of nutrient input. Amoebas have been reported to be an amplifier for several species of pathogenic bacteria and their potential involvement to explain the large amount of B. thuringiensis and B. cereus spores in soil has been frequently proposed. Here, we studied the fate of Bacillus and amoebas when cultured together. We show that the virulence factors produced by B. thuringiensis and B. cereus do not affect the amoeba Acanthamoeba castellanii, which, on the contrary, can phagocytose and effectively digest vegetative Bacillus cells to grow and prevent the formation of cysts. Bacterial spores can germinate in the amoeba environment and the vegetative cells can then form chains or aggregates that appear to be less efficiently phagocyted by the amoeba. The use of transcriptional fusions between fluorescent reporter genes and stationary phase- and sporulation-specific promoters showed that the sporulation process occurs more efficiently in the presence of amoebas than in their absence. Moreover, our results showed the amoeba environment to promote spore germination and allow the bacteria to complete their developmental cycle. Overall, this study suggests that the amoeba-Bacillus interaction creates a virtuous circle in which each protagonist helps the other to develop.

Differentiation of Bacillus cereus and Bacillus thuringiensis Using Genome-Guided MALDI-TOF MS Based on Variations in Ribosomal Proteins

Bacillus cereus and B. thuringiensis are closely related species that are relevant to foodborne diseases and biopesticides, respectively. Unambiguous differentiation of these two species is crucial for bacterial taxonomy. As genome analysis offers an objective but time-consuming classification of B. cereus and B. thuringiensis, in the present study, matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) was used to accelerate this process. By combining in silico genome analysis and MALDI-TOF MS measurements, four species-specific peaks of B. cereus and B. thuringiensis were screened and identified. The species-specific peaks of B. cereus were m/z 3211, 6427, 9188, and 9214, and the species-specific peaks of B. thuringiensis were m/z 3218, 6441, 9160, and 9229. All the above peaks represent ribosomal proteins, which are conserved and consistent with the phylogenetic relationship between B. cereus and B. thuringiensis. The specificity of the peaks was robustly verified using common foodborne pathogens. Thus, we concluded that genome-guided MALDI-TOF MS allows high-throughput differentiation of B. cereus and B. thuringiensis and provides a framework for differentiating other closely related species.

Draft Genome Sequences of 6 Actinobacterial Strains Isolated from Rock Surfaces Obtained from Indian Stone Ruins in Tamil Nadu, India, and Rocks from New England, United States

Here, we report the draft genome sequences obtained for 6 actinobacterial strains isolated from stone surfaces acquired from New England and Indian ruins. These strains were sequenced to determine their potential functional roles in the stone microbiome. The strains belong to the genera Allobranchiibius, Agrococcus, Dermococcus, Leifsonia, and Mycobacterium.

Development of a high resolution melting method based on a novel molecular target for discrimination between Bacillus cereus and Bacillus thuringiensis

Delimitation within the Bacillus cereus group is confusing due to the highly similar genetic background of its constituent bacteria. This study aimed to develop a rapid and efficient method for the identification of Bacillus cereus and Bacillus thuringiensis, two closely related species within the B. cereus group. Using average nucleotide identity analysis (ANI) and ribosomal multilocus sequence typing (rMLST), the authenticity of the genomes of B. cereus and B. thuringiensis was determined. Emetic B. cereus and Bacillus bombysepticus were also included to provide novel genomic insights into the boundaries within the B. cereus group. Using pan-genome analysis, ispD, a novel core and single-copy molecular target, was identified for the differentiation between B. cereus and B. thuringiensis. Based on the single nucleotide polymorphism within ispD, a high resolution melting (HRM) method for the determination of B. cereus and B. thuringiensis was developed. This method can not only distinguish B. cereus and B. thuringiensis, but can also separate B. cereus from other foodborne pathogenic bacteria. The detection limit of this method could reach 1 pg of pure genomic DNA and 3.7 × 10² cfu/mL of pure culture. Moreover, this new method could effectively differentiate B. cereus and B. thuringiensis in spiked, mixed, and real food samples. Collectively, the established HRM method can provide a new reference paradigm for the sensitive and specific nucleic acid detection of pathogens with identical genomes.

Biofilm-associated heat resistance of Bacillus cereus spores in vitro and in a food model, Cheonggukjang jjigae

Bacillus cereus, a foodborne pathogen, is capable of forming spores and biofilms as methods to withstand environmental stresses. These bacterial structures are an issue for food safety as they aid the bacteria survive heat sterilisation processes of foods and food contact surfaces. This study was conducted to investigate the role of the biofilm structure in providing an extra layer of protection to spores against heat treatments. For this, heat resistance of B. cereus spores in intact biofilms was compared to that of planktonic spores in vitro and in a Cheonggukjang jjigae food model. Using methods developed in this study to measure the wet and dry heat resistance of spores in intact biofilms, it was found that B. cereus spores have significantly higher heat resistances when present in biofilms rather than as planktonic spores, and that dry heat is less effective than wet heat at killing spores in biofilms. In further detail, for wet heat treatments, spores in biofilms of the strain isolated from Cheonggukjang (Korean fermented whole soybean), B. cereus CH3, had generally higher wet heat resistances than the reference strain, B. cereus ATCC 10987, both in vitro and in the Cheonggukjang jjigae food model. However, the spores in biofilms of the two strains showed similar heat resistance to dry heat, with some exceptions, when biofilms were formed in vitro or in Cheonggukjang jjigae broth. Meanwhile, B. cereus ATCC 10987 spores in biofilms had higher or similar wet heat resistances in vitro compared to in Cheonggukjang jjigae broth. Wet heat resistances of B. cereus CH3 spores in biofilms were all statistically similar regardless of biofilm formation media (brain heart infusion and Cheonggukjang jjigae broths). For dry heat, spores in biofilms of both B. cereus strains were more heat resistant when biofilms were formed in the Cheonggukjang jjigae food model rather than in vitro. Altogether, heat resistances of spores in biofilms formed in vitro and in the food environment were found to be different depending on the tested B. cereus strain, but higher than planktonic spores in any case. This is the first study examining the heat resistance of B. cereus spores in intact biofilms matrices attached to the surface, both in vitro and in a food model. Therefore, this research is valuable to understand the protective effects of biofilms formed in food environments and to reduce the food safety risks associated with B. cereus.

Profiling of Bacillus cereus on Canadian grain

Microorganisms that cause foodborne illnesses challenge the food industry; however, environmental studies of these microorganisms on raw grain, prior to food processing, are uncommon. Bacillus cereus sensu lato is a diverse group of bacteria that is common in our everyday environment and occupy a wide array of niches. While some of these bacteria are beneficial to agriculture due to their entomopathogenic properties, others can cause foodborne illness; therefore, characterization of these bacteria is important from both agricultural and food safety standpoints. We performed a survey of wheat and flax grain samples in 2018 (n = 508) and 2017 (n = 636) and discovered that B. cereus was present in the majority of grain samples, as 56.3% and 85.2%, in two years respectively. Whole genome sequencing and comparative genomics of 109 presumptive B. cereus isolates indicates that most of the isolates were closely related and formed two genetically distinct groups. Comparisons to the available genomes of reference strains suggested that the members of these two groups are not closely related to strains previously reported to cause foodborne illness. From the same data set, another, genetically more diverse group of B. cereus was inferred, which had varying levels of similarity to previously reported strains that caused disease. Genomic analysis and PCR amplification of genes linked to toxin production indicated that most of the isolates carry the genes nheA and hbID, while other toxin genes and gene clusters, such as ces, were infrequent. This report of B. cereus on grain from Canada is the first of its kind and demonstrates the value of surveillance of bacteria naturally associated with raw agricultural commodities such as cereal grain and oilseeds.

Trichococcus shcherbakoviae subsp. psychrophilus subsp. nov., a psychrotolerant facultative anaerobe isolated from a cold spring

A psychrotolerant facultative anaerobe, strain SKBG^T, was isolated from the bottom sediments of the cold mineral spring Buxichen (Buryatia, Russia). Gram-positive non-motile cocci with a diameter of 1.75-2.5 µm were observed singly or in long chains. Cells grew in the temperature range from ̶ 5-35 °C. Growth was observed within the pH range of 7.0-9.5, with the optimum growth at pH 7.6 and at a NaCl concentration from 0-1.0 % (optimum 0.1 % (w/v)). Strain SKBG^T was a chemoorganoheterotroph that used sugars and some organic acids as substrates. The predominant fatty acids in cell walls were С_16:1ω9, С_18:1ω9, and С_16 : 0. The 16S rRNA gene sequence of strain SKBG^T shared high similarity (>99 %) with those of the type strains of the genus Trichococcus. Digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values between strain SKBG^T and Trichococcus shcherbakoviae Art^T (=DSM 107162^T=VKM B-3260^T) were 70.1 and 95.4 %, respectively. The genomic DNA G+C content of strain SKBG^T was 47.1 mol%. Compared with the type strain of T. shcherbakoviae, the new strain was characterized by a temperature optimum for growth (10 °C) significantly lower than that of T. shcherbakoviae DSM 107162^T (20-30 °C). Based on phenotypic and genomic characteristics, the isolate SKBG^T was classified as T. shcherbakoviae subsp. psychrophilus subsp. nov. The type strain is SKBG^T (=VKM B-3241^Т=JCM 33326^T).

Comprehensive Genomic Analysis of the Endophytic Bacillus altitudinis Strain GLB197, a Potential Biocontrol Agent of Grape Downy Mildew

Bacillus has been extensively studied for agricultural application as a biocontrol agent. B. altitudinis GLB197, an endophytic bacterium isolated from grape leaves, exhibits distinctive inhibition to grape downy mildew based on unknown mechanisms. To determine the genetic traits involved in the mechanism of biocontrol and host-interaction traits, the genome sequence of GLB197 was obtained and further analyzed. The genome of B. altitudinis GLB197 consisted of one plasmid and a 3,733,835-bp circular chromosome with 41.56% G + C content, containing 3,770 protein-coding genes. Phylogenetic analysis of 17 Bacillus strains using the concatenated 1,226 single-copy core genes divided into different clusters was conducted. In addition, average nucleotide identity (ANI) values indicate that the current taxonomy of some B. pumilus group strains is incorrect. Comparative analysis of B. altitudinis GLB197 proteins with other B. altitudinis strains identified 3,157 core genes. Furthermore, we found that the pan-genome of B. altitudinis is open. The genome of B. altitudinis GLB197 contains one nonribosomal peptide synthetase (NRPS) gene cluster which was annotated as lichenysin. Interestingly, the cluster in B. altitudinis has two more genes than other Bacillus strains (lgrD and lgrB). The two genes were probably obtained via horizontal gene transfer (HGT) during the evolutionary process from Brevibacillus. Taken together, these observations enable the future application of B. altitudinis GLB197 as a biocontrol agent for control of grape downy mildew and promote our understanding of the beneficial interactions between B. altitudinis GLB197 and plants.

Genotype and virulence gene analyses of Bacillus cereus group clinical isolates from the Chinese softshell turtle (Pelodiscus sinensis) in Taiwan

Chinese softshell turtles (Pelodiscus sinensis) (CST) are susceptible to infections by bacteria belonging to the Bacillus cereus group (Bcg). Bcg includes several closely related species, two of which, B. cereus and B. thuringiensis, are pathogens of aquatic animals or insects. In the present study, we collected 57 Bcg isolates obtained from diseased CST from 2016 to 2019 in Kaohsiung and Pingtung, the areas with the most CST farms in Taiwan. All isolates were divided into four genotypes with two restriction enzymes, SmaI and NotI, by pulsed-field gel electrophoresis and enterobacterial repetitive intergenic consensus-polymerase chain reaction (ERIC-PCR). Representative isolates from each genotype were subjected to phylogenetic tree analysis using 16S rDNA and pyruvate carboxylase genes as phylogenetic markers, and these CST isolates appeared in different clades. PCR was performed targeting six selected virulence genes, four of which were detected in CST isolates, including cytotoxin K (1/57), hblC of the haemolysin BL complex (46/57), nheA of the non-haemolytic enterotoxin complex (52/57) and enterotoxin FM (57/57), whereas cereulide synthetase and cereulide peptide synthase-like genes were not detected in any isolates.

Phenotypic plasticity and evolution of thermal tolerance in bacteria from temperate and hot spring environments

Phenotypic plasticity allows individuals to respond to the selective forces of a new environment, followed by adaptive evolution. We do not know to what extent phenotypic plasticity allows thermal tolerance evolution in bacteria at the border of their physiological limits. We analyzed growth and reaction norms to temperature of strains of two bacterial lineages, Bacillus cereus sensu lato and Bacillus subtilis sensu lato, that evolved in two contrasting environments, a temperate lagoon (T) and a hot spring (H). Our results showed that despite the co-occurrence of members of both lineages in the two contrasting environments, norms of reactions to temperature exhibited a similar pattern only in strains within the lineages, suggesting fixed phenotypic plasticity. Additionally, strains from the H environment showed only two to three degrees centigrade more heat tolerance than strains from the T environment. Their viability decreased at temperatures above their optimal for growth, particularly for the B. cereus lineage. However, sporulation occurred at all temperatures, consistent with the known cell population heterogeneity that allows the Bacillus to anticipate adversity. We suggest that these mesophilic strains survive in the hot-spring as spores and complete their life cycle of germination and growth during intermittent opportunities of moderate temperatures. The limited evolutionary changes towards an increase in heat tolerance in bacteria should alert us of the negative impact of climate change on all biological cycles in the planet, which at its most basic level depends on microorganisms.

Pan-Genome Portrait of Bacillus mycoides Provides Insights into the Species Ecology and Evolution

Bacillus mycoides is poorly known despite its frequent occurrence in a wide variety of environments. To provide direct insight into its ecology and evolutionary history, a comparative investigation of the species pan-genome and the functional gene categorization of 35 isolates obtained from soil samples from northeastern Poland was performed. The pan-genome of these isolates is composed of 20,175 genes and is characterized by a strong predominance of adaptive genes (∼83%), a significant amount of plasmid genes (∼37%), and a great contribution of prophages and insertion sequences. The pan-genome structure and phylodynamic studies had suggested a wide genomic diversity among the isolates, but no correlation between lineages and the bacillus origin was found. Nevertheless, the two B. mycoides populations, one from Białowieża National Park, the last European natural primeval forest with soil classified as organic, and the second from mineral soil samples taken in a farm in Jasienówka, a place with strong anthropogenic pressure, differ significantly in the frequency of genes encoding proteins enabling bacillus adaptation to specific stress conditions and production of a set of compounds, thus facilitating their colonization of various ecological niches. Furthermore, differences in the prevalence of essential stress sigma factors might be an important trail of this process. Due to these numerous adaptive genes, B. mycoides is able to quickly adapt to changing environmental conditions.

IMPORTANCE This research allows deeper understanding of the genetic organization of natural bacterial populations, specifically, Bacillus mycoides, a psychrotrophic member of the Bacillus cereus group that is widely distributed worldwide, especially in areas with continental cold climates. These thorough analyses made it possible to describe, for the first time, the B. mycoides pan-genome, phylogenetic relationship within this species, and the mechanisms behind the species ecology and evolutionary history. Our study indicates a set of functional properties and adaptive genes, in particular, those encoding sigma factors, associated with B. mycoides acclimatization to specific ecological niches and changing environmental conditions.

Bacillus cereus in Packaging Material: Molecular and Phenotypical Diversity Revealed

The Bacillus cereus group has been isolated from soils, water, plants and numerous food products. These species can produce a variety of toxins including several enterotoxins [non-hemolytic enterotoxin (Nhe), hemolysin BL (Hbl), cytotoxin K, and enterotoxin FM], the emetic toxin cereulide and insecticidal Bt toxins. This is the first study evaluating the presence of B. cereus in packaging material. Among 75 different isolates, four phylogenetic groups were detected (II, III, IV, and VI), of which the groups III and IV were the most abundant with 46.7 and 41.3%, respectively. One isolate was affiliated to psychrotolerant group VI. Growth experiments showed a mesophilic predominance. Based on PCR analysis, nhe genes were detectable in 100% of the isolates, while hbl genes were only found in 50.7%. The cereulide encoding gene was found in four out of 75 isolates, no isolate carried a crystal toxin gene. In total, thirteen different toxin gene profiles were identified. We showed that a variety of B. cereus group strains can be found in packaging material. Here, this variety lies in the presence of four phylogenetic groups, thirteen toxin gene profiles, and different growth temperatures. The results suggest that packaging material does not contain significant amounts of highly virulent strains, and the low number of cereulide producing strains is in accordance with other results.