Linking Bacillus cereus Genotypes and Carbohydrate Utilization Capacity

Comparative Genomic and Functional Analyses for Insights into Pantoea agglomerans Strains Adaptability in Diverse Ecological Niches

Pantoea agglomerans inhabit diverse ecological niches, ranging from epiphytes and endophytes in plants, body of animals, and occasionally in the human system. This multifaceted bacterium contributes substantially to plant growth promotion, stress resilience, and biocontrol but can also act as a pathogen to its host. The genetic determinants underlying these diverse functions remain largely unfathomed and to uncover this phenomenon, nineteen strains of Pantoea agglomerans were selected and analyzed. Genome-to-Genome Distance Calculator (GGDC) which uses the Genome Blast Distance Phylogeny (GBDP) technique to calculate digital DDH values. Phylogenetic analysis via Genome-to-Genome distance, Average Nucleotide Identity, and Amino Acid Identity calculation revealed that all strains belonged to the genus Pantoea. However, strain 33.1 had a lower value than the threshold for the same species delineation. Bacterial Pan Genome Analysis (BPGA) Pipeline and MinPath analysis revealed genetic traits associated with environmental resilience, such as oxidative stress, UV radiation, temperature extremes, and metabolism of distinct host-specific carbohydrates. Protein-protein interactome analysis illustrated osmotic stress proteins closely linked with core proteins, while heavy metal tolerance, nitrogen metabolism, and Type III and VI secretion systems proteins generally associated with pathogenicity formed a separate network, indicating strain-specific characteristics. These findings shed new light on the intricate genetic architecture of Pantoea agglomerans, revealing its adaptability to inhabit diverse niches and thrive in varied environments.

The potential of lactose to inhibit cereulide biosynthesis of emetic Bacillus cereus in milk

This study aims to investigate the potential role of lactose on cereulide biosynthesis by emetic Bacillus cereus in dairy matrices. The cereulide yields in whole milk and lactose-free milk were investigated using the emetic reference strain F4810/72. To eliminate the influence of complex food substrates, the LB medium model was further used to characterize the effect of lactose on cereulide produced by F4810/72 and five other emetic B. cereus strains. Results showed that the lactose-free milk displayed a 13-fold higher amount of cereulide than whole milk, but the cereulide level could be reduced by 91 % when the lactose content was restored. The significant inhibition of lactose on cereulide yields of all tested B. cereus strains was observed in LB medium, showing a dose-dependent manner with inhibition rates ranging of 89-98 %. The growth curves and lactose utilization patterns of all strains demonstrated that B. cereus cannot utilize lactose as a carbon source and lactose might act as a signal molecule to regulate cereulide production. Moreover, lactose strongly repressed the expression of cereulide synthetase genes (ces), possibly by inhibiting the key regulator Spo0A at the transcriptional level. Our findings highlight the potential of lactose as an effective strategy to control cereulide production in food.

Bifidobacterium longum subsp. infantis regulates Th1/Th2 balance through the JAK-STAT pathway in growing mice

Objectives: Bifidobacterium longum subsp. infantis is a dominant bacterium in infant gut, which plays a critical role in maintaining the health and development of infants. This study investigated the abilities of eight different strains of B. longum subsp. infantis to regulate the T helper (Th)1/Th2 balance. Methods: Eight B. longum subsp. infantis strains, including I2MI (FJSWXI2MIM1), I4MI [FJSWXI4MI (CCFM1270)], I4MNI (FJSWXI4MNIM1), I5TI (FJSWXI5TIM1), I6TI (FJSWXI6TIM1), I8TI [FJSWXI8TI (CCFM1271)], I10TI [FJSWXI10TI (CCFM1272)], and B6MNI [BJSWXB6MNIM1 (CCFM1269)], were gavaged to BALB/C pups in both female (n = 8) and male (n = 8) mice starting from 1 to 3 weeks old (1 × 109 CFU/day/mice). Selected immune cells were assessed by immunofluorescence and flow cytometry. Cytokines and immunoglobulins were determined by ELISA. Bacterial and bifidobacterial communities were determined by 16S rRNA gene sequencing and bifidobacterial groEL sequencing. Results: B. longum subsp. infantis I4MI and I8TI were shown to increase the ration of colonic IgG2a/IgE in male mice (P < 0.05). B6MNI was demonstrated to significantly increase the levels of colonic IFN-γ and IgG2a, as well as the ratio of IgG2a/IgE in female mice (P < 0.05). It was also shown to significantly increase the ratio of colonic IgG2a/IgE (P < 0.05) and reduce the level of colonic IL-4 in male mice (P < 0.05). Furthermore, B6MNI was demonstrated to regulate colonic JAK/STAT pathway in both male and female mice. I4MI, I5TI, and B6MNI were shown to increase the relative abundance of Bifidobacterium and B. longum subsp. infantis in both male and female mice, whereas I8TI was only shown to increase the relative abundance of Bifidobacterium and B. longum subsp. infantis in male mice (P < 0.05). Conclusion: These results indicated supplementation with B. longum subsp. infantis in early infancy may regulate the Th1/Th2 immune balance, which may prevent the development of related diseases.

Targeting the Impossible: A Review of New Strategies against Endospores

Endospore-forming bacteria are ubiquitous, and their endospores can be present in food, in domestic animals, and on contaminated surfaces. Many spore-forming bacteria have been used in biotechnological applications, while others are human pathogens responsible for a wide range of critical clinical infections. Due to their resistant properties, it is challenging to eliminate spores and avoid the reactivation of latent spores that may lead to active infections. Furthermore, endospores play an essential role in the survival, transmission, and pathogenesis of some harmful strains that put human and animal health at risk. Thus, different methods have been applied for their eradication. Nevertheless, natural products are still a significant source for discovering and developing new antibiotics. Moreover, targeting the spore for clinical pathogens such as Clostridioides difficile is essential to disease prevention and therapeutics. These strategies could directly aim at the structural components of the spore or their germination process. This work summarizes the current advances in upcoming strategies and the development of natural products against endospores. This review also intends to highlight future perspectives in research and applications.

MsmR1, a global transcription factor, regulates polymyxin synthesis and carbohydrate metabolism in Paenibacillus polymyxa SC2

The multiple-sugar metabolism regulator (MsmR), a transcription factor belonging to the AraC/XylS family, participates in polysaccharide metabolism and virulence. However, the transcriptional regulatory mechanisms of MsmR1 in Paenibacillus polymyxa remain unclear. In this study, knocking out msmR1 was found to reduce polymyxin synthesis by the SC2-M1 strain. Chromatin immunoprecipitation assay with sequencing (ChIP-seq) revealed that most enriched pathway was that of carbohydrate metabolism. Additionally, electromobility shift assays (EMSA) confirmed the direct interaction between MsmR1 and the promoter regions of oppC3, sucA, sdr3, pepF, yycN, PPSC2_23180, pppL, and ydfp. MsmR1 stimulates polymyxin biosynthesis by directly binding to the promoter regions of oppC3 and sdr3, while also directly regulating sucA and influencing the citrate cycle (TCA cycle). In addition, MsmR1 directly activates pepF and was beneficial for spore and biofilm formation. These results indicated that MsmR1 could regulate carbohydrate and amino acid metabolism, and indirectly affect biological processes such as polymyxin synthesis, biofilm formation, and motility. Moreover, MsmR1 could be autoregulated. Hence, this study expand the current knowledge of MsmR1 and will be beneficial for the application of P. polymyxa SC2 in the biological control against the certain pathogens in pepper.

Characterization and Genetic Diversity of Bacillus cereus Strains Isolated from Baby Wipes

Bacillus cereus, a ubiquitous environmental microorganism known to cause foodborne illness, was isolated from samples taken from imported baby wipes from two different countries. These strains were characterized using a comprehensive molecular approach involving endpoint PCR, whole genome sequencing (WGS), comparative genomics, and biochemical analyses. A multiplex endpoint PCR assay was used to identify the enterotoxins: hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Phylogenetically, the strains clustered into two major groups according to sequence type (ST) and singleton. We used the Center for Food Safety and Applied Nutrition (CFSAN) GalaxyTrakr BTyper computational tool to characterize the strains further. As an additional means of characterization, we investigated the possible role of carbohydrate transport systems and their role in nutrient uptake by performing a BLAST analysis of the 40 B. cereus genomes recovered from baby wipes. This study outlines a multifaceted workflow that uses the analysis of enterotoxigenic potential, bioinformatics, genomic diversity, genotype, phenotype, and carbohydrate utilization as a comprehensive strategy to characterize these B. cereus strains isolated from baby wipes and further our understanding of the phylogenetic relatedness of strains associated with baby wipe production facilities that could potentially pose an infection risk to a vulnerable infant population.

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

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

Metabolic profiling reveals nutrient preferences during carbon utilization in Bacillus species

The genus Bacillus includes species with diverse natural histories, including free-living nonpathogenic heterotrophs such as B. subtilis and host-dependent pathogens such as B. anthracis (the etiological agent of the disease anthrax) and B. cereus, a cause of food poisoning. Although highly similar genotypically, the ecological niches of these three species are mutually exclusive, which raises the untested hypothesis that their metabolism has speciated along a nutritional tract. Here, we developed a pipeline for quantitative total assessment of the use of diverse sources of carbon for general metabolism to better appreciate the "culinary preferences" of three distinct Bacillus species, as well as related Staphylococcus aureus. We show that each species has widely varying metabolic ability to utilize diverse sources of carbon that correlated to their ecological niches. This approach was applied to the growth and survival of B. anthracis in a blood-like environment and find metabolism shifts from sugar to amino acids as the preferred source of energy. Finally, various nutrients in broth and host-like environments are identified that may promote or interfere with bacterial metabolism during infection.

Adaptation of Bacillus thuringiensis to Plant Colonization Affects Differentiation and Toxicity

The Bacillus cereus group (Bacillus cereus sensu lato) has a diverse ecology, including various species that are vertebrate or invertebrate pathogens. Few isolates from the B. cereus group have however been demonstrated to benefit plant growth. Therefore, it is crucial to explore how bacterial development and pathogenesis evolve during plant colonization. Herein, we investigated Bacillus thuringiensis (Cry-) adaptation to the colonization of Arabidopsis thaliana roots and monitored changes in cellular differentiation in experimentally evolved isolates. Isolates from two populations displayed improved iterative ecesis on roots and increased virulence against insect larvae. Molecular dissection and recreation of a causative mutation revealed the importance of a nonsense mutation in the rho transcription terminator gene. Transcriptome analysis revealed how Rho impacts various B. thuringiensis genes involved in carbohydrate metabolism and virulence. Our work suggests that evolved multicellular aggregates have a fitness advantage over single cells when colonizing plants, creating a trade-off between swimming and multicellularity in evolved lineages, in addition to unrelated alterations in pathogenicity. IMPORTANCE Biologicals-based plant protection relies on the use of safe microbial strains. During application of biologicals to the rhizosphere, microbes adapt to the niche, including genetic mutations shaping the physiology of the cells. Here, the experimental evolution of Bacillus thuringiensis lacking the insecticide crystal toxins was examined on the plant root to reveal how adaptation shapes the differentiation of this bacterium. Interestingly, evolution of certain lineages led to increased hemolysis and insect larva pathogenesis in B. thuringiensis driven by transcriptional rewiring. Further, our detailed study reveals how inactivation of the transcription termination protein Rho promotes aggregation on the plant root in addition to altered differentiation and pathogenesis in B. thuringiensis.

A native conjugative plasmid confers potential selective advantages to plant growth-promoting Bacillus velezensis strain GH1-13

The conjugative plasmid (pBV71) possibly confers a selective advantage to Bacillus velezensis strain GH1-13, although a selective marker gene is yet to be identified. Here we show that few non-mucoid wild-type GH1-13 cells are spontaneously converted to mucoid variants with or without the loss of pBV71. Mucoid phenotypes, which contain or lack the plasmid, become sensitive to bacitracin, gramicidin, selenite, and tellurite. Using the differences in antibiotic resistance and phenotype, we isolated a reverse complement (COM) and a transconjugant of strain FZB42 with the native pBV71. Transformed COM and FZB42p cells were similar to the wild-type strain GH1-13 with high antibiotic resistance and slow growth rates on lactose compared to those of mucoid phenotypes. RT-PCR analysis revealed that the expression of plasmid-encoded orphan aspartate phosphatase (pRapD) was coordinated with a new quorum-sensing (QS) cassette of RapF2-PhrF2 present in the chromosome of strain GH1-13, but not in strain FZB42. Multi-omics analysis on wild-type and plasmid-cured cells of strain GH1-13 suggested that the conjugative plasmid expression has a crucial role in induction of early envelope stress response that promotes cell morphogenesis, biofilm formation, catabolite repression, and biosynthesis of extracellular-matrix components and antibiotics for protection of host cell during exponential phase.

Characterization of Bacillus cereus Group Isolates From Human Bacteremia by Whole-Genome Sequencing

Members of the Bacillus cereus group are spore-forming organisms commonly associated with food poisoning and intestinal infections. Moreover, some strains of the group (i.e., B. cereus sensu stricto and Bacillus thuringiensis) can cause bacteremia in humans, mainly in immunocompromised individuals. Here we performed the genetic characterization of 17 human clinical strains belonging to B. cereus group isolated from blood culture. The whole-genome sequencing (WGS) revealed that the isolates were closely related to B. cereus sensu stricto and B. thuringiensis-type strain. Multilocus sequence typing analysis performed on the draft genome revealed the genetic diversity of our isolates, which were assigned to different sequence types. Based on panC nucleotide sequence, the isolates were grouped in the phylogenetic groups III and IV. The NHE, cer, and inhA gene cluster, entA, entFM, plcA, and plcB, were the most commonly detected virulence genes. Although we did not assess the ability to generate biofilm by phenotypic tests, we verified the prevalence of biofilm associated genes using an in silico approach. A high prevalence of pur gene cluster, xerC, clpY, codY, tasA, sipW, sinI, and sigB genes, was found. Genes related to the resistance to penicillin, trimethoprim, and ceftriaxone were identified in most of the isolates. Intriguingly, the majority of these virulence and AMR genes appeared to be evenly distributed among B. cereus s.s. isolates, as well as closely related to B. thuringiensis isolates. We showed the WGS represents a good approach to rapidly characterize B. cereus group strains, being able to give useful information about genetic epidemiology, the presence of virulence and antimicrobial genes, and finally about the potential hazard related to this underestimated risk.

X-ray Structure, Bioinformatics Analysis, and Substrate Specificity of a 6-Phospho-β-glucosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis

In bacteria, mono- and disaccharides are phosphorylated during the uptake processes through the vastly spread transport system phosphoenolpyruvate-dependent phosphotransferase. As an initial step in the phosphorylated disaccharide metabolism pathway, 6-phospho-β-glucosidases and 6-phospho-β-galactosidases play a crucial role by releasing phosphorylated and nonphosphorylated monosaccharides. However, structural determinants for the specificity of these enzymes still need to be clarified. Here, an X-ray structure of a glycoside hydrolase family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglH, was determined at 2.2 Å resolution, and its substrate specificity was investigated. The sequence of BlBglH was compared to the sequences of 58 other GH1 enzymes using sequence alignments, sequence identity calculations, phylogenetic analysis, and motif discovery. Through these various analyses, BlBglH was found to have sequence features characteristic of the 6-phospho-β-glucosidase activity enzymes. Motif and structural observations highlighted the importance of loop L8 in 6-phospho-β-glucosidase activity enzymes. To further affirm enzyme specificity, molecular docking and molecular dynamics simulations were performed using the crystallographic structure of BlBglH. Docking was carried out with a 6-phospho-β-glucosidase enzyme activity positive and negative control ligand, followed by 400 ns of MD simulations. The positive and negative control ligands were PNP6Pglc and PNP6Pgal, respectively. PNP6Pglc maintained favorable interactions within the active site until the end of the MD simulation, while PNP6Pgal exhibited instability. The favorable binding of substrate stabilized the loops that surround the active site. Binding free energy calculations showed that the PNP6Pglc complex had a substantially lower binding energy compared to the PNP6Pgal complex. Altogether, the findings of this study suggest that BlBglH possesses 6-phospho-β-glucosidase enzymatic activity and revealed sequence and structural differences between bacterial GH1 enzymes of various activities.

High-level production of poly-γ-glutamic acid from untreated molasses by Bacillus siamensis IR10

Background

Poly-γ-glutamic acid (γ-PGA) is a promising biopolymer and has been applied in many fields. Bacillus siamensis SB1001 was a newly isolated poly-γ-glutamic acid producer with sucrose as its optimal carbon source. To improve the utilization of carbon source, and then molasses can be effectively used for γ-PGA production, ⁶⁰cobalt gamma rays was used to mutate the genes of B. siamensis SB1001.

Results

Bacillus siamensis IR10 was screened for the production of γ-PGA from untreated molasses. In batch fermentation, 17.86 ± 0.97 g/L γ-PGA was obtained after 15 h, which is 52.51% higher than that of its parent strain. Fed-batch fermentation was performed to further improve the yield of γ-PGA with untreated molasses, yielding 41.40 ± 2.01 g/L of γ-PGA with a productivity of 1.73 ± 0.08 g/L/h. An average γ-PGA productivity of 1.85 g/L/h was achieved in the repeated fed-batch fermentation. This is the first report of such a high γ-PGA productivity. The analysis of the enzyme activities showed that they were affected by the carbon sources, enhanced ICDH and GDH, and decreased ODHC, which are important for γ-PGA production.

Conclusion

These results suggest that untreated molasses can be used for economical and industrial-scale production of γ-PGA by B. siamensis IR10.

Genes Associated With Psychrotolerant Bacillus cereus Group Isolates

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

Porcine Gastric Mucin Triggers Toxin Production of Enteropathogenic Bacillus cereus

Enteropathogenic Bacillus cereus causes foodborne infections due to the production of pore-forming enterotoxins in the intestine. Before that, spores have to be ingested, survive the stomach passage, and germinate. Thus, before reaching epithelial cells, B. cereus comes in contact with the intestinal mucus layer. In the present study, different aspects of this interaction were analyzed. Total RNA sequencing revealed major transcriptional changes of B. cereus strain F837/76 upon incubation with porcine gastric mucin (PGM), comprising genes encoding enterotoxins and further putative virulence factors, as well as proteins involved in adhesion to and degradation of mucin. Indeed, PGM was partially degraded by B. cereus via secreted, EDTA-sensitive proteases. The amount of enterotoxins detectable in culture media supplemented with PGM was also clearly increased. Tests of further strains revealed that enhancement of enterotoxin production upon contact with PGM is broadly distributed among B. cereus strains. Interestingly, evidence was found that PGM can also strain-specifically trigger germination of B. cereus spores and that vegetative cells actively move toward mucin. Overall, our data suggest that B. cereus is well adapted to the host environment due to massive transcriptome changes upon contact with PGM, attributing mucin an important and, thus far, neglected role in pathogenesis.

Salvage of the 5-deoxyribose byproduct of radical SAM enzymes

5-Deoxyribose is formed from 5′-deoxyadenosine, a toxic byproduct of radical S-adenosylmethionine (SAM) enzymes. The degradative fate of 5-deoxyribose is unknown. Here, we define a salvage pathway for 5-deoxyribose in bacteria, consisting of phosphorylation, isomerization, and aldol cleavage steps. Analysis of bacterial genomes uncovers widespread, unassigned three-gene clusters specifying a putative kinase, isomerase, and sugar phosphate aldolase. We show that the enzymes encoded by the Bacillus thuringiensis cluster, acting together in vitro, convert 5-deoxyribose successively to 5-deoxyribose 1-phosphate, 5-deoxyribulose 1-phosphate, and dihydroxyacetone phosphate plus acetaldehyde. Deleting the isomerase decreases the 5-deoxyribulose 1-phosphate pool size, and deleting either the isomerase or the aldolase increases susceptibility to 5-deoxyribose. The substrate preference of the aldolase is unique among family members, and the X-ray structure reveals an unusual manganese-dependent enzyme. This work defines a salvage pathway for 5-deoxyribose, a near-universal metabolite.

5-Deoxyribose is formed from 5′-deoxyadenosine, a toxic byproduct of radical S-adenosylmethionine enzymes. Here, the authors identify and biochemically characterize a bacterial salvage pathway for 5-deoxyribose, consisting of three enzymes, and solve the crystal structure of the key aldolase.

Rapid, High-Throughput Identification of Anthrax-Causing and Emetic Bacillus cereus Group Genome Assemblies via BTyper, a Computational Tool for Virulence-Based Classification of Bacillus cereus Group Isolates by Using Nucleotide Sequencing Data

The Bacillus cereus group comprises nine species, several of which are pathogenic. Differentiating between isolates that may cause disease and those that do not is a matter of public health and economic importance, but it can be particularly challenging due to the high genomic similarity within the group. To this end, we have developed BTyper, a computational tool that employs a combination of (i) virulence gene-based typing, (ii) multilocus sequence typing (MLST), (iii) panC clade typing, and (iv) rpoB allelic typing to rapidly classify B. cereus group isolates using nucleotide sequencing data. BTyper was applied to a set of 662 B. cereus group genome assemblies to (i) identify anthrax-associated genes in non-B. anthracis members of the B. cereus group, and (ii) identify assemblies from B. cereus group strains with emetic potential. With BTyper, the anthrax toxin genes cya, lef, and pagA were detected in 8 genomes classified by the NCBI as B. cereus that clustered into two distinct groups using k-medoids clustering, while either the B. anthracis poly-γ-d-glutamate capsule biosynthesis genes capABCDE or the hyaluronic acid capsule hasA gene was detected in an additional 16 assemblies classified as either B. cereus or Bacillus thuringiensis isolated from clinical, environmental, and food sources. The emetic toxin genes cesABCD were detected in 24 assemblies belonging to panC clades III and VI that had been isolated from food, clinical, and environmental settings. The command line version of BTyper is available at https://github.com/lmc297/BTyper. In addition, BMiner, a companion application for analyzing multiple BTyper output files in aggregate, can be found at https://github.com/lmc297/BMiner.

IMPORTANCEBacillus cereus is a foodborne pathogen that is estimated to cause tens of thousands of illnesses each year in the United States alone. Even with molecular methods, it can be difficult to distinguish nonpathogenic B. cereus group isolates from their pathogenic counterparts, including the human pathogen Bacillus anthracis, which is responsible for anthrax, as well as the insect pathogen B. thuringiensis. By using the variety of typing schemes employed by BTyper, users can rapidly classify, characterize, and assess the virulence potential of any isolate using its nucleotide sequencing data.

Genomic Reconstruction of Carbohydrate Utilization Capacities in Microbial-Mat Derived Consortia

Two nearly identical unicyanobacterial consortia (UCC) were previously isolated from benthic microbial mats that occur in a heliothermal saline lake in northern Washington State. Carbohydrates are a primary source of carbon and energy for most heterotrophic bacteria. Since CO₂ is the only carbon source provided, the cyanobacterium must provide a source of carbon to the heterotrophs. Available genomic sequences for all members of the UCC provide opportunity to investigate the metabolic routes of carbon transfer between autotroph and heterotrophs. Here, we applied a subsystem-based comparative genomics approach to reconstruct carbohydrate utilization pathways and identify glycohydrolytic enzymes, carbohydrate transporters and pathway-specific transcriptional regulators in 17 heterotrophic members of the UCC. The reconstructed metabolic pathways include 800 genes, near a one-fourth of which encode enzymes, transporters and regulators with newly assigned metabolic functions resulting in discovery of novel functional variants of carbohydrate utilization pathways. The in silico analysis revealed the utilization capabilities for 40 carbohydrates and their derivatives. Two Halomonas species demonstrated the largest number of sugar catabolic pathways. Trehalose, sucrose, maltose, glucose, and beta-glucosides are the most commonly utilized saccharides in this community. Reconstructed regulons for global regulators HexR and CceR include central carbohydrate metabolism genes in the members of Gammaproteobacteria and Alphaproteobacteria, respectively. Genomics analyses were supplemented by experimental characterization of metabolic phenotypes in four isolates derived from the consortia. Measurements of isolate growth on the defined medium supplied with individual carbohydrates confirmed most of the predicted catabolic phenotypes. Not all consortia members use carbohydrates and only a few use complex polysaccharides suggesting a hierarchical carbon flow from cyanobacteria to each heterotroph. In summary, the genomics-based identification of carbohydrate utilization capabilities provides a basis for future experimental studies of carbon flow in UCC.

Analysis of Germination Capacity and Germinant Receptor (Sub)clusters of Genome-Sequenced Bacillus cereus Environmental Isolates and Model Strains

Spore germination of 17 Bacillus cereus food isolates and reference strains was evaluated using flow cytometry analysis in combination with fluorescent staining at a single-spore level. This approach allowed for rapid collection of germination data under more than 20 conditions, including heat activation of spores, germination in complex media (brain heart infusion [BHI] and tryptone soy broth [TSB]), and exposure to saturating concentrations of single amino acids and the combination of alanine and inosine. Whole-genome sequence comparison revealed a total of 11 clusters of operons encoding germinant receptors (GRs): GerK, GerI, and GerL were present in all strains, whereas GerR, GerS, GerG, GerQ, GerX, GerF, GerW, and GerZ (sub)clusters showed a more diverse presence/absence in different strains. The spores of tested strains displayed high diversity with regard to their sensitivity and responsiveness to selected germinants and heat activation. The two laboratory strains, B. cereus ATCC 14579 and ATCC 10987, and 11 food isolates showed a good germination response under a range of conditions, whereas four other strains (B. cereus B4085, B4086, B4116, and B4153) belonging to phylogenetic group IIIA showed a very weak germination response even in BHI and TSB media. Germination responses could not be linked to specific (combinations of) GRs, but it was noted that the four group IIIA strains contained pseudogenes or variants of subunit C in their gerL cluster. Additionally, two of those strains (B4086 and B4153) carried pseudogenes in the gerK and gerRI (sub)clusters that possibly affected the functionality of these GRs.

IMPORTANCE

Germination of bacterial spores is a critical step before vegetative growth can resume. Food products may contain nutrient germinants that trigger germination and outgrowth of Bacillus species spores, possibly leading to food spoilage or foodborne illness. Prediction of spore germination behavior is, however, very challenging, especially for spores of natural isolates that tend to show more diverse germination responses than laboratory strains. The approach used has provided information on the genetic diversity in GRs and corresponding subclusters encoded by B. cereus strains, as well as their germination behavior and possible associations with GRs, and it provides a basis for further extension of knowledge on the role of GRs in B. cereus (group member) ecology and transmission to the host.

Adaptation in Bacillus cereus: From Stress to Disease

Bacillus cereus is a food-borne pathogen that causes diarrheal disease in humans. After ingestion, B. cereus experiences in the human gastro-intestinal tract abiotic physical variables encountered in food, such as acidic pH in the stomach and changing oxygen conditions in the human intestine. B. cereus responds to environmental changing conditions (stress) by reversibly adjusting its physiology to maximize resource utilization while maintaining structural and genetic integrity by repairing and minimizing damage to cellular infrastructure. As reviewed in this article, B. cereus adapts to acidic pH and changing oxygen conditions through diverse regulatory mechanisms and then exploits its metabolic flexibility to grow and produce enterotoxins. We then focus on the intricate link between metabolism, redox homeostasis, and enterotoxins, which are recognized as important contributors of food-borne disease.

Comparative Genomics of Iron-Transporting Systems in Bacillus cereus Strains and Impact of Iron Sources on Growth and Biofilm Formation

Iron is an important element for bacterial viability, however it is not readily available in most environments. We studied the ability of 20 undomesticated food isolates of Bacillus cereus and two reference strains for capacity to use different (complex) iron sources for growth and biofilm formation. Studies were performed in media containing the iron scavenger 2,2-Bipyridine. Transcriptome analysis using B. cereus ATCC 10987 indeed showed upregulation of predicted iron transporters in the presence of 2,2-Bipyridine, confirming that iron was depleted upon its addition. Next, the impact of iron sources on growth performance of the 22 strains was assessed and correlations between growth stimulation and presence of putative iron transporter systems in the genome sequences were analyzed. All 22 strains effectively used Fe citrate and FeCl₃ for growth, and possessed genes for biosynthesis of the siderophore bacillibactin, whereas seven strains lacked genes for synthesis of petrobactin. Hemoglobin could be used by all strains with the exception of one strain that lacked functional petrobactin and IlsA systems. Hemin could be used by the majority of the tested strains (19 of 22). Notably, transferrin, ferritin, and lactoferrin were not commonly used by B. cereus for growth, as these iron sources could be used by 6, 3, and 2 strains, respectively. Furthermore, biofilm formation was found to be affected by the type of iron source used, including stimulation of biofilms at liquid-air interphase (FeCl₃ and Fe citrate) and formation of submerged type biofilms (hemin and lactoferrin). Our results show strain variability in the genome-encoded repertoire of iron-transporting systems and differences in efficacy to use complex iron sources for growth and biofilm formation. These features may affect B. cereus survival and persistence in specific niches.