Multilocus sequence typing reveals thatBacillus cereusstrains isolated from clinical infections have distinct phylogenetic origins

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

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

Bacillus cereus-Attributable Primary Cutaneous Anthrax-Like Infection in Newborn Infants, India

During March 13-June 23, 2018, anthrax-like cutaneous lesions attributed to the Bacillus cereus group of organisms developed in 12 newborns in India. We traced the source of infection to the healthcare kits used for newborn care. We used multilocus sequence typing to characterize the 19 selected strains from various sources in hospital settings, including the healthcare kits. This analysis revealed the existence of a genetically diverse population comprising mostly new sequence types. Phylogenetic analysis clustered most strains into the previously defined clade I, composed primarily of pathogenic bacilli. We suggest that the synergistic interaction of nonhemolytic enterotoxin and sphingomyelinase might have a role in the development of cutaneous lesions. The infection was controlled by removing the healthcare kits and by implementing an ideal housekeeping program. All the newborns recovered after treatment with ciprofloxacin and amikacin.

Genetic diversity and population structure of the Bacillus cereus group bacteria from diverse marine environments

The phylogenetic diversity of marine bacteria belonged to the Bacillus cereus group has not been well investigated. Here, we present the genetic diversity and population structure of 71 bacteria from diverse marine environments, using a multilocus sequence typing (MLST) approach and the analyses of digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) based on some representative genomic sequences. The MLST analysis demonstrated that these isolates were highly diverse and a wide distribution in marine environments and some of them showed niche specificity to some extent. They were assigned to 27 sequence types (STs) with 23 novel STs. Phylogenetic analysis of 82 bacteria containing 11 type strains based on MLST discriminated them as 20 clusters including 10 new ones. Both the dDDH and ANI results supported the proposition that each of 20 clusters represented one independent species, including 10 putative novel species. Values of 98.3% of MLST similarity and 96.2% of ANI were proposed as the standard for the species definition of this group. In summary, the first insight into the phylogenetic diversity of the group bacteria from marine environments will contribute to better understanding of their ecological role and evolution in contrast with terrestrial environments.

Genotypic heterogeneity of emetic toxin producing Bacillus cereus isolates from China

Emetic toxin-producing Bacillus cereus (emetic B. cereus) is the third member of B. cereus group whose toxins are encoded by megaplasmids, beside anthrax and insecticidal toxins of B. anthracis and B. thuringiensis, respectively. A total of 18 emetic isolates collected from food poisoning events, clinical and non-random food samples in Zhejiang province of China, were analyzed by plasmid screening, pulse field gel electrophoresis, multilocus sequence typing, and toxic gene identification to investigate their genotypic diversity. In this study, 13 plasmid profile types, 14 pulse types and 6 different STs from emetic isolates were detected, in which ST 1035,1038,1053,1054 and 1065 were first assigned and reported. The toxic gene ces existed on its own, or coexisted with other toxic genes bceT, cytk, entFM and nhe, but never with hbl in emetic isolates. The results demonstrated that the emetic B. cereus strains from China were heterologous at genotypic level.

MALDI-TOF MS portrait of emetic and non-emetic Bacillus cereus group members

The number of foodborne intoxications caused by emetic Bacillus cereus isolates has increased significantly. As such, rapid and reliable methods to identify emetic strains appear to be clinically relevant. In this study, intact cell matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) was used to differentiate emetic and non-emetic bacilli. The phyloproteomic clustering of 34 B. cereus emetic and 88 non-emetic isolates classified as B. cereus, Bacillus thuringiensis, Bacillus weihenstephanensis, and Bacillus mycoides, showed (i) a clear separation of both groups at a similarity level of 43%, and (ii) a high relatedness among the emetic isolates (similarity of 78%). Specifically, 83 mass peak classes were recognized in the spectral window range between m/z 4000 and 12 000 that were tentatively assigned to 41 protein variants based on a bioinformatic approach. Mass variation between the emetic and the non-emetic subsets was recorded for 27 of them, including ten ribosomal subunit proteins, for which inter-strain polymorphism was confirmed by gene sequencing. Additional peaks were assigned to other proteins such as small acid soluble proteins, cold shock proteins and hypothetical proteins, e.g., carbohydrate kinase. Moreover, the results were supported by in silico analysis of the biomarkers in 259 members of B. cereus group, including Bacillus anthracis, based on their whole-genome sequences. In conclusion, the proteomic profiling by MALDI-TOF MS is a promising and rapid method for pre-screening B. cereus to identify medically relevant isolates and for epidemiologic purposes.

Is Cytotoxin K from Bacillus cereus a bona fide enterotoxin?

Cytotoxin K (CytK) produced by Bacillus cereus s.l. has generally been considered to be associated with the foodborne diarrhoeal syndrome. Two distinct variants of CytK have been reported: CytK-1 from Bacillus cytotoxicus and CytK-2 from B. cereus. In order to determine whether CytK plays a significant role in the diarrhoeal disease, the occurrence of cytK genes was assessed among 390 B. cereus isolates with different origins including clinical and food poisoning samples and was found to be 46%. Interestingly, the cytK occurrence was slightly lower in food poisoning and clinical isolates than in environmental samples. Seventy cytK-2 positive strains (including 28 isolates from foodborne outbreaks) were then selected in order to assess their genetic diversity. A genetic dendrogram based on the cytK-2 sequences of these 70 strains and on two cytK-1 sequences from strains NVH 391-98 and 883-00 showed an important diversity. However, no strain clustering according to the origin or source of isolation was observed. These observations were confirmed by Multi-Locus Sequences Typing (MLST) based on five different loci of housekeeping genes (ccpA, recF, sucC, purF and gdpD) for which no grouping of foodborne outbreak strains could be identified. Therefore, the choice of cytK as virulence factor for the diarrhoeal pathotype does not seem to be relevant per se, even though the involvement of CytK in the diarrhoeal syndrome cannot be fully excluded. Potential synergistic effects between CytK and other virulence factors, together with their potential variable expression levels should be further investigated.

Characteristics and phylogeny of Bacillus cereus strains isolated from Maari, a traditional West African food condiment

Maari is a spontaneously fermented food condiment made from baobab tree seeds in West African countries. This type of product is considered to be safe, being consumed by millions of people on a daily basis. However, due to the spontaneous nature of the fermentation the human pathogen Bacillus cereus occasionally occurs in Maari. This study characterizes succession patterns and pathogenic potential of B. cereus isolated from the raw materials (ash, water from a drilled well (DW) and potash), seed mash throughout fermentation (0-96h), after steam cooking and sun drying (final product) from two production sites of Maari. Aerobic mesophilic bacterial (AMB) counts in raw materials were of 10(5)cfu/ml in DW, and ranged between 6.5×10(3) and 1.2×10(4)cfu/g in potash, 10(9)-10(10)cfu/g in seed mash during fermentation and 10(7) - 10(9) after sun drying. Fifty three out of total 290 AMB isolates were identified as B. cereus sensu lato by use of ITS-PCR and grouped into 3 groups using PCR fingerprinting based on Escherichia coli phage-M13 primer (M13-PCR). As determined by panC gene sequencing, the isolates of B. cereus belonged to PanC types III and IV with potential for high cytotoxicity. Phylogenetic analysis of concatenated sequences of glpF, gmk, ilvD, pta, pur, pycA and tpi revealed that the M13-PCR group 1 isolates were related to B. cereus biovar anthracis CI, while the M13-PCR group 2 isolates were identical to cereulide (emetic toxin) producing B. cereus strains. The M13-PCR group 1 isolates harboured poly-γ-D-glutamic acid capsule biosynthesis genes capA, capB and capC showing 99-100% identity with the environmental B. cereus isolate 03BB108. Presence of cesB of the cereulide synthetase gene cluster was confirmed by PCR in M13-PCR group 2 isolates. The B. cereus harbouring the cap genes were found in potash, DW, cooking water and at 8h fermentation. The "emetic" type B. cereus were present in DW, the seed mash at 48-72h of fermentation and in the final product, while the remaining isolates (PanC type IV) were detected in ash, at 48-72h fermentation and in the final product. This work sheds light on the succession and pathogenic potential of B. cereus species in traditional West African food condiment and clarifies their phylogenetic relatedness to B. cereus biovar anthracis. Future implementation of GMP and HACCP and development of starter cultures for controlled Maari fermentations will help to ensure a safe product.

Identification and validation of specific markers of Bacillus anthracis spores by proteomics and genomics approaches

Bacillus anthracis is the causative bacteria of anthrax, an acute and often fatal disease in humans. The infectious agent, the spore, represents a real bioterrorism threat and its specific identification is crucial. However, because of the high genomic relatedness within the Bacillus cereus group, it is still a real challenge to identify B. anthracis spores confidently. Mass spectrometry-based tools represent a powerful approach to the efficient discovery and identification of such protein markers. Here we undertook comparative proteomics analyses of Bacillus anthracis, cereus and thuringiensis spores to identify proteoforms unique to B. anthracis. The marker discovery pipeline developed combined peptide- and protein-centric approaches using liquid chromatography coupled to tandem mass spectrometry experiments using a high resolution/high mass accuracy LTQ-Orbitrap instrument. By combining these data with those from complementary bioinformatics approaches, we were able to highlight a dozen novel proteins consistently observed across all the investigated B. anthracis spores while being absent in B. cereus/thuringiensis spores. To further demonstrate the relevance of these markers and their strict specificity to B. anthracis, the number of strains studied was extended to 55, by including closely related strains such as B. thuringiensis 9727, and above all the B. cereus biovar anthracis CI, CA strains that possess pXO1- and pXO2-like plasmids. Under these conditions, the combination of proteomics and genomics approaches confirms the pertinence of 11 markers. Genes encoding these 11 markers are located on the chromosome, which provides additional targets complementary to the commonly used plasmid-encoded markers. Last but not least, we also report the development of a targeted liquid chromatography coupled to tandem mass spectrometry method involving the selection reaction monitoring mode for the monitoring of the 4 most suitable protein markers. Within a proof-of-concept study, we demonstrate the value of this approach for the further high throughput and specific detection of B. anthracis spores within complex samples.

Phylogenetic Analysis ofBacillus cereusIsolates from Severe Systemic Infections Using Multilocus Sequence Typing Scheme

Bacillus cereus strains from cases of severe or lethal systemic infections, including respiratory symptoms cases, were analyzed using multilocus sequence typing scheme of B. cereus MLST database. The isolates were evenly distributed between the two main clades, and 60% of them had allele profiles new to the database. Half of the collection's strains clustered in a lineage neighboring Bacillus anthracis phylogenetic origin. Strains from lethal cases with respiratory symptoms were allocated in both main clades. This is the first report of strains causing respiratory symptoms to be identified as genetically distant from B. anthracis. The phylogenetic location of the presented here strains was compared with all previously submitted to the database isolates from systemic infections, and were found to appear in the same clusters where clinical isolates from other studies had been assigned. It seems that the pathogenic strains are forming clusters on the phylogenetic tree.

Extended genetic analysis of Brazilian isolates of Bacillus cereus and Bacillus thuringiensis

Multiple locus sequence typing (MLST) was undertaken to extend the genetic characterization of 29 isolates of Bacillus cereus and Bacillus thuringiensis previously characterized in terms of presence/absence of sequences encoding virulence factors and via variable number tandem repeat (VNTR). Additional analysis involved polymerase chain reaction for the presence of sequences (be, cytK, inA, pag, lef, cya and cap), encoding putative virulence factors, not investigated in the earlier study. MLST analysis ascribed novel and unique sequence types to each of the isolates. A phylogenetic tree was constructed from a single sequence of 2,838 bp of concatenated loci sequences. The strains were not monophyletic by analysis of any specific housekeeping gene or virulence characteristic. No clear association in relation to source of isolation or to genotypic profile based on the presence or absence of putative virulence genes could be identified. Comparison of VNTR profiling with MLST data suggested a correlation between these two methods of genetic analysis. In common with the majority of previous studies, MLST was unable to provide clarification of the basis for pathogenicity among members of the B. cereus complex. Nevertheless, our application of MLST served to reinforce the notion that B. cereus and B. thuringiensis should be considered as the same species.

Contaminated ventilator air flow sensor linked to Bacillus cereus colonization of newborns

We investigated Bacillus cereus–positive tracheal aspirates from infants on ventilators in a neonatal intensive care unit. Multilocus sequence typing determined a genetic match between strains isolated from samples from a case-patient and from the air flow sensor in the ventilator. Changing the sterilization method for sensors to steam autoclaving stopped transmission.

Identification of a Bacillus anthracis specific indel in the yeaC gene and development of a rapid pyrosequencing assay for distinguishing B. anthracis from the B. cereus group

Bacillus anthracis, the causative agent of anthrax, is a potential source of bioterrorism. The existing assays for its identification lack specificity due to the close genetic relationship it exhibits to other members of the B. cereus group. Our comparative analyses of protein sequences from Bacillus species have identified a 24 amino acid deletion in a conserved region of the YeaC protein that is uniquely present in B. anthracis. PCR primers based on conserved regions flanking this indel in the Bacillus cereus group of species (viz. Bacillus cereus, B. anthracis, B. thuringiensis, B. mycoides, B. weihenstephnensis and B. pseudomycoides) specifically amplified a 282 bp fragment from all six reference B. anthracis strains, whereas a 354 bp fragment was amplified from 15 other B. cereus group of species/strains. These fragments, due to large size difference, are readily distinguished by means of agarose gel electrophoresis. In contrast to the B. cereus group, no PCR amplification was observed with any of the non-B. cereus group of species/strains. This indel was also used for developing a rapid pyrosequencing assay for the identification of B. anthracis. Its performance was evaluated by examining the presence or absence of this indel in a panel of 81 B. cereus-like isolates from various sources that included 39 B. anthracis strains. Based upon the sequence data from the pyrograms, the yeaC indel was found to be a distinctive characteristic of various B. anthracis strains tested and not found in any other species/strains from these samples. Therefore, this B. anthracis specific indel provides a robust and highly-specific chromosomal marker for the identification of this high-risk pathogen from other members of the B. cereus group independent of a strain's virulence. The pyrosequencing platform also allows for the rapid and simultaneous screening of multiple samples for the presence of this B. anthracis-specific marker.

Rapidly progressive, fatal, inhalation anthrax-like infection in a human: case report, pathogen genome sequencing, pathology, and coordinated response

CONTEXT

Ten years ago a bioterrorism event involving Bacillus anthracis spores captured the nation's interest, stimulated extensive new research on this pathogen, and heightened concern about illegitimate release of infectious agents. Sporadic reports have described rare, fulminant, and sometimes fatal cases of pneumonia in humans and nonhuman primates caused by strains of Bacillus cereus , a species closely related to Bacillus anthracis.

OBJECTIVES

To describe and investigate a case of rapidly progressive, fatal, anthrax-like pneumonia and the overwhelming infection caused by a Bacillus species of uncertain provenance in a patient residing in rural Texas.

DESIGN

We characterized the genome of the causative strain within days of its recovery from antemortem cultures using next-generation sequencing and performed immunohistochemistry on tissues obtained at autopsy with antibodies directed against virulence proteins of B anthracis and B cereus.

RESULTS

We discovered that the infection was caused by a previously unknown strain of B cereus that was closely related to, but genetically distinct from, B anthracis . The strain contains a plasmid similar to pXO1, a genetic element encoding anthrax toxin and other known virulence factors. Immunohistochemistry demonstrated that several homologs of B anthracis virulence proteins were made in infected tissues, likely contributing to the patient's death.

CONCLUSIONS

Rapid genome sequence analysis permitted us to genetically define this strain, rule out the likelihood of bioterrorism, and contribute effectively to the institutional response to this event. Our experience strongly reinforced the critical value of deploying a well-integrated, anatomic, clinical, and genomic strategy to respond rapidly to a potential emerging, infectious threat to public health.

Extended and global phylogenetic view of the Bacillus cereus group population by combination of MLST, AFLP, and MLEE genotyping data

The Bacillus cereus group of bacteria includes species that can cause food-poisoning or spoilage, such as B. cereus, as well as Bacillus anthracis, the cause of anthrax. In the present report we have conducted a multi-datatype analysis using tools from the HyperCAT database (http://mlstoslo.uio.no/) that we recently developed, combining data from multilocus sequence typing (Tourasse et al., 2010), amplified fragment length polymorphism, and multilocus enzyme electrophoresis typing techniques. We provide a comprehensive snapshot of the B. cereus group population, incorporating 2213 isolates including 450 from food and dairy products, in the form of both phylogenetic supertrees and superclusters of genetically closely related isolates. Our main findings include the detection of phylogenetically separated groups of isolates possibly representing novel evolutionary lineages within the B. cereus group, a putative new branch of B. anthracis, as well as new groups of related strains containing both environmental and clinical isolates. In addition, the multi-datatype analysis revealed to a larger extent than previously recognized that food-borne isolates can share identical genotyping profiles with strains from various other origins. Altogether, the global analysis confirms and extends the results underlining the opportunistic nature of B. cereus group organisms, and the fact that isolates responsible for disease outbreaks and contamination of foodstuffs can originate from various genetic backgrounds.

Rapid multi-locus sequence typing using microfluidic biochips

Background

Multiple locus sequence typing (MLST) has become a central genotyping strategy for analysis of bacterial populations. The scheme involves de novo sequencing of 6–8 housekeeping loci to assign unique sequence types. In this work we adapted MLST to a rapid microfluidics platform in order to enhance speed and reduce laboratory labor time.

Methodology/Principal Findings

Using two integrated microfluidic devices, DNA was purified from 100 Bacillus cereus soil isolates, used as a template for multiplex amplification of 7 loci and sequenced on forward and reverse strands. The time on instrument from loading genomic DNA to generation of electropherograms was only 1.5 hours. We obtained full-length sequence of all seven MLST alleles from 84 representing 46 different Sequence Types. At least one allele could be sequenced from a further 15 strains. The nucleotide diversity of B. cereus isolated in this study from one location in Rockville, Maryland (0.04 substitutions per site) was found to be as great as the global collection of isolates.

Conclusions/Significance

Biogeographical investigation of pathogens is only one of a panoply of possible applications of microfluidics based MLST; others include microbiologic forensics, biothreat identification, and rapid characterization of human clinical samples.

What sets Bacillus anthracis apart from other Bacillus species?

Bacillus anthracis is the cause of anthrax, and two large plasmids are essential for toxicity: pXO1, which contains the toxin genes, and pXO2, which encodes a capsule. B. anthracis forms a highly monomorphic lineage within the B. cereus group, but strains of Bacillus thuringiensis and B. cereus exist that are genetically closely related to the B. anthracis cluster. During the past five years B. cereus strains that contain the pXO1 virulence plasmid were discovered, and strains with both pXO1 and pXO2 have been isolated from great apes in Africa. Therefore, the presence of pXO1 and pXO2 no longer principally separates B. anthracis from other Bacilli. The B. anthracis lineage carries a specific mutation in the global regulator PlcR, which controls the transcription of secreted virulence factors in B. cereus and B. thuringiensis. Coevolution of the B. anthracis chromosome with its plasmids may be the basis for the successful development and uniqueness of the B. anthracis lineage.

Evolution of pathogenicity in the Bacillus cereus group

The Bacillus cereus group of bacteria comprises soil-dwelling saprophytes but on occasion these bacteria can cause a wide range of diseases in humans, including food poisoning, systemic infections and highly lethal forms of anthrax. While anthrax is almost invariably caused by strains from a single evolutionary lineage, Bacillus anthracis, variation in the virulence properties of strains from other lineages has not been fully addressed. Using multi-locus sequence data from 667 strains, we reconstructed the evolutionary history of the B. cereus group in terms of both clonal inheritance and recombination. The strains included 155 clinical isolates representing B. anthracis, and isolates from emetic and diarrhoeal food poisoning, septicaemia and related infections, wound, and lung infections. We confirmed the existence of three major clades and found that clinical isolates of B. cereus (with the exception of emetic toxin-producing strains) are evenly distributed between and within clades 1 and 2. B. anthracis in particular and emetic toxin-producing B. cereus show more clonal structure and are restricted to clade 1. Our characterization of the patterns of genetic exchange showed that there exist partial barriers to gene flow between the three clades. The pathogenic strains do not exhibit atypically high or low rates of recombination, consistent with the opportunistic nature of most pathogenic infections. However, there have been a large number of recent imports in clade 1 of strains from external origins, which is indicative of an on-going shift in gene-flow boundaries for this clade.

Genetic diversity of clinical isolates of Bacillus cereus using multilocus sequence typing

Background

Bacillus cereus is most commonly associated with foodborne illness (diarrheal and emetic) but is also an opportunistic pathogen that can cause severe and fatal infections. Several multilocus sequence typing (MLST) schemes have recently been developed to genotype B. cereus and analysis has suggested a clonal or weakly clonal population structure for B. cereus and its close relatives B. anthracis and B. thuringiensis. In this study we used MLST to determine if B. cereus isolates associated with illnesses of varying severity (e.g., severe, systemic vs. gastrointestinal (GI) illness) were clonal or formed clonal complexes.

Results

A retrospective analysis of 55 clinical B. cereus isolates submitted to the Centers for Disease Control and Prevention between 1954 and 2004 was conducted. Clinical isolates from severe infections (n = 27), gastrointestinal (GI) illness (n = 18), and associated isolates from food (n = 10) were selected for analysis using MLST. The 55 isolates were diverse and comprised 38 sequence types (ST) in two distinct clades. Of the 27 isolates associated with serious illness, 13 clustered in clade 1 while 14 were in clade 2. Isolates associated with GI illness were also found throughout clades 1 and 2, while no isolates in this study belonged to clade 3. All the isolates from this study belonging to the clade 1/cereus III lineage were associated with severe disease while isolates belonging to clade1/cereus II contained isolates primarily associated with severe disease and emetic illness. Only three STs were observed more than once for epidemiologically distinct isolates.

Conclusion

STs of clinical B. cereus isolates were phylogenetically diverse and distributed among two of three previously described clades. Greater numbers of strains will need to be analyzed to confirm if specific lineages or clonal complexes are more likely to contain clinical isolates or be associated with specific illness, similar to B. anthracis and emetic B. cereus isolates.

The Bacillus cereus containing sub-branch most closely related to Bacillus anthracis, have single amino acid substitutions in small acid-soluble proteins, while remaining sub-branches are more variable

Hoffmaster et al. [Hoffmaster AR, Ravel J, Rasko DA, Chapman GD, Chute MD, Marston CK, et al. Identification of anthrax toxin genes in Bacillus cereus associated with illness resembling inhalation anthrax. Proc Natl Acad Sci U S A 2004;101:8449-54; Hoffmaster AR, Hill KK, Gee JE, Marston CK, De BK, Popovic T, et al. Characterization of Bacillus cereus isolates associated with fatal pneumonias: strains are closely related to Bacillus anthracis and harbor B. anthracis virulence genes. J Clin Microbiol 2006;44:3352-60] phylogenetically divided Bacillus cereus strains into 10 branches by amplified fragment length polymorphism (AFLP) with Branch F including all Bacillus anthracis strains and pneumonia-causing strains of B. cereus. There are four sub-branches within Branch F, referred to here as F1-A, F1-B, F2-A and F2-B. The B. anthracis strains are found within sub-branch F1-B. Concerning, the currently available B. cereus pneumonia-causing isolates, one was found to categorize within sub-branch F1-B and two within F2-B. In the following work the sequence variation between B. cereus strains was determined by MALDI-TOF MS and MS-MS for each strain of B. cereus in Branch F. ESI-MS was performed on selected strains for confirmation. Small acid-soluble proteins (SASPs) of B. cereus strains found in F1-B showed a single amino acid substitution, while strains in the other three sub-branches were more variable generally showing one or two amino acid substitutions. The single substitutions always occurred in the C-terminus. Double substitutions occurred in both N and C termini. Of the pneumonia-causing strains, one exhibited a single amino acid substitution, while the other two exhibited a two amino acid substitution.

Multiple-locus sequence typing and analysis of toxin genes in Bacillus cereus food-borne isolates

In the present study we characterized 47 food-borne isolates of Bacillus cereus using multilocus sequence typing (MLST). Newly determined sequences were combined with sequences available in public data banks in order to produce the largest data set possible. Phylogenetic analysis was performed on a total of 296 strains for which MLST sequence information is available, and three main lineages--I, II, and III--within the B. cereus complex were identified. With few exceptions, all food-borne isolates were in group I. The occurrence of horizontal gene transfer (HGT) among various strains was analyzed by several statistical methods, providing evidence of widespread lateral gene transfer within B. cereus. We also investigated the occurrence of toxin-encoding genes, focusing on their evolutionary history within B. cereus. Several patterns were identified, indicating a pivotal role of HGT in the evolution of toxin-encoding genes. Our results indicate that HGT is an important element in shaping the population structure of the B. cereus complex. The results presented here also provide strong evidence of reticulate evolution within the B. cereus complex.

SuperCAT: a supertree database for combined and integrative multilocus sequence typing analysis of the Bacillus cereus group of bacteria (including B. cereus, B. anthracis and B. thuringiensis)

The Bacillus cereus group of bacteria is an important group including mammalian and insect pathogens, such as B. anthracis, the anthrax bacterium, B. thuringiensis, used as a biological pesticide and B. cereus, often involved in food poisoning incidents. To characterize the population structure and epidemiology of these bacteria, five separate multilocus sequence typing (MLST) schemes have been developed, which makes results difficult to compare. Therefore, we have developed a database that compiles and integrates MLST data from all five schemes for the B. cereus group, accessible at http://mlstoslo.uio.no/. Supertree techniques were used to combine the phylogenetic information from analysis of all schemes and datasets, in order to produce an integrated view of the B. cereus group population. The database currently contains strain information and sequence data for 1029 isolates and 26 housekeeping gene fragments, which can be searched by keywords, MLST scheme, or sequence similarity. Supertrees can be browsed according to various criteria such as species, isolate source, or genetic distance, and subtrees containing strains of interest can be extracted. Besides analysis of the available data, the user has the possibility to enter her/his own sequences and compare them to the database and/or include them into the supertree reconstructions.

Cell wall carbohydrate compositions of strains from the Bacillus cereus group of species correlate with phylogenetic relatedness

Members of the Bacillus cereus group contain cell wall carbohydrates that vary in their glycosyl compositions. Recent multilocus sequence typing (MLST) refined the relatedness of B. cereus group members by separating them into clades and lineages. Based on MLST, we selected several B. anthracis, B. cereus, and B. thuringiensis strains and compared their cell wall carbohydrates. The cell walls of different B. anthracis strains (clade 1/Anthracis) were composed of glucose (Glc), galactose (Gal), N-acetyl mannosamine (ManNAc), and N-acetylglucosamine (GlcNAc). In contrast, the cell walls from clade 2 strains (B. cereus type strain ATCC 14579 and B. thuringiensis strains) lacked Gal and contained N-acetylgalactosamine (GalNAc). The B. cereus clade 1 strains had cell walls that were similar in composition to B. anthracis in that they all contained Gal. However, the cell walls from some clade 1 strains also contained GalNAc, which was not present in B. anthracis cell walls. Three recently identified clade 1 strains of B. cereus that caused severe pneumonia, i.e., strains 03BB102, 03BB87, and G9241, had cell wall compositions that closely resembled those of the B. anthracis strains. It was also observed that B. anthracis strains cell wall glycosyl compositions differed from one another in a plasmid-dependent manner. When plasmid pXO2 was absent, the ManNAc/Gal ratio decreased, while the Glc/Gal ratio increased. Also, deletion of atxA, a global regulatory gene, from a pXO2- strain resulted in cell walls with an even greater level of Glc.

Genomic diversity and relationship of Bacillus thuringiensis and Bacillus cereus by multi-REP-PCR fingerprinting

The genomic diversity and relationship among 56 Bacillus thuringiensis and Bacillus cereus type strains were investigated by multi-REP-PCR fingerprinting consisting of three PCR reactions targeting the enterobacterial ERIC1 and ERIC2 and the streptococcal BOXA1R consensus sequences. A total of 113 polymorphic bands were generated in the REP-PCR profiles that allowed tracing of a single dendrogram with three major groups. Bacillus cereus strains clustered together in the A and B groups. Most of the B. thuringiensis strains clustered in group C, which included groups of serovars with a within-group similarity higher than 40% as follows: darmstadiensis, israelensis, and morrisoni; aizawai, kenyae, pakistani, and thompsoni; canadensis, entomocidus, galleriae, kurstaki, and tolworthi; alesti, dendrolimus, and kurstaki; and finitimus, sotto, and thuringiensis. Multi-REP-PCR fingerprinting clustered B. thuringiensis serovars in agreement with previously developed multilocus sequence typing schemes, indicating that it represents a rapid shortcut for addressing the genetic relationship of unknown strains with the major known serovars.

The Bacillus cereus group: novel aspects of population structure and genome dynamics

AIMS

To provide new insights into the population and genomic structure of the Bacillus cereus group of bacteria.

METHODS AND RESULTS

The genetic relatedness among B. cereus group strains was assessed by multilocus sequence typing (MLST) using an optimized scheme based on seven chromosomal housekeeping genes. A set of 48 strains from different clinical sources was included, and six clonal complexes containing several genetically similar isolates from unrelated patients were identified. Interestingly, several clonal groups contained strains that were isolated from similar human sources. Furthermore, comparative whole genome sequence analysis of 16 strains led to the discovery of novel ubiquitous genome features of the B. cereus group, such as atypical group II introns, IStrons, and hitherto uncharacterized repeated elements.

CONCLUSIONS

The B. cereus group constitutes a coherent population unified by the presence of ubiquitous and specific genetic elements which do not show any pattern, either in their sequences or genomic locations, which allows to differentiate between the member species of the group. Nevertheless, the population is very dynamic, as particular lineages of clinical origin can evolve to form clonal complexes. At the genome level, the dynamic behaviour is indicated by the presence of numerous mobile and repeated elements.

SIGNIFICANCE AND IMPACT OF THE STUDY

The B. cereus group of bacteria comprises species that are of medical and economic importance. The MLST data, along with the primers and protocols used, will be available in a public, web-accessible database (http://mlstoslo.uio.no).

Extending the Bacillus cereus group genomics to putative food-borne pathogens of different toxicity

The Bacillus cereus group represents sporulating soil bacteria containing pathogenic strains which may cause diarrheic or emetic food poisoning outbreaks. Multiple locus sequence typing revealed a presence in natural samples of these bacteria of about 30 clonal complexes. Application of genomic methods to this group was however biased due to the major interest for representatives closely related to Bacillus anthracis. Albeit the most important food-borne pathogens were not yet defined, existing data indicate that they are scattered all over the phylogenetic tree. The preliminary analysis of the sequences of three genomes discussed in this paper narrows down the gaps in our knowledge of the B. cereus group. The strain NVH391-98 is a rare but particularly severe food-borne pathogen. Sequencing revealed that the strain should be a representative of a novel bacterial species, for which the name Bacillus cytotoxis or Bacillus cytotoxicus is proposed. This strain has a reduced genome size compared to other B. cereus group strains. Genome analysis revealed absence of sigma B factor and the presence of genes encoding diarrheic Nhe toxin, not detected earlier. The strain B. cereus F837/76 represents a clonal complex close to that of B. anthracis. Including F837/76, three such B. cereus strains had been sequenced. Alignment of genomes suggests that B. anthracis is their common ancestor. Since such strains often emerge from clinical cases, they merit a special attention. The third strain, KBAB4, is a typical facultative psychrophile generally found in soil. Phylogenic studies show that in nature it is the most active group in terms of gene exchange. Genomic sequence revealed high presence of extra-chromosomal genetic material (about 530kb) that may account for this phenomenon. Genes coding Nhe-like toxin were found on a big plasmid in this strain. This may indicate a potential mechanism of toxicity spread from the psychrophile strain community. The results of this genomic work and ecological compartments of different strains incite to consider a necessity of creating prophylactic vaccines against bacteria closely related to NVH391-98 and F837/76. Presumably developing of such vaccines can be based on the properties of non-pathogenic strains such as KBAB4 or ATCC14579 reported here or earlier. By comparing the protein coding genes of strains being sequenced in this project to others we estimate the shared proteome, or core genome, in the B. cereus group to be 3000+/-200 genes and the total proteome, or pan-genome, to be 20-25,000 genes.

Complete sequence analysis of novel plasmids from emetic and periodontal Bacillus cereus isolates reveals a common evolutionary history among the B. cereus-group plasmids, including Bacillus anthracis pXO1

The plasmids of the members of the Bacillus cereus sensu lato group of organisms are essential in defining the phenotypic traits associated with pathogenesis and ecology. For example, Bacillus anthracis contains two plasmids, pXO1 and pXO2, encoding toxin production and encapsulation, respectively, that define this species pathogenic potential, whereas the presence of a Bt toxin-encoding plasmid defines Bacillus thuringiensis isolates. In this study the plasmids from B. cereus isolates that produce emetic toxin or are linked to periodontal disease were sequenced and analyzed. Two periodontal isolates examined contained almost identical approximately 272-kb plasmids, named pPER272. The emetic toxin-producing isolate contained one approximately 270-kb plasmid, named pCER270, encoding the cereulide biosynthesis gene cluster. Comparative sequence analyses of these B. cereus plasmids revealed a high degree of sequence similarity to the B. anthracis pXO1 plasmid, especially in a putative replication region. These plasmids form a newly defined group of pXO1-like plasmids. However, these novel plasmids do not contain the pXO1 pathogenicity island, which in each instance is replaced by plasmid specific DNA. Plasmids pCER270 and pPER272 share regions that are not found in any other pXO1-like plasmids. Evolutionary studies suggest that these plasmids are more closely related to each other than to other identified B. cereus plasmids. Screening of a population of B. cereus group isolates revealed that pXO1-like plasmids are more often found in association with clinical isolates. This study demonstrates that the pXO1-like plasmids may define pathogenic B. cereus isolates in the same way that pXO1 and pXO2 define the B. anthracis species.

Molecular approaches to identify and differentiate Bacillus anthracis from phenotypically similar Bacillus species isolates

Background

Bacillus anthracis and Bacillus cereus can usually be distinguished by standard microbiological methods (e.g., motility, hemolysis, penicillin susceptibility and susceptibility to gamma phage) and PCR. However, we have identified 23 Bacillus spp. isolates that gave discrepant results when assayed by standard microbiological methods and PCR. We used multiple-locus variable-number tandem repeat analysis (MLVA), multiple-locus sequence typing (MLST), and phenotypic analysis to characterize these isolates, determine if they cluster phylogenetically and establish whether standard microbiological identification or PCR were associated with false positive/negative results.

Results

Six isolates were LRN real-time PCR-positive but resistant to gamma phage; MLVA data supported the identification of these isolates as gamma phage-resistant B. anthracis. Seventeen isolates were LRN real-time PCR-negative but susceptible to gamma phage lysis; these isolates appear to be a group of unusual gamma phage-susceptible B. cereus isolates that are closely related to each other and to B. anthracis. All six B. anthracis MLVA chromosomal loci were amplified from one unusual gamma phage-susceptible, motile, B. cereus isolate (although the amplicons were atypical sizes), and when analyzed phylogenetically, clustered with B. anthracis by MLST.

Conclusion

MLVA and MLST aided in the identification of these isolates when standard microbiological methods and PCR could not definitely identify or rule out B. anthracis. This study emphasized the need to perform multiple tests when attempting to identify B. anthracis since relying on a single assay remains problematic due to the diverse nature of bacteria.