Current knowledge on the microbiota of edible insects intended for human consumption: A state-of-the-art review

Because of their positive nutritional characteristics and low environmental impact, edible insects might be considered a 'food of the future'. However, there are safety concerns associated with the consumption of insects, such as contaminating chemical and biological agents. The possible presence of pathogenic and toxigenic microorganisms is one of the main biological hazards associated with edible insects. This review presents an overview of the microbiota of edible insects, highlighting the potential risks for human health. Detailed information on the microbiota of edible insects from literature published in 2000-2019 is presented. These data show complex ecosystems, with marked variations in microbial load and diversity, among edible insects as well as stable and species-specific microbiota for some of the most popular edible insect species, such as mealworm larvae (Tenebrio molitor) and grasshoppers (Locusta migratoria). Raw edible insects generally contain high numbers of mesophilic aerobes, bacterial endospores or spore-forming bacteria, Enterobacteriaceae, lactic acid bacteria, psychrotrophic aerobes, and fungi, and potentially harmful species (i.e. pathogenic, mycotoxigenic, and spoilage microbes) may be present. Several studies have focused on reducing the microbial contamination of edible insects by applying treatments such as starvation, rinsing, thermal treatments, chilling, drying, fermentation, and marination, both alone and, sometimes, in combination. Although these studies show that various heat treatments were the most efficient methods for reducing microbial numbers, they also highlight the need for species-specific mitigation strategies. The feasibility of using edible insects as ingredients in the food industry in the development of innovative insect-based products has been explored; although, in some cases, the presence of spore-forming bacteria and other food-borne pathogens is a concern. Recent studies have shown that a risk assessment of edible insects should also include an evaluation of the incidence of antibiotic-resistance (AR) genes and antibiotic-resistant microorganisms in the production chain. Finally, as proposed in the literature, microbial hazards should be limited through the implementation of good hygienic practices during rearing, handling, processing, and storage, as well as the implementation of an appropriate HACCP system for edible insect supply chains. Another issue frequently reported in the literature is the need for a legislative framework for edible insect production, commercialisation, and trading, as well as the need for microbiological criteria specifically tailored for edible insects. Microbiological criteria like those already been established for the food safety and hygiene (e.g. those in the European Union food law) of different food categories (e.g. ready-to-eat products) could be applied to edible insect-based products.

Genome Analysis of Entomopathogenic Bacillus sp. ABP14 Isolated from a Lignocellulosic Compost

We report the complete genome sequence of Bacillus sp. strain ABP14 isolated from lignocellulosic compost and selected by its ability in hydrolyzing carboxymethyl cellulose. This strain does not produce a Cry-like protein but showed an insecticidal activity against larvae of Anticarsia gemmatalis (Lepidoptera). Genome-based taxonomic analysis revealed that the ABP14 chromosome is genetically close to Bacillus thuringiensis serovar finitimus YBT020. ABP14 also carries one plasmid which showed no similarity with those from YBT020. Genome analysis of ABP14 identified unique genes related to cell surface structures, cell wall, metabolic competence, and virulence factors that may contribute for its survival and environmental adaptation, as well as its entomopathogenic activity.

Diversity of Bacillus cereus sensu lato mobilome

BACKGROUND

Bacillus cereus sensu lato s.l.) is a group of bacteria displaying close phylogenetic relationships but a high ecological diversity. The three most studied species are Bacillus anthracis, Bacillus cereus sensu stricto and Bacillus thuringiensis. While some species are pathogenic to mammals or associated with food poisoning, Bacillus thuringiensis is a well-known entomopathogenic bacterium used as biopesticide worldwide. B. cereus s.l. also contains a large variety of mobile genetic elements (MGEs).

RESULTS

In this study, we detail the occurrence and plasmid vs. chromosome distribution of several MGEs in 102 complete and annotated genomes of B. cereus s.l. These MGEs include 16 Insertion Sequence (IS) families, the Tn3 family, 18 different Bacillus cereus repeats (BCRs) and 30 known group II introns.

CONCLUSIONS

Our analysis not only shows the diversity of these MGEs among strains of the same species and between different species within the B. cereus s.l. group, but also highlights the potential impact of these elements on the plasticity of the plasmid pool, and the TEs (Transposable Elements) - species relationship within B. cereus s.l.

Multifaceted toxin profile, an approach toward a better understanding of probiotic Bacillus cereus

Strains of the Bacillus cereus group have been widely used as probiotics for human beings, food animals, plants, and environmental remediation. Paradoxically, B. cereus is responsible for both gastrointestinal and nongastrointestinal syndromes and represents an important opportunistic food-borne pathogen. Toxicity assessment is a fundamental issue to evaluate safety of probiotics. Here, we summarize the state of our current knowledge about the toxins of B. cereus sensu lato to be considered for safety assessment of probiotic candidates. Surfactin-like emetic toxin (cereulide) and various enterotoxins including nonhemolytic enterotoxin, hemolysin BL, and cytotoxin K are responsible for food poisoning outbreaks characterized by emesis and diarrhea. In addition, other factors, such as hemolysin II, Certhrax, immune inhibitor A1, and sphingomyelinase, contribute to toxicity and overall virulence of B. cereus.

The Bacillus cereus Group: Bacillus Species with Pathogenic Potential

The Bacillus cereus group includes several Bacillus species with closely related phylogeny. The most well-studied members of the group, B. anthracis, B. cereus, and B. thuringiensis, are known for their pathogenic potential. Here, we present the historical rationale for speciation and discuss shared and unique features of these bacteria. Aspects of cell morphology and physiology, and genome sequence similarity and gene synteny support close evolutionary relationships for these three species. For many strains, distinct differences in virulence factor synthesis provide facile means for species assignment. B. anthracis is the causative agent of anthrax. Some B. cereus strains are commonly recognized as food poisoning agents, but strains can also cause localized wound and eye infections as well as systemic disease. Certain B. thuringiensis strains are entomopathogens and have been commercialized for use as biopesticides, while some strains have been reported to cause infection in immunocompromised individuals. In this article we compare and contrast B. anthracis, B. cereus, and B. thuringiensis, including ecology, cell structure and development, virulence attributes, gene regulation and genetic exchange systems, and experimental models of disease.

Detection and Identification of Bacillus cereus, Bacillus cytotoxicus, Bacillus thuringiensis, Bacillus mycoides and Bacillus weihenstephanensis via Machine Learning Based FTIR Spectroscopy

The Bacillus cereus group comprises genetical closely related species with variable toxigenic characteristics. However, detection and differentiation of the B. cereus group species in routine diagnostics can be difficult, expensive and laborious since current species designation is linked to specific phenotypic characteristic or the presence of species-specific genes. Especially the differentiation of Bacillus cereus and Bacillus thuringiensis, the identification of psychrotolerant Bacillus mycoides and Bacillus weihenstephanensis, as well as the identification of emetic B. cereus and Bacillus cytotoxicus, which are both producing highly potent toxins, is of high importance in food microbiology. Thus, we investigated the use of a machine learning approach, based on artificial neural network (ANN) assisted Fourier transform infrared (FTIR) spectroscopy, for discrimination of B. cereus group members. The deep learning tool box of Matlab was employed to construct a one-level ANN, allowing the discrimination of the aforementioned B. cereus group members. This model resulted in 100% correct identification for the training set and 99.5% correct identification overall. The established ANN was applied to investigate the composition of B. cereus group members in soil, as a natural habitat of B. cereus, and in food samples originating from foodborne outbreaks. These analyses revealed a high complexity of B. cereus group populations, not only in soil samples but also in the samples from the foodborne outbreaks, highlighting the importance of taking multiple isolates from samples implicated in food poisonings. Notable, in contrast to the soil samples, no bacteria belonging to the psychrotolerant B. cereus group members were detected in the food samples linked to foodborne outbreaks, while the overall abundancy of B. thuringiensis did not significantly differ between the sample categories. None of the isolates was classified as B. cytotoxicus, fostering the hypothesis that the latter species is linked to very specific ecological niches. Overall, our work shows that machine learning assisted (FTIR) spectroscopy is suitable for identification of B. cereus group members in routine diagnostics and outbreak investigations. In addition, it is a promising tool to explore the natural habitats of B. cereus group, such as soil.

Bactéries sporulées et biofilms : un problème récurrent dans les lignes de production de lait reconstitué ou recombiné pasteurisé

In the dairy industry, bacterial contaminants persist on equipment surfaces due to spore and biofilm formation. These are involved in cross-contamination problems that affect the quality of processed products and limit their life. The pasteurization process, in which milk is submitted to moderate heat treatment, is inefficient against bacterial spores. The most prevalent sporulated bacteria belong to Bacillus and related genera. The situation is more complicated in countries where pasteurized milk is derived from imported milk powder originally contaminated by bacterial spores. Studies have shown biofilm formation on dairy equipment by mesophilic strains from the group Bacillus cereus and thermophilic strains from the genus Geobacillus. These biofilms are resistant to cleaning procedures and are sources of chronic contamination of pasteurized milk. This review analyzes the dairy situation in Algeria exposed to sporulated flora and derived biofilm problems, with the aim of proposing efficient solutions in the light of current knowledge. [Journal translation].

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.

LysPBC2, a Novel Endolysin Harboring a Bacillus cereus Spore Binding Domain

To control the spore-forming human pathogen Bacillus cereus, we isolated and characterized a novel endolysin, LysPBC2, from a newly isolated B. cereus phage, PBC2. Compared to the narrow host range of phage PBC2, LysPBC2 showed very broad lytic activity against all Bacillus, Listeria, and Clostridium species tested. In addition to a catalytic domain and a cell wall binding domain, LysPBC2 has a spore binding domain (SBD) partially overlapping its catalytic domain, which specifically binds to B. cereus spores but not to vegetative cells of B. cereus Both immunogold electron microscopy and a binding assay indicated that the SBD binds the external region of the spore cortex layer. Several amino acid residues required for catalytic or spore binding activity of LysPBC2 were determined by mutagenesis studies. Interestingly, LysPBC2 derivatives with impaired spore binding activity showed an increased lytic activity against vegetative cells of B. cereus compared with that of wild-type LysPBC2. Further biochemical studies revealed that these LysPBC2 derivatives have lower thermal stability, suggesting a stabilizing role of SBD in LysPBC2 structure.IMPORTANCE Bacteriophages produce highly evolved lytic enzymes, called endolysins, to lyse peptidoglycan and release their progeny from bacterial cells. Due to their potent lytic activity and specificity, the use of endolysins has gained increasing attention as a natural alternative to antibiotics. Since most endolysins from Gram-positive-bacterium-infecting phages have a modular structure, understanding the function of each domain is crucial to make effective endolysin-based therapeutics. Here, we report the functional and biochemical characterization of a Bacillus cereus phage endolysin, LysPBC2, which has an unusual spore binding domain and a cell wall binding domain. A single point mutation in the spore binding domain greatly enhanced the lytic activity of endolysin at the cost of reduced thermostability. This work contributes to the understanding of the role of each domain in LysPBC2 and will provide insight for the rational design of efficient antimicrobials or diagnostic tools for controlling B. cereus.

Whole-genome-based phylogeny of Bacillus cytotoxicus reveals different clades within the species and provides clues on ecology and evolution

Bacillus cytotoxicus is a member of the Bacillus cereus group linked to fatal cases of diarrheal disease. Information on B. cytotoxicus is very limited; in particular comprehensive genomic data is lacking. Thus, we applied a genomic approach to characterize B. cytotoxicus and decipher its population structure. To this end, complete genomes of ten B. cytotoxicus were sequenced and compared to the four publicly available full B. cytotoxicus genomes and genomes of other B. cereus group members. Average nucleotide identity, core genome, and pan genome clustering resulted in clear distinction of B. cytotoxicus strains from other strains of the B. cereus group. Genomic content analyses showed that a hydroxyphenylalanine operon is present in B. cytotoxicus, but absent in all other members of the B. cereus group. It enables degradation of aromatic compounds to succinate and pyruvate and was likely acquired from another Bacillus species. It allows for utilization of tyrosine and might have given a B. cytotoxicus ancestor an evolutionary advantage resulting in species differentiation. Plasmid content showed that B. cytotoxicus is flexible in exchanging genes, allowing for quick adaptation to the environment. Genome-based phylogenetic analyses divided the B. cytotoxicus strains into four clades that also differed in virulence gene content.

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

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

The cytochrome P450BM-1 of Bacillus megaterium A14K is induced by 2,3,7,8-Tetrachlorinated dibenzo-p-dioxin: Biophysical, molecular and biochemical determinants

The current study describes biological changes in Bacillus megaterium A14K cells growing in the presence of 2,3,7,8-Tetrachlorinated dibenzo-p-dioxin (TCDD), the most potent congener of dioxins. The results indicate that the metabolizing of 2,3,7,8-TCDD by BmA14K was accompanied with a novel morphological and biophysical profile typified by the growth of single cells with high levels of biosurfactant production, surface hydrophobicity and cell membrane permeability. Moreover, the TCDD-grown bacteria exhibited a specific fatty acid profile characterized by low ratios of branched/straight chain fatty acids (BCFAs/SCFAs) and saturated/unsaturated fatty acids (SFAs/USFAs) with a specific "signature" due to the presence of branched chain unsaturated fatty acids (BCUFAs). This was synchronized with a significant induction of P450_BM-1, an unsaturated fatty acid-metabolizing enzyme in B. megaterium. Subsequently, the profile of oxygenated fatty acids in the TCDD-grown bacteria was typified by the presence of 5,6-epoxy derived from unsaturated C15, C16 and C17 fatty acids, that were absent in control bacteria. A net increase was also detected in both hydroxylated and epoxidized fatty acids, especially those derived from C15:0 and C16:1, respectively, suggesting a specific TCDD-induced "signature" of oxygenated fatty acids in BmA14K. Overall, this study sheds light on the use of B. megaterium A14K as a promising bioindicator/biodegrader of dioxins.

Underrecognized niche of spore-forming bacilli as a nitrite-producer isolated from the processing lines and end-products of powdered infant formula

Although nitrite in powdered milk formula (PIF) is a recognized health risk for infants, the presence of nitrite in PIF has only been investigated as a chemical contaminant during the inspection of end-products. The risk posed by microbial sources of nitrite during the PIF manufacturing process has not been considered. This is the first study to report the taxonomy and physiological characteristics of nitrite-producing bacteria isolated from PIF processing environments. All isolates identified as nitrite-producers (133 out of 501 strains collected over four years) from work-in-process and end-products of PIF were spore-forming bacilli. Nitrite-producing metabolism under PIF processing conditions was found in not only thermophilic isolates (3 Bacillus, 60 Geobacillus from 63 strains; 100%) but also in mesophilic isolates (65 Bacillus, 1 Anoxybacillus from 70 strains; 65.7%). Geobacillus was the only highly heat-resistant sporeformer and vigorous nitrite-producer exhibiting dramatic increases in nitrite over short periods of incubation (a maximum value within 3 h). High conversions of nitrate to nitrite (up to 88.8%) was also observed, highlighting bacteria as a key source of nitrite in PIF processing lines. Further research into the diversity of metabolic activity observed in this study can facilitate specialized management of nitrite-producers in PIF processing lines.

Analysis of enterotoxigenic Bacillus cereus strains from dried foods using whole genome sequencing, multi-locus sequence analysis and toxin gene prevalence and distribution using endpoint PCR analysis

Bacillus cereus strains were isolated from dried foods, which included international brands of spices from South East Asia, Mexico and India purchased from several retail stores, samples of powdered infant formula (PIF), medicated fish feed and dietary supplements. The genetic diversity of 64 strains from spices and PIF was determined using a multiplex endpoint PCR assay designed to identify hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Thirteen different B. cereus toxigenic gene patterns or profiles were identified among the strains. Randomly selected B. cereus strains were sequenced and compared with reference Genomic Groups from National Center Biotechnology Information using bioinformatics tools. A comprehensive multi-loci sequence analysis (MLSA) was designed using alleles from 25 known MLST genes specifically tailored for use with whole genome assemblies. A cohort of representative genomes of strains from a few FDA regulated commodities like dry foods and medicated fish feed was used to demonstrate the utility of the 25-MLSA approach for rapid clustering and identification of Genome Groups. The analysis clustered the strains from medicated fish feed, dry foods, and dietary supplements into phylogenetically-related groups. 25-MLSA also pointed to a greater diversity of B. cereus strains from foods and feed than previously recognized. Our integrated approach of toxin gene PCR, and to our knowledge, whole genome sequencing (WGS) based sequence analysis, may be the first of its kind that demonstrates enterotoxigenic potential and genomic diversity in parallel.

Bacillus cereus group: genetic aspects related to food safety and dairy processing

ABSTRACT: Bacillus cereus group includes not pathogenic and high pathogenic species. They are considered as a risk to public health due to foodborne diseases and as an important cause of economic losses to industries due to production of spoilage enzymes. Some researches have been performed in order to assess the possible factors that contribute to put public health into risk because of consumption of food contaminated with viable cells or toxins which have complex mechanisms of production. The control of these bacteria in food is difficult because they are resistant to several processes used in industries. Thus, in this way, this review focused on highlighting the risk due to toxins production by bacteria from B. cereus group in food and the consequences for food safety and dairy industries.

Genomic Characterization and Probiotic Potency of Bacillus sp. DU-106, a Highly Effective Producer of L-Lactic Acid Isolated From Fermented Yogurt

Bacillus sp. DU-106, a newly isolated member of Bacillus cereus group, exhibits the predominant ability to produce L-lactic acid. The probiotic potency of test strain revealed its survivability at acidic pH, bile salts and viability in simulated gastric juice in vitro. The acute oral toxicity test indicated its no toxicity to laboratory mice in vivo. We further determined the complete genome of strain DU-106 to understand genetic basis as a potential probiotic. It has a circular chromosome and three plasmids for a total genome 5,758,208 bp in size with a G + C content of 35.10%. Genes associated with lactate synthesis were found in the DU-106 genome. We also annotated various stress-related, bile salt resistance, and adhesion-related domains in this strain, which likely provide support in exerting probiotic action by enabling adhesion to host epithelial cells and survival under gastrointestinal tract. Moreover, strain DU-106 genome lacks the virulence genes encodes cereulide synthetase, enterotoxin FM, and cytotoxin K. These phenotypic and genomic probiotic potencies facilitate its potential candidate as probiotic starter in food industry.

Novel Bacillus cereus strain from electrokinetically remediated saline soil towards the remediation of crude oil

A new strain SWH-15 was successfully isolated after initial electrokinetic remediation experiment using the same saline soil sampled from Shengli Oilfield, China. Four methods (morphological and biochemical characteristics, whole-cell fatty acid methyl esters (FAMEs) analysis, 16S rRNA sequence analysis and DNA G + C content and DNA-DNA hybridization analysis) were used to identify the taxonomic status of SWH-15 and confirmed that SWH-15 was a novel species of the Bacillus (B.) cereus group. Then, we assessed the degrading ability of the novel strain SWH-15 to crude oil through a microcosm experiment with four treatments, including control (CK), bioremediation using SWH-15 (Bio), electrokinetic remediation (EK), and combined bioremediation and electrokinetic remediation (Bio + EK). The results showed that the Bio + EK combined remediation treatment was more effective than the CK, Bio, and EK treatments in degrading crude oil contaminants. Bioaugmentation, by addition of the strain SWH-15 had synergistic effect with EK in Bio + EK treatment. Bacterial community analysis showed that electrokinetic remediation alone significantly altered the bacterial community of the saline soil. The addition of the strain SWH-15 alone had a weak effect on the bacterial community. However, the strain SWH-15 boosted the growth of other bacterial species in the metabolic network and weakened the impact of electrical field on the whole bacterial community structure in the Bio + EK treatment.

Bacteria from the Bacillus cereus group as contaminants in requeijão curd cheeses and especialidade láctea tipo requeijão

ABSTRACT: This study aimed to evaluate the occurrence of Bacillus cereus group in requeijões and especialidades lácteas tipo requeijão (regular and light) and to verify if there is differences in relation to this occurrence among different categories of these products. A set of 14 (35%) lots was contaminated with this bacterial group from the 40 lots with low counts (maximum 3.1 × 10 CFU/g), and no significant difference regarding counts or presence/absence were observed among the categories of the products. It can be concluded that contamination by B. cereus group in these products is unable to consist in risk to consumers, regarding adequate refrigeration during selling. This study was the first one to report this bacteria group for these dairy products and highlights the needs of further investigations to evaluate the impact of its spoilage during shelf life.

Pathogenicity, population genetics and dissemination of Bacillus anthracis

Bacillus anthracis, the etiological agent of anthrax, procures its particular virulence by a capsule and two AB type toxins: the lethal factor LF and the edema factor EF. These toxins primarily disable immune cells. Both toxins are translocated to the host cell by the adhesin-internalin subunit called protective antigen PA. PA enables LF to reach intra-luminal vesicles, where it remains active for long periods. Subsequently, LF translocates to non-infected cells, leading to inefficient late therapy of anthrax. B. anthracis undergoes slow evolution because it alternates between vegetative and long spore phases. Full genome sequence analysis of a large number of worldwide strains resulted in a robust evolutionary reconstruction of this bacterium, showing that B. anthracis is split in three main clades: A, B and C. Clade A efficiently disseminated worldwide underpinned by human activities including heavy intercontinental trade of goat and sheep hair. Subclade A.Br.WNA, which is widespread in the Northern American continent, is estimated to have split from clade A reaching the Northern American continent in the late Pleistocene epoch via the former Bering Land Bridge and further spread from Northwest southwards. An alternative hypothesis is that subclade A.Br.WNA. evolved from clade A.Br.TEA tracing it back to strains from Northern France that were assumingly dispatched by European explorers that settled along the St. Lawrence River. Clade B established mostly in Europe along the alpine axis where it evolved in association with local cattle breeds and hence displays specific geographic subclusters. Sequencing technologies are also used for forensic applications to trace unintended or criminal acts of release of B. anthracis. Under natural conditions, B. anthracis generally affects domesticated and wild ruminants in arid ecosystems. The more recently discovered B. cereus biovar anthracis spreads in tropical forests, where it threatens particularly endangered primate populations.

Sensitive and Specific Recombinase Polymerase Amplification Assays for Fast Screening, Detection, and Identification of Bacillus anthracis in a Field Setting

Four isothermal recombinase polymerase amplification (RPA) assays were developed for fast in-field identification of Bacillus anthracis The RPA assays targeted three specific sequences (i.e., the BA_5345 chromosomal marker, the lethal factor lef [from pXO1], and the capsule-biosynthesis-related capA [from pXO2]) and a conserved sequence in the adenylate cyclase gene (adk) for the Bacillus cereus group. B. anthracis-specific RPA assays were tested first with purified genomic DNAs (n = 60), including 11 representatives of B. anthracis, and then with soil (n = 8) and white powder (n = 8) samples spiked with inactivated B. anthracis spores and/or other biological agents. The RPA assays were also tested in another laboratory facility, which blindly provided DNA and lysate samples (n = 30, including 20 B. anthracis strains). RPA assays displayed 100% specificity and sensitivity. The hands-off turnaround times at 42°C ranged from 5 to 6 min for 10² genomic copies. The analytical sensitivity of each RPA assay was ∼10 molecules per reaction. In addition, the BA_5345 and adk RPA assays were assessed under field conditions with a series of surface swabs (n = 13, including 11 swabs contaminated with B. thuringiensis spores) that were blindly brought to the field laboratory by a chemical, biological, radiological, and nuclear (CBRN) sampling team. None of the 13 samples, except the control, tested positive for B. anthracis, and all samples that had been harvested from spore-contaminated surfaces tested positive with the adk RPA assay. All three B. anthracis-specific RPA assays proved suitable for rapid and reliable identification of B. anthracis and therefore could easily be used by first responders under field conditions to quickly discriminate between a deliberate release of B. anthracis spores and a hoax attack involving white powder.IMPORTANCE In recent decades, particularly following the 11 September 2001 and Amerithrax attacks, the world has experienced attempts to sow panic and chaos in society through thousands of white-powder copycats using household powders to mimic real bioterrorism attacks. In such circumstances, field-deployable detection methods are particularly needed to screen samples collected from the scene. The aim is to test the samples directly using a fast and reliable assay for detection of the presence of B. anthracis While this would not preclude further confirmatory tests from being performed in reference laboratories, it would bring useful, timely, and relevant information to local crisis managers and help them make appropriate decisions without having to wait for quantitative PCR results (with turnaround times of a few hours) or phenotypic identification and sequencing (with turnaround times of a few days). In the current investigation, we developed a set of isothermal RPA assays for the rapid screening and identification of B. anthracis in powders and soil samples, with the purpose of discriminating a deliberate release of B. anthracis spores from a hoax attack involving white powder; this would also apply to dispersion by spraying of aerosolized forms of B. anthracis Further work is now ongoing to confirm the first observations and validate the on-site use of these assays by first responders.

Genome sequence and comparative analysis of Bacillus cereus BC04, reveals genetic diversity and alterations for antimicrobial resistance

In this study, we delineated the genome sequence of a Bacillus cereus strain BC04 isolated from a stool sample in India. The draft genome is 5.1 Mb in size and consists of total 109 scaffolds, GC content is 35.2% with 5182 coding genes. The comparative analysis with other completely sequenced genomes highlights the unique presence of genomic islands, hemolysin, capsular synthetic protein, modifying enzymes accC7 and catA15, regulators of antibiotic resistance MarR and LysR with annotated functions related to virulence, stress response, and antimicrobial resistance. Overall, this study not only signifies the genetic diversity in gut isolate BC04 in particular, but also pinpoints the presence of unique genes possessed by B. cereus which can be pertinently exploited to design novel drugs and intervention strategies for the treatment of food borne diseases.

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

While some species in the Bacillus cereus group are well-characterized human pathogens (e.g., B. anthracis and B. cereus sensu stricto), the pathogenicity of other species (e.g., B. pseudomycoides) either has not been characterized or is presently not well understood. To provide an updated characterization of the pathogenic potential of species in the B. cereus group, we classified a set of 52 isolates, including 8 type strains and 44 isolates from dairy-associated sources, into 7 phylogenetic clades and characterized them for (i) the presence of toxin genes, (ii) phenotypic characteristics used for identification, and (iii) cytotoxicity to human epithelial cells. Overall, we found that B. cereus toxin genes are broadly distributed but are not consistently present within individual species and/or clades. After growth at 37°C, isolates within a clade did not typically show a consistent cytotoxicity phenotype, except for isolates in clade VI (B. weihenstephanensis/B. mycoides), where none of the isolates were cytotoxic, and isolates in clade I (B. pseudomycoides), which consistently displayed cytotoxic activity. Importantly, our study highlights that B. pseudomycoides is cytotoxic toward human cells. Our results indicate that the detection of toxin genes does not provide a reliable approach to predict the pathogenic potential of B. cereus group isolates, as the presence of toxin genes is not always consistent with cytotoxicity phenotype. Overall, our results suggest that isolates from multiple B. cereus group clades have the potential to cause foodborne illness, although cytotoxicity is not always consistently found among isolates within each clade.IMPORTANCE Despite the importance of the Bacillus cereus group as a foodborne pathogen, characterizations of the pathogenic potential of all B. cereus group species were lacking. We show here that B. pseudomycoides (clade I), which has been considered a harmless environmental microorganism, produces toxins and exhibits a phenotype consistent with the production of pore-forming toxins. Furthermore, B. mycoides/B. weihenstephanensis isolates (clade VI) did not show cytotoxicity when grown at 37°C, despite carrying multiple toxin genes. Overall, we show that the current standard methods to characterize B. cereus group isolates and to detect the presence of toxin genes are not reliable indicators of species, phylogenetic clades, or an isolate's cytotoxic capacity, suggesting that novel methods are still needed for differentiating pathogenic from nonpathogenic species within the B. cereus group. Our results also contribute data that are necessary to facilitate risk assessments and a better understanding as to which B. cereus group species are likely to cause foodborne illness.

Isolation, Identification, Prevalence, and Genetic Diversity of Bacillus cereus Group Bacteria From Different Foodstuffs in Tunisia

Bacillus cereus group is widespread in nature and foods. Several members of this group are recognized as causing food spoilage and/or health issues. This study was designed to determine the prevalence and genetic diversity of the B. cereus group strains isolated in Tunisia from different foods (cereals, spices, cooked food, fresh-cut vegetables, raw and cooked poultry meats, seafood, canned, pastry, and dairy products). In total, 687 different samples were collected and searched for the presence of the B. cereus group after selective plating on MYP agar and enumeration of each sample. The typical pink-orange uniform colonies surrounded by a zone of precipitate were assumed to belong to the B. cereus group. One typical colony from each sample was subcultured and preserved as cryoculture. Overall, 191 (27.8%) food samples were found positive, giving rise to a collection of 191 B. cereus-like isolates. The concentration of B. cereus-like bacteria were below 10³ cfu/g or ml in 77.5% of the tested samples. Higher counts (>10⁴ cfu/g or ml) were found in 6.8% of samples including fresh-cut vegetables, cooked foods, cereals, and pastry products. To verify whether B. cereus-like isolates belonged to the B. cereus group, a PCR test targeting the sspE gene sequence specific of the group was carried out. Therefore, 174 isolates were found to be positive. Food samples were contaminated as follows: cereals (67.6%), pastry products (46.2%), cooked food (40.8%), cooked poultry meat (32.7%), seafood products (32.3%), spices (28.8%), canned products (16.7%), raw poultry meat (9.4%), fresh-cut vegetables (5.0%), and dairy products (4.8%). The 174 B. cereus isolates were characterized by partial sequencing of the panC gene, using a Sym'Previous software tool to assign them to different phylogenetic groups. Strains were distributed as follows: 61.3, 29.5, 7.5, and 1.7% in the group III, IV, II, and V, respectively. The genetic diversity was further assessed by ERIC-PCR and PFGE typing methods. PFGE and ERIC-PCR patterns analysis allowed discriminating 143 and 99 different profiles, respectivey. These findings, associated to a relatively higher prevalence of B. cereus group in different foods, could be a significant etiological agent of food in Tunisia.

Application of a novel amplicon-based sequencing approach reveals the diversity of the Bacillus cereus group in stored raw and pasteurized milk

Members of the Bacillus cereus sensu lato (B. cereus group) are spore-forming organisms commonly associated with spoilage of milk and dairy products. Previous studies have shown, by using 16S marker gene sequencing, that the genus Bacillus is part of the core microbiota of raw bovine milk and that some members of this genus are able to grow during sub-optimal storage (8 °C) of pasteurized consumption milk. Here, the composition of this genus in pasteurized consumption milk samples, collected from two dairies, over a one-year period and stored at 4 or 8 °C up to the end of shelf life is uncovered. Our results show that the B. cereus group is the dominant Bacillus group in stored consumption milk. By applying a new marker gene sequencing approach, several dominating phylogenetic clusters were identified within the B. cereus group populations from the milk samples. There was a higher phylogenetic diversity among bacteria from milk stored at 8 °C compared to milk stored at 4 °C. Sampling period and the dairy the samples were collected from, also significantly influenced the diversity, which shows that the B. cereus group population in consumption milk is heterogeneous and subjected to temporal and spatial changes. The new approach applied in this study will facilitate the identification of isolates within the B. cereus group, of which some are potential spoilage bacteria and pathogenic contaminants of milk and dairy products.

Diversity and enzymatic potentialities of Bacillus sp. strains isolated from a polluted freshwater ecosystem in Cuba

Genotypic and phenotypic characterization of Bacillus spp. from polluted freshwater has been poorly addressed. The objective of this research was to determine the diversity and enzymatic potentialities of Bacillus spp. strains isolated from the Almendares River. Bacilli strains from a polluted river were characterized by considering the production of extracellular enzymes using API ZYM. 14 strains were selected and identified using 16S rRNA, gyrB and aroE genes. Genotypic diversity of the Bacillus spp. strains was evaluated using pulsed field gel electrophoresis. Furthermore, the presence of genetic determinants of potential virulence toxins of the Bacillus cereus group and proteinaceous crystal inclusions of Bacillus thuringiensis was determined. 10 strains were identified as B. thuringiensis, two as Bacillus megaterium, one as Bacillus pumilus and one as Bacillus subtilis. Most strains produced proteases, amylases, phosphatases, esterases, aminopeptidases and glucanases, which reflect the abundance of biopolymeric matter in Almendares River. Comparison of the typing results revealed a spatio-temporal distribution among B. thuringiensis strains along the river. The results of the present study highlight the genotypic and phenotypic diversity of Bacillus spp. strains from a polluted river, which contributes to the knowledge of genetic diversity of Bacilli from tropical polluted freshwater ecosystems.

Unique biomarkers as a potential predictive tool for differentiation of Bacillus cereus group based on real-time PCR

The aim of the study was to develop unique biomarkers for qPCR detection of Bacillus cereus group. Clinical and soil isolates were identified by specifically designed biomarkers - Lipoprotein (OPL-114-lipo), Methyltransferase (MT-17) and S-layer homology domain protein (151-1BC). In order to design biomarkers, we used 120 bacterial strains grouped into B. cereus and non-Bacillus group. The B. cereus group was confirmed by 108 strains of B. cereus and B. thuringiensis (30 reference and 78 wild), along with 3 strains of B. mycoides, B. pseudomycoides, and B. weihenstephanensis; while the non-Bacillus group was composed of 9 Gram-positive and Gram-negative strains. Direct analysis of samples revealed specificity towards identification and characterization of B. cereus group. The newly developed markers OPL-114-lipo and MT-17 showed specificity of 95% and 81%, respectively in identification of B. cereus. They are efficient tools to identify contaminated sources and the degree of bacterial contamination. Environmental and food samples do not require band isolation, re-amplification, sequencing or sequence identification. Thus, reducing the time and cost of analysis. Hence, it will be an alternative approach to traditional culture methods. Commercial food processing industries will be able to employ these biomarkers specific for B. cereus group as a detection tool to reduce economic loss due to B. cereus contamination.

A minireplicon of plasmid pBMB26 represents a new typical replicon in the megaplasmids of Bacillus cereus group

A new minireplicon (rep26 minireplicon) from pBMB26, the 188 kb indigenous plasmid related to spore-crystal association (SCA) phenotype in Bacillus thuringiensis strain YBT-020, was characterized. A 12 kb EcoRI fragment, which encoded 10 putative open reading frames (ORFs), was capable of supporting replication when cloned in a replication probe vector. Deletion and frame shift mutation analysis showed that a 4.1 kb region encompassing two putative ORFs (orf21 and orf22) was essential for the plasmid replication in B. thuringiensis. Gene orf21 encoding a 49.8 kDa protein (named Rep26) with a helix-turn-helix motif showed no homology with known replication proteins and gene orf22 encoding a protein of 82.6 kDa showed homology to bacterial PcrA helicase. The replication origin of rep26 minireplicon was proved to be located in the coding region of orf21. Plasmid stability experiments indicated that the recombinant plasmid containing rep26 minireplicon has excellent segregational stability. BLASTP analysis revealed that amino acid sequences of ORF21 and ORF22 were well conserved among Bacillus cereus group strains. The rep26 minireplicon was widely distributed and could be defined as a new typical replicon in the megaplasmids of B. cereus group.

Subterranean Mammals: Reservoirs of Infection or Overlooked Sentinels of Anthropogenic Environmental Soiling?

Global reports of emergent pathogens in humans have intensified efforts to identify wildlife reservoirs. Subterranean mammals, such as bathyergid mole rats, are largely overlooked, despite their high-level exposure to soil-dwelling microbes. Initial assessment of bathyergid reservoir potential was determined using a broad-range 16S rRNA PCR approach, which revealed an 83% PCR-positivity for the 234 bathyergid lung samples evaluated. The presence of the Bacillus cereus complex, a ubiquitous bacterial assemblage, containing pathogenic and zoonotic species, was confirmed through nucleotide sequencing, prior to group- and species-specific PCR sequencing. The latter allowed for enhanced placement and prevalence estimations of Bacillus in four bathyergid species sampled across a range of transformed landscapes in the Western Cape Province, South Africa. Two novel Bacillus strains (1 and 2) identified on the basis of the concatenated 16S rRNA-groEL-yeaC data set (2066 nucleotides in length), clustered with B. mycoides (ATCC 6462) and B. weihenstephanensis (WSBC 10204), within a well-supported monophyletic lineage. The levels of co-infection, evaluated with a groEL strain-specific assay, developed specifically for this purpose, were high (71%). The overall Bacillus presence of 17.95% (ranging from 0% for Georychus capensis to 45.35% for Bathyergus suillus) differed significantly between host species (χ² = 69.643; df = 3; P < 0.05), being significantly higher in bathyergids sampled near an urban informal settlement (χ² = 70.245; df = 3; P < 0.05). The results highlight the sentinel potential of soil-dwelling mammals for monitoring anthropogenically introduced, opportunistic pathogens and the threats they pose to vulnerable communities, particularly in the developing world.

Spores and vegetative cells of phenotypically and genetically diverse Bacillus cereus sensu lato are common bacteria in fresh water of northeastern Poland

Gram-positive rods Bacillus cereus sensu lato (sl) are common in natural habitats and food products. It is believed that they are restricted to spores; however, their ecology in aquatic habitats is still poorly investigated. Thus, the aim of the study was to assess the rain-dependent fluctuations in the concentration of B. cereus sl vegetative cells and spores, with evaluation of their phylogenetic and population structure in relation to the toxicity and psychrotolerance. We proved that vegetative cells of B. cereus sl are widely distributed in fresh water of rivers and lakes, being as common as spores. Moreover, heavy rain has a huge impact on their concentration in undisturbed environments. The diversity of B. cereus sl reflects the multiple sources of bacteria and the differences between their distinct environments. Next, their diverse genetic structure and phenotypes better fit their ecological properties than their taxonomic affiliation.

Evidence for Complex Formation of the Bacillus cereus Haemolysin BL Components in Solution

Haemolysin BL is an important virulence factor regarding the diarrheal type of food poisoning caused by Bacillus cereus. However, the pathogenic importance of this three-component enterotoxin is difficult to access, as nearly all natural B. cereus culture supernatants additionally contain the highly cytotoxic Nhe, the second three-component toxin involved in the aetiology of B. cereus-induced food-borne diseases. To better address the toxic properties of the Hbl complex, a system for overexpression and purification of functional, cytotoxic, recombinant (r)Hbl components L₂, L₁ and B from E. coli was established and an nheABC deletion mutant was constructed from B. cereus reference strain F837/76. Furthermore, 35 hybridoma cell lines producing monoclonal antibodies (mAbs) against Hbl L₂, L₁ and B were generated. While mAbs 1H9 and 1D8 neutralized Hbl toxicity and thus, represent important tools for future investigations of the mode-of-action of Hbl on the target cell surface, mAb 1D7, in contrast, even enhanced Hbl toxicity by supporting the binding of Hbl B to the cell surface. By using the specific mAbs in Dot blots, indirect and hybrid sandwich enzyme immuno assays (EIAs), complex formation between Hbl L₁ and B, as well as L₁ and L₂ in solution could be shown for the first time. Surface plasmon resonance experiments with the rHbl components confirmed these results with K_D values of 4.7 × 10⁻⁷ M and 1.5 × 10⁻⁷ M, respectively. These findings together with the newly created tools lay the foundation for the detailed elucidation of the molecular mode-of-action of the highly complex three-component Hbl toxin.

Identification and genomic comparison of temperate bacteriophages derived from emetic Bacillus cereus

Cereulide-producing Bacillus cereus isolates can cause serious emetic (vomiting) syndrome and even acute lethality. As mobile genetic elements, the exploration of prophages derived from emetic B. cereus isolates will help in our understanding of the genetic diversity and evolution of these pathogens. In this study, five temperate phages derived from cereulide-producing B. cereus strains were induced, with four of them undergoing genomic sequencing. Sequencing revealed that they all belong to the Siphoviridae family, but presented in different forms in their hosts. PfNC7401 and PfIS075 have typical icosahedral heads, probably existing alone as phagemids in the host with self-replicating capability in the lysogenic state. PfEFR-4, PfEFR-5, and PfATCC7953 have elongated heads, with the genomes of the former two identified as linear dsDNA, which could be integrated into the host genome during the lysogenic state. Genomic comparison of the four phages with others also derived from emetic B. cereus isolates showed similar genome structures and core genes, thus displaying host spectrum specificity. In addition, phylogenic analysis based on the complete genome and conserved tail fiber proteins of 36 Bacillus species-derived phages confirmed that the phages derived from emetic B. cereus strains were highly similar. Furthermore, one endolysin LysPfEFR-4 was cloned and showed lytic activity against all tested emetic B. cereus strains and cross-lytic activity against some other pathogenic bacteria, implying a potential to control bacterial contamination in the food supply.

Complete Genome sequence of the nematicidal Bacillus thuringiensis MYBT18246

10.1601/nm.5000 is a rod-shaped facultative anaerobic spore forming bacterium of the genus 10.1601/nm.4857. The defining feature of the species is the ability to produce parasporal crystal inclusion bodies, consisting of δ-endotoxins, encoded by cry-genes. Here we present the complete annotated genome sequence of the nematicidal 10.1601/nm.5000 strain MYBT18246. The genome comprises one 5,867,749 bp chromosome and 11 plasmids which vary in size from 6330 bp to 150,790 bp. The chromosome contains 6092 protein-coding and 150 RNA genes, including 36 rRNA genes. The plasmids encode 997 proteins and 4 t-RNA's. Analysis of the genome revealed a large number of mobile elements involved in genome plasticity including 11 plasmids and 16 chromosomal prophages. Three different nematicidal toxin genes were identified and classified according to the Cry toxin naming committee as cry13Aa2, cry13Ba1, and cry13Ab1. Strikingly, these genes are located on the chromosome in close proximity to three separate prophages. Moreover, four putative toxin genes of different toxin classes were identified on the plasmids p120510 (Vip-like toxin), p120416 (Cry-like toxin) and p109822 (two Bin-like toxins). A comparative genome analysis of 10.1601/nm.5000 MYBT18246 with three closely related 10.1601/nm.5000 strains enabled determination of the pan-genome of 10.1601/nm.5000 MYBT18246, revealing a large number of singletons, mostly represented by phage genes, morons and cryptic genes.

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.

Proposal of nine novel species of the Bacillus cereus group

Nine novel Gram-stain-positive bacteria were investigated by a polyphasic taxonomic approach. Based on the analysis of 16S rRNA gene sequences, these strains belonged to the Bacillus cereus group, sharing over 97 % similarity with the known species of this group, and less than 95 % similarity with other species of the genus Bacillus. Multilocus sequence typing analysis showed that they formed nine robust and well-separated branches from the known species. The digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values between the nine strains were, respectively, below the 70 and 96 % threshold values for species definition, and between each strain and the known type strains of this group were also below the two threshold values. On the basis of the phenotypic and phylogenetic data, along with low dDDH and ANI values among these strains, these bacteria are assigned to the following nine novel species of the B. cereus group: Bacillus paranthracis sp. nov., type strain Mn5T (=MCCC 1A00395T=KCTC 33714T=LMG 28873T); Bacillus pacificus sp. nov., type strain EB422T (=MCCC 1A06182T=KCTC 33858T); Bacillus tropicus sp. nov., type strain N24T (=MCCC 1A01406T=KCTC 33711T=LMG 28874T); Bacillus albus sp. nov., type strain N35-10-2T (=MCCC 1A02146T=KCTC 33710T=LMG 28875T); Bacillus mobilis sp. nov., type strain 0711P9-1T (=MCCC 1A05942T=KCTC 33717T=LMG 28877T); Bacillus luti sp. nov., type strain TD41T (=MCCC 1A00359T=KCTC 33716T=LMG 28872T); Bacillus proteolyticus sp. nov., type strain TD42T (=MCCC 1A00365T=KCTC 33715T=LMG 28870T); Bacillus nitratireducens sp. nov., type strain 4049T (=MCCC 1A00732T=KCTC 33713T=LMG 28871T); and Bacillus paramycoides sp. nov., type strain NH24A2T (=MCCC 1A04098T=KCTC 33709T=LMG 28876T).

Pan-genome and phylogeny of Bacillus cereus sensu lato

BACKGROUND

Bacillus cereus sensu lato (s. l.) is an ecologically diverse bacterial group of medical and agricultural significance. In this study, I use publicly available genomes and novel bioinformatic workflows to characterize the B. cereus s. l. pan-genome and perform the largest phylogenetic and population genetic analyses of this group to date in terms of the number of genes and taxa included. With these fundamental data in hand, I identify genes associated with particular phenotypic traits (i.e., "pan-GWAS" analysis), and quantify the degree to which taxa sharing common attributes are phylogenetically clustered.

METHODS

A rapid k-mer based approach (Mash) was used to create reduced representations of selected Bacillus genomes, and a fast distance-based phylogenetic analysis of this data (FastME) was performed to determine which species should be included in B. cereus s. l. The complete genomes of eight B. cereus s. l. species were annotated de novo with Prokka, and these annotations were used by Roary to produce the B. cereus s. l. pan-genome. Scoary was used to associate gene presence and absence patterns with various phenotypes. The orthologous protein sequence clusters produced by Roary were filtered and used to build HaMStR databases of gene models that were used in turn to construct phylogenetic data matrices. Phylogenetic analyses used RAxML, DendroPy, ClonalFrameML, PAUP*, and SplitsTree. Bayesian model-based population genetic analysis assigned taxa to clusters using hierBAPS. The genealogical sorting index was used to quantify the phylogenetic clustering of taxa sharing common attributes.

RESULTS

The B. cereus s. l. pan-genome currently consists of ≈60,000 genes, ≈600 of which are "core" (common to at least 99% of taxa sampled). Pan-GWAS analysis revealed genes associated with phenotypes such as isolation source, oxygen requirement, and ability to cause diseases such as anthrax or food poisoning. Extensive phylogenetic analyses using an unprecedented amount of data produced phylogenies that were largely concordant with each other and with previous studies. Phylogenetic support as measured by bootstrap probabilities increased markedly when all suitable pan-genome data was included in phylogenetic analyses, as opposed to when only core genes were used. Bayesian population genetic analysis recommended subdividing the three major clades of B. cereus s. l. into nine clusters. Taxa sharing common traits and species designations exhibited varying degrees of phylogenetic clustering.

CONCLUSIONS

All phylogenetic analyses recapitulated two previously used classification systems, and taxa were consistently assigned to the same major clade and group. By including accessory genes from the pan-genome in the phylogenetic analyses, I produced an exceptionally well-supported phylogeny of 114 complete B. cereus s. l. genomes. The best-performing methods were used to produce a phylogeny of all 498 publicly available B. cereus s. l. genomes, which was in turn used to compare three different classification systems and to test the monophyly status of various B. cereus s. l. species. The majority of the methodology used in this study is generic and could be leveraged to produce pan-genome estimates and similarly robust phylogenetic hypotheses for other bacterial groups.

Characterization of Enterotoxigenic Bacillus cereus sensu lato and Staphylococcus aureus Isolates and Associated Enterotoxin Production Dynamics in Milk or Meat-Based Broth

Bacillus cereus sensu lato species, as well as Staphylococcus aureus, are important pathogenic bacteria which can cause foodborne illness through the production of enterotoxins. This study characterised enterotoxin genes of these species and examined growth and enterotoxin production dynamics of isolates when grown in milk or meat-based broth. All B. cereus s. l. isolates harboured nheA, hblA and entFM toxin genes, with lower prevalence of bceT and hlyII. When grown at 16 °C, toxin production by individual B. cereus s. l. isolates varied depending on the food matrix; toxin was detected at cell densities below 5 log₁₀(CFU/mL). At 16 °C no staphylococcal enterotoxin C (SEC) production was detected by S. aureus isolates, although low levels of SED production was noted. At 30 °C all S. aureus isolates produced detectable enterotoxin in the simulated meat matrix, whereas SEC production was significantly reduced in milk. Relative to B. cereus s. l. toxin production, S. aureus typically required reaching higher cell numbers to produce detectable levels of enterotoxin. Phylogenetic analysis of the sec and sel genes suggested population evolution which correlated with animal host adaptation, with subgroups of bovine isolates or caprine/ovine isolates noted, which were distinct from human isolates. Taken together, this study highlights the marked differences in the production of enterotoxins both associated with different growth matrices themselves, but also in the behaviour of individual strains when exposed to different food matrices.

Surviving Between Hosts: Sporulation and Transmission

To survive adverse conditions, some bacterial species are capable of developing into a cell type, the "spore," which exhibits minimal metabolic activity and remains viable in the presence of multiple environmental challenges. For some pathogenic bacteria, this developmental state serves as a means of survival during transmission from one host to another. Spores are the highly infectious form of these bacteria. Upon entrance into a host, specific signals facilitate germination into metabolically active replicating organisms, resulting in disease pathogenesis. In this article, we will review spore structure and function in well-studied pathogens of two genera, Bacillus and Clostridium, focusing on Bacillus anthracis and Clostridium difficile, and explore current data regarding the lifestyles of these bacteria outside the host and transmission from one host to another.