Draft Whole-Genome Sequences of 11 Bacillus cereus Food Isolates

Cereulide Synthetase Acquisition and Loss Events within the Evolutionary History of Group III Bacillus cereus Sensu Lato Facilitate the Transition between Emetic and Diarrheal Foodborne Pathogens

Cereulide-producing members of Bacillus cereussensu lato group III (also known as emetic B. cereus) possess cereulide synthetase, a plasmid-encoded, nonribosomal peptide synthetase encoded by the ces gene cluster. Despite the documented risks that cereulide-producing strains pose to public health, the level of genomic diversity encompassed by emetic B. cereus has never been evaluated at a whole-genome scale. Here, we employ a phylogenomic approach to characterize group III B. cereussensu lato genomes which possess ces (ces positive) alongside their closely related, ces-negative counterparts (i) to assess the genomic diversity encompassed by emetic B. cereus and (ii) to identify potential ces loss and/or gain events within the evolutionary history of the high-risk and medically relevant sequence type (ST) 26 lineage often associated with emetic foodborne illness. Using all publicly available ces-positive group III B. cereussensu lato genomes and the ces-negative genomes interspersed among them (n = 159), we show that emetic B. cereus is not clonal; rather, multiple lineages within group III harbor cereulide-producing strains, all of which share an ancestor incapable of producing cereulide (posterior probability = 0.86 to 0.89). Members of ST 26 share an ancestor that existed circa 1748 (95% highest posterior density [HPD] interval = 1246.89 to 1915.64) and first acquired the ability to produce cereulide before 1876 (95% HPD = 1641.43 to 1946.70). Within ST 26 alone, two subsequent ces gain events were observed, as well as three ces loss events, including among isolates responsible for B. cereussensu lato toxicoinfection (i.e., "diarrheal" illness).IMPORTANCEB. cereus is responsible for thousands of cases of foodborne disease each year worldwide, causing two distinct forms of illness: (i) intoxication via cereulide (i.e., emetic syndrome) or (ii) toxicoinfection via multiple enterotoxins (i.e., diarrheal syndrome). Here, we show that emetic B. cereus is not a clonal, homogenous unit that resulted from a single cereulide synthetase gain event followed by subsequent proliferation; rather, cereulide synthetase acquisition and loss is a dynamic, ongoing process that occurs across lineages, allowing some group III B. cereussensu lato populations to oscillate between diarrheal and emetic foodborne pathogens over the course of their evolutionary histories. We also highlight the care that must be taken when selecting a reference genome for whole-genome sequencing-based investigation of emetic B. cereussensu lato outbreaks, since some reference genome selections can lead to a confounding loss of resolution and potentially hinder epidemiological investigations.

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.

Draft Genome Sequences of Enterotoxigenic Bacillus cereus Strains Obtained from Powdered Infant Formula

We introduce the draft genome sequences of five enterotoxigenic Bacillus cereus strains: Bc 12, Bc 67, Bc 111, Bc 112, and Bc 113, which were obtained from powdered infant formula. The genome sizes of the strains ranged from 5.5 to 5.8 Mb, and the G+C contents were ~35.2%.

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.

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.

Linking Bacillus cereus Genotypes and Carbohydrate Utilization Capacity

We characterised carbohydrate utilisation of 20 newly sequenced Bacillus cereus strains isolated from food products and food processing environments and two laboratory strains, B. cereus ATCC 10987 and B. cereus ATCC 14579. Subsequently, genome sequences of these strains were analysed together with 11 additional B. cereus reference genomes to provide an overview of the different types of carbohydrate transporters and utilization systems found in B. cereus strains. The combined application of API tests, defined growth media experiments and comparative genomics enabled us to link the carbohydrate utilisation capacity of 22 B. cereus strains with their genome content and in some cases to the panC phylogenetic grouping. A core set of carbohydrates including glucose, fructose, maltose, trehalose, N-acetyl-glucosamine, and ribose could be used by all strains, whereas utilisation of other carbohydrates like xylose, galactose, and lactose, and typical host-derived carbohydrates such as fucose, mannose, N-acetyl-galactosamine and inositol is limited to a subset of strains. Finally, the roles of selected carbohydrate transporters and utilisation systems in specific niches such as soil, foods and the human host are discussed.