A comparative study of enzyme variation in Bacillus cereus and Bacillus thuringiensis

Thirty-two strains of Bacillus spp. were examined in a multilocus enzyme study by agarose gel electrophoresis. The organisms were Bacillus thuringiensis (21 strains, B. cereus (8), including two of var. mycoides, and B. megaterium (3). Strains having similar enzyme variants were grouped into zymovars. A total of 10 of 11 enzyme loci studied were polymorphic and 27 zymovars were distinguished among the 32 strains. The results were subjected to numerical analysis, phenetic affinities and genetic distances between the strains were calculated. The numerical analysis was unable to differentiate between B. thuringiensis and B. cereus. Our results indicated that based on this multilocus enzyme study these zymovars should be considered as belonging to the same species. A mycoides variant of B. cereus was the most distinctive strain studied and clearly belonged to a separate species, B. mycoides. The technique also allowed for identification of contamination and mislabelling of strains.     

A randomly amplified polymorphic DNA marker specific for the Bacillus cereus group is diagnostic for Bacillus anthracis

Aiming to develop a DNA marker specific for Bacillus anthracis and able to discriminate this species from Bacillus cereus, Bacillus thuringiensis, and Bacillus mycoides, we applied the randomly amplified polymorphic DNA (RAPD) fingerprinting technique to a collection of 101 strains of the genus Bacillus, including 61 strains of the B. cereus group. An 838-bp RAPD marker (SG-850) specific for B. cereus, B. thuringiensis, B. anthracis, and B. mycoides was identified. This fragment included a putative (366-nucleotide) open reading frame highly homologous to the ypuA gene of Bacillus subtilis. The restriction analysis of the SG-850 fragment with AluI distinguished B. anthracis from the other species of the B. cereus group.     

The hemolytic enterotoxin HBL is broadly distributed among species of the Bacillus cereus group

The prevalence of the hemolytic enterotoxin complex HBL was determined in all species of the Bacillus cereus group with the exception of Bacillus anthracis. hblA, encoding the binding subunit B, was detected by PCR and Southern analysis and was confirmed by partial sequencing of 18 strains. The sequences formed two clusters, one including B. cereus and Bacillus thuringiensis strains and the other one consisting of Bacillus mycoides, Bacillus pseudomycoides, and Bacillus weihenstephanensis strains. From eight B. thuringiensis strains, the enterotoxin gene hblA could be amplified. Seven of them also expressed the complete HBL complex as determined with specific antibodies against the L(1), L(2), and B components. Eleven of 16 B. mycoides strains, all 3 B. pseudomyoides strains, 9 of 15 B. weihenstephanensis strains, and 10 of 23 B. cereus strains carried hblA. While HBL was not expressed in the B. pseudomycoides strains, the molecular assays were in accordance with the immunological assays for the majority of the remaining strains. In summary, the hemolytic enterotoxin HBL seems to be broadly distributed among strains of the B. cereus group and relates neither to a certain species nor to a specific environment. The consequences of this finding for food safety considerations need to be evaluated.     

The autolytic phenotype of the Bacillus cereus group

AIM: To determine the autolytic phenotype of five species in the Bacillus cereus group. METHODS AND RESULTS: The autolytic rate of 96 strains belonging to five species in the B. cereus group was examined under starvation conditions at pH 6, 6.5 and 8.5 in different buffers. The autolytic rate was strain-dependent with a wide variability at pH 6, but higher and more uniform at pH 6.5. At pH 8.5, and respect to the extent of autolysis at pH 6.5, it was relatively low for most of the strains with the lowest values between 13 and 52% in Bacillus mycoides and Bacillus pseudomycoides. Peptidoglycan hydrolase patterns evaluated by renaturing sodium dodecyl sulfate-polyacrylamide gel electrophoresis using cells of Bacillus thuringiensis ssp. tolworthi HD125 as an indicator, revealed complex profiles with lytic bands of about 90, 63, 46, 41, 38, 32, 28 and 25 kDa in B. cereus, B. thuringiensis and Bacillus weihenstephanensis. Bacillus mycoides and B. pseudomycoides had simpler profiles with lytic bands of 63, 46 and 38 kDa. Changes in the autolytic pattern were observed for cells harvested at the stationary phase of growth (72 h) showing an increase in the intensity of the 25 kDa band in the case of B. cereus, B. thuringiensis and B. weihenstephanensis, while no changes were observed for B. mycoides. Using Micrococcus lysodeicticus and Listeria monocytogenes as indicators lytic activity was retained by proteins of 63, 46, 38, 32 and 25 kDa and a new one of about 20 kDa in B. mycoides. Growth in the different media did not affect the autolytic pattern. NaCl abolished the activity of all the peptidoglycan hydrolases except for those of B. mycoides and B. weihenstephanensis. Lytic activity was retained in the presence of MgCl(2), MnCl(2) and EDTA and increased at basic pH. CONCLUSIONS: Bacillus cereus/B. thuringiensis/B. weihenstephanensis showed a high extent of autolysis around neutral pH, even though they presented relatively complex autolysin profiles at alkaline pH. Bacillus mycoides/B. pseudomycoides had a higher extent of autolysis at acidic pH and a simpler autolysin pattern. SIGNIFICANCE AND IMPACT OF THE STUDY: Information on the autolytic phenotype expand the phenotypic characterization of the different species in the B. cereus group.     

Homoduplex and heteroduplex polymorphisms of the amplified ribosomal 16S-23S internal transcribed spacers describe genetic relationships in the "Bacillus cereus group"

Bacillus anthracis, Bacillus cereus, Bacillus mycoides, Bacillus pseudomycoides, Bacillus thuringiensis, and Bacillus weihenstephanensis are closely related in phenotype and genotype, and their genetic relationship is still open to debate. The present work uses amplified 16S-23S internal transcribed spacers (ITS) to discriminate between the strains and species and to describe the genetic relationships within the "B. cereus group," advantage being taken of homoduplex-heteroduplex polymorphisms (HHP) resolved by polyacrylamide gel electrophoresis and silver staining. One hundred forty-one strains belonging to the six species were investigated, and 73 ITS-HHP pattern types were distinguished by MDE, a polyacrylamide matrix specifically designed to resolve heteroduplex and single-strand conformation polymorphisms. The discriminating bands were confirmed as ITS by Southern hybridization, and the homoduplex or heteroduplex nature was identified by single-stranded DNA mung bean nuclease digestion. Several of the ITS-HHP types corresponded to specific phenotypes such as B. anthracis or serotypes of B. thuringiensis. Unweighted pair group method arithmetic average cluster analysis revealed two main groups. One included B. mycoides, B. weihenstephanensis, and B. pseudomycoides. The second included B. cereus and B. thuringiensis, B. anthracis appeared as a lineage of B. cereus.     

Identification of quorum-quenching N-acyl homoserine lactonases from Bacillus species

A range of gram-negative bacterial species use N-acyl homoserine lactone (AHL) molecules as quorum-sensing signals to regulate different biological functions, including production of virulence factors. AHL is also known as an autoinducer. An autoinducer inactivation gene, aiiA, coding for an AHL lactonase, was cloned from a bacterial isolate, Bacillus sp. strain 240B1. Here we report identification of more than 20 bacterial isolates capable of enzymatic inactivation of AHLs from different sources. Eight isolates showing strong AHL-inactivating enzyme activity were selected for a preliminary taxonomic analysis. Morphological phenotypes and 16S ribosomal DNA sequence analysis indicated that these isolates probably belong to the species Bacillus thuringiensis. Enzymatic analysis with known Bacillus strains confirmed that all of the strains of B. thuringiensis and the closely related species B. cereus and B. mycoides tested produced AHL-inactivating enzymes but B. fusiformis and B. sphaericus strains did not. Nine genes coding for AHL inactivation were cloned either by functional cloning or by a PCR procedure from selected bacterial isolates and strains. Sequence comparison of the gene products and motif analysis showed that the gene products belong to the same family of AHL lactonases.     

Population structure and evolution of the Bacillus cereus group

Representative strains of the Bacillus cereus group of bacteria, including Bacillus anthracis (11 isolates), B. cereus (38 isolates), Bacillus mycoides (1 isolate), Bacillus thuringiensis (53 isolates from 17 serovars), and Bacillus weihenstephanensis (2 isolates) were assigned to 59 sequence types (STs) derived from the nucleotide sequences of seven alleles, glpF, gmk, ilvD, pta, pur, pycA, and tpi. Comparisons of the maximum likelihood (ML) tree of the concatenated sequences with individual gene trees showed more congruence than expected by chance, indicating a generally clonal structure to the population. The STs followed two major lines of descent. Clade 1 comprised B. anthracis strains, numerous B. cereus strains, and rare B. thuringiensis strains, while clade 2 included the majority of the B. thuringiensis strains together with some B. cereus strains. Other species were allocated to a third, heterogeneous clade. The ML trees and split decomposition analysis were used to assign STs to eight lineages within clades 1 and 2. These lineages were defined by bootstrap analysis and by a preponderance of fixed differences over shared polymorphisms among the STs. Lineages were named with reference to existing designations: Anthracis, Cereus I, Cereus II, Cereus III, Kurstaki, Sotto, Thuringiensis, and Tolworthi. Strains from some B. thuringiensis serovars were wholly or largely assigned to a single ST, for example, serovar aizawai isolates were assigned to ST-15, serovar kenyae isolates were assigned to ST-13, and serovar tolworthi isolates were assigned to ST-23, while other serovars, such as serovar canadensis, were genetically heterogeneous. We suggest a revision of the nomenclature in which the lineage and clone are recognized through name and ST designations in accordance with the clonal structure of the population.     

Improved cytotoxicity assay for Bacillus cereus diarrhoeal enterotoxin

An improved McCoy cell cytotoxicity assay for Bacillus cereus diarrhoeal toxin, which includes a staining procedure in addition to visual examination, was developed and the method was compared with two commercially available kits (Oxoid BCET-RPLA and Tecra BDE-VIA). A total of 71 strains of 15 different Bacillus, Brevibacillus and Paenibacillus species, including 16 strains of B. cereus from outbreaks of food-borne illness, were evaluated for toxin production. Eleven of the outbreak strains exhibited cytotoxicity, including all six B. cereus strains associated with diarrhoeal illness. Several other isolates of B. cereus, and its relatives B. anthracis, B. mycoides and B. thuringiensis, exhibited similar cytotoxicity. The other species showed no cytotoxicity, with the exception of certain B. subtilis strains. The cytotoxicity assay was more sensitive than the Oxoid kit and unlike the Tecra kit, did not give false positive results with supernatant fluids heat-treated to destroy the toxin.     

PCR assay of the groEL gene for detection and differentiation of Bacillus cereus group cells

Strains of species in the Bacillus cereus group are potentially enterotoxic. Thus, the detection of all B. cereus group strains is important. As 16S ribosomal DNA sequence analysis cannot adequately differentiate species of the B. cereus group, we explored the potential of the groEL gene as a phylogenetic marker. A phylogenetic analysis of the groEL sequences of 78 B. cereus group strains revealed that the B. cereus group strains were split into two major clusters, one including six B. mycoides and one B. pseudomycoides (cluster II) and the other including two B. mycoides and the rest of the B. cereus group strains (cluster I). Cluster I was further differentiated into two subclusters, Ia and Ib. The sodA gene sequences of representative strains from different clusters were also compared. The phylogenetic tree constructed from the sodA sequences showed substantial similarity to the tree constructed from the groEL sequences. Based on the groEL sequences, a PCR assay for detection and identification of B. cereus group strains was developed. Subsequent restriction fragment length polymorphism (RFLP) analysis verified the PCR amplicons and the differentiation of the B. cereus group strains. RFLP with MboI was identical for all the B. cereus group strains analyzed, while RFLP with MfeI or PstI classified all B. cereus and B. thuringiensis strains into two groups. All cluster II B. mycoides and B. pseudomycoides strains could be discriminated from other B. cereus group bacteria by restriction analysis with TspRI.     

Differentiation of dairy strains of the Bacillus cereus group by phage typing, minimum growth temperature, and fatty acid analysis

A total of 130 Bacillus strains were isolated from dairy products, the dairy environment and from packaging boards and board-producing machines. Ninety-eight of these were members of the B. cereus group ( B. cereus, B. mycoides and B. thuringiensis ) as determined by whole cell fatty acid composition. Fatty acid composition did not differentiate between the three species. Of the 98 strains, which were indistinguishable by biochemical tests, 87 could be assigned into 21 different phage types (11 strains remained untypable) when tested with 12 B. cereus, B. mycoides and B. thuringiensis phages. The distribution of phage types between strains from different sources showed that the source of contamination of the dairy products was of milk origin and not from the packaging materials. Most strains isolated from the dairy products were able to grow below 10°C, whereas strains from the dairy environment and from board mills had higher minimum growth temperatures.     

Sequence diversity of the Bacillus thuringiensis and B. cereus sensu lato flagellin (H antigen) protein: comparison with H serotype diversity

We set out to analyze the sequence diversity of the Bacillus thuringiensis flagellin (H antigen [Hag]) protein and compare it with H serotype diversity. Some other Bacillus cereus sensu lato species and strains were added for comparison. The internal sequences of the flagellin (hag) alleles from 80 Bacillus thuringiensis strains and 16 strains from the B. cereus sensu lato group were amplified and cloned, and their nucleotide sequences were determined and translated into amino acids. The flagellin allele nucleotide sequences for 10 additional strains were retrieved from GenBank for a total of 106 Bacillus species and strains used in this study. These included 82 B. thuringiensis strains from 67 H serotypes, 5 B. cereus strains, 3 Bacillus anthracis strains, 3 Bacillus mycoides strains, 11 Bacillus weihenstephanensis strains, 1 Bacillus halodurans strain, and 1 Bacillus subtilis strain. The first 111 and the last 66 amino acids were conserved. They were referred to as the C1 and C2 regions, respectively. The central region, however, was highly variable and is referred to as the V region. Two bootstrapped neighbor-joining trees were generated: a first one from the alignment of the translated amino acid sequences of the amplified internal sequences of the hag alleles and a second one from the alignment of the V region amino acid sequences, respectively. Of the eight clusters revealed in the tree inferred from the entire C1-V-C2 region amino acid sequences, seven were present in corresponding clusters in the tree inferred from the V region amino acid sequences. With regard to B. thuringiensis, in most cases, different serovars had different flagellin amino acid sequences, as might have been expected. Surprisingly, however, some different B. thuringiensis serovars shared identical flagellin amino acid sequences. Likewise, serovars from the same H serotypes were most often found clustered together, with exceptions. Indeed, some serovars from the same H serotype carried flagellins with sufficiently different amino acid sequences as to be located on distant clusters. Species-wise, B. halodurans, B. subtilis, and B. anthracis formed specific branches, whereas the other four species, all in the B. cereus sensu lato group, B. mycoides, B. weihenstephanensis, B. cereus, and B. thuringiensis, did not form four specific clusters as might have been expected. Rather, strains from any of these four species were placed side by side with strains from the other species. In the B. cereus sensu lato group, B. anthracis excepted, the distribution of strains was not species specific.     

Differentiation of dairy strains of the Bacillus cereus group by phage typing, minimum growth temperature, and fatty acid analysis

A total of 130 Bacillus strains were isolated from dairy products, the dairy environment and from packaging boards and board-producing machines. Ninety-eight of these were members of the B. cereus group (B. cereus, B. mycoides and B. thuringiensis) as determined by whole cell fatty acid composition. Fatty acid composition did not differentiate between the three species. Of the 98 strains, which were indistinguishable by biochemical tests, 87 could be assigned into 21 different phage types (11 strains remained untypable) when tested with 12 B. cereus, B. mycoides and B. thuringiensis phages. The distribution of phage types between strains from different sources showed that the source of contamination of the dairy products was of milk origin and not from the packaging materials. Most strains isolated from the dairy products were able to grow below 10 degrees C, whereas strains from the dairy environment and from board mills had higher minimum growth temperatures.     

Cloning and nucleotide sequence analysis of gyrB of Bacillus cereus, B. thuringiensis, B. mycoides, and B. anthracis and their application to the detection of B. cereus in rice

As 16S rRNA sequence analysis has proven inadequate for the differentiation of Bacillus cereus from closely related species, we employed the gyrase B gene (gyrB) as a molecular diagnostic marker. The gyrB genes of B. cereus JCM 2152(T), Bacillus thuringiensis IAM 12077(T), Bacillus mycoides ATCC 6462(T), and Bacillus anthracis Pasteur #2H were cloned and sequenced. Oligonucleotide PCR primer sets were designed from within gyrB sequences of the respective bacteria for the specific amplification and differentiation of B. cereus, B. thuringiensis, and B. anthracis. The results from the amplification of gyrB sequences correlated well with results obtained with the 16S rDNA-based hybridization study but not with the results of their phenotypic characterization. Some of the reference strains of both B. cereus (three serovars) and B. thuringiensis (two serovars) were not positive in PCR amplification assays with gyrB primers. However, complete sequencing of 1.2-kb gyrB fragments of these reference strains showed that these serovars had, in fact, lower homology than their originally designated species. We developed and tested a procedure for the specific detection of the target organism in boiled rice that entailed 15 h of preenrichment followed by PCR amplification of the B. cereus-specific fragment. This method enabled us to detect an initial inoculum of 0.24 CFU of B. cereus cells per g of boiled rice food homogenate without extracting DNA. However, a simple two-step filtration step is required to remove PCR inhibitory substances.     

Application of the multiplex PCR and PCR-RFLP method in the identification of the Bacillus anthracis

The aim of this study was to apply the multiplex PCR and PCR-RFLP method for the identification of the B. anthracis strains and to distinguish those bacteria from other members of the Bacillus cereus group. The multiplex PCR method enables to detect the virulence factors, i.e. the toxin and the capsule in B. anthracis strains. To do that, the authors have used 5 primer pairs specific for the fragments of lef, cya, pag genes which are present in the pXO1 plasmid and encode the toxin, the cap gene, which is present in the pXO2 plasmid and encodes the capsule, and the Ba813 chromosomal sequence. Among the four B. anthracis strains examined, three contained two plasmids and the Ba813 chromosomal sequence, while the fourth one contained the pXO1 plasmid only, together and the Ba813 chromosomal sequence. Other bacterial species, belonging to the B. cereus group, were also examined: 6 strains of B. cereus, 4 strains of B. thuringiensis and one strain of B. mycoides. The presence of Ba813 chromosomal sequence has been detected in two B. cereus strains. Neither plasmids nor Ba813 chromosomal sequence have been discovered in other B. cereus, B. thuringiensis and B. mycoides strains. The results of the survey indicate that the Ba813 chromosomal sequence does not occur solely in B. anthracis strains. The PCR-RFLP method with the use of SG-749f and SG-749r primers enabled to demonstrate the presence of DNA sequence (SG-749) in B. anthracis, B. cereus, B. thuringiensis and B. mycoides strains. Restriction analysis with enzyme AluI of the SG-749 sequence, has shown the presence of two DNA fragments at the size of about 90 and 660 bp in all B. anthracis strains. The restriction profile obtained was characteristic for B. anthracis strains and it did not occur in other investigated bacterial species belonging to the B. cereus group. It was not observed even in such B. cereus strains in which the presence of Ba813 sequence was discovered and it enabled to differentiate between B. anthracis strains and other closely related species of the B. cereus group.     

Rapid characterization of spores of Bacillus cereus group bacteria by matrix-assisted laser desorption-ionization time-of-flight mass spectrometry

Matrix-assisted laser desorption-ionization (MALDI) time-of-flight mass spectrometry was used to characterize the spores of 14 microorganisms of the Bacillus cereus group. This group includes the four Bacillus species B. anthracis, B. cereus, B. mycoides, and B. thuringiensis. MALDI mass spectra obtained from whole bacterial spores showed many similarities between the species, except for B. mycoides. At the same time, unique mass spectra could be obtained for the different B. cereus and B. thuringiensis strains, allowing for differentiation at the strain level. To increase the number of detectable biomarkers in the usually peak-poor MALDI spectra of spores, the spores were treated by corona plasma discharge (CPD) or sonicated prior to MALDI analysis. Spectra of sonicated or CPD-treated spores displayed an ensemble of biomarkers common for B. cereus group bacteria. Based on the spectra available, these biomarkers differentiate B. cereus group spores from those of Bacillus subtilis and Bacillus globigii. The effect of growth medium on MALDI spectra of spores was also explored.     

Efficacy of a Fluorescent-Antibody Procedure for Identifying Bacillus cereus in Foods

One hundred and seventeen strains of Bacillus were examined by the fluorescent-antibody technique by using the globulin fraction of serum prepared against spores of B. cereus T. All but one strain of the 59 B. cereus tested fluoresced at the exosporium surface. Fluorescent staining of B. anthracis, B. thuringiensis, and B. mycoides was also observed. Absorption of the globulin fraction with B. anthracis and B. mycoides resulted in the elimination of staining of these organisms. Absorption with B. thuringiensis ATCC 10792 removed antibodies reacting with 6 of the strains of B. thuringiensis tested. Absorption with B. thuringiensis var. galleriae removed antibodies against B. cereus to such a degree that the globulin fraction was unusable.     

Genetic differentiation between sympatric populations of Bacillus cereus and Bacillus thuringiensis

Little is known about genetic exchanges in natural populations of bacteria of the spore-forming Bacillus cereus group, because no population genetics studies have been performed with local sympatric populations. We isolated strains of Bacillus thuringiensis and B. cereus from small samples of soil collected at the same time from two separate geographical sites, one within the forest and the other at the edge of the forest. A total of 100 B. cereus and 98 B. thuringiensis strains were isolated and characterized by electrophoresis to determine allelic composition at nine enzymatic loci. We observed genetic differentiation between populations of B. cereus and B. thuringiensis. Populations of a given Bacillus species--B. thuringiensis or B. cereus--were genetically more similar to each other than to populations of the other Bacillus species. Hemolytic activity provided further evidence of this genetic divergence, which remained evident even if putative clones were removed from the data set. Our results suggest that the rate of gene flow was higher between strains of the same species, but that exchanges between B. cereus and B. thuringiensis were nonetheless possible. Linkage disequilibrium analysis revealed sufficient recombination for B. cereus populations to be considered panmictic units. In B. thuringiensis, the balance between clonal proliferation and recombination seemed to depend on location. Overall, our data indicate that it is not important for risk assessment purposes to determine whether B. cereus and B. thuringiensis belong to a single or two species. Assessment of the biosafety of pest control based on B. thuringiensis requires evaluation of the extent of genetic exchange between strains in realistic natural conditions.     

Genetic relationship in the 'Bacillus cereus group' by rep-PCR fingerprinting and sequencing of a Bacillus anthracis-specific rep-PCR fragment

AIMS: To evaluate the genetic relationship in the Bacillus cereus group by rep-PCR fingerprinting. METHODS AND RESULTS: A collection of 112 strains of the six species of the B. cereus group was analysed by rep-PCR fingerprinting using the BOX-A1R primer. A relative genetic distinctness was found among the species. Cluster analysis of the rep-PCR patterns showed clusters of B. thuringiensis strains quite separate from those of B. cereus strains. The B. anthracis strains represented an independent lineage in a B. cereus cluster. The B. mycoides, B. pseudomycoides and B. weihenstephanensis strains were clustered into three groups at some distance from the other species. Comparison of sequences of AC-390, a typical B. anthracis rep-PCR fragment, from 27 strains of B. anthracis, B. cereus, B. thuringiensis and B. weihenstephanensis, representative of different clusters identified by rep-PCR fingerprinting, confirmed that B. anthracis diverges from its related species. CONCLUSIONS: The genetic relationship deduced from the rep-PCR patterns indicates a relatively clear separation of the six species, suggesting that they can indeed be considered as separate units. SIGNIFICANCE AND IMPACT OF THE STUDY: rep-PCR fingerprinting can make a contribution in the clarification of the genetic relationships between the species of the B. cereus group.     

Bacillus anthracis diverges from related clades of the Bacillus cereus group in 16S-23S ribosomal DNA intergenic transcribed spacers containing tRNA genes

Mung bean nuclease treatment of 16S-23S ribosomal DNA intergenic transcribed spacers (ITS) amplified from several strains of the six species of the Bacillus cereus group showed that B. anthracis Davis TE702 and B. mycoides G2 have other intermediate fragments in addition to the 220- and 550-bp homoduplex fragments typical of the B. cereus group. Long and intermediate homoduplex ITS fragments from strains Davis TE702 and G2 and from another 19 strains of the six species were sequenced. Two main types of ITS were found, either with two tRNA genes (tRNA(Ile) and tRNA(Ala)) or without any at all. Strain Davis TE702 harbors an additional ITS with a single tRNA gene, a hybrid between the tRNA(Ile) and tRNA(Ala) genes, suggesting that a recombination event rather than a deletion generated the single tDNA-containing ITS. Strain G2 showed an additional ITS of intermediate length with no tDNA and no similarity to other known sequences. Neighbor-joining analysis of tDNA-containing long ITS indicated that B. cereus and B. thuringiensis represent a single clade. Three signature sequences discriminated B. anthracis from B. cereus and B. thuringiensis, indicating that the anthrax agent started evolving separately from the related clades of the B. cereus group. B. mycoides and B. weienstephanensis were very closely related, while B. pseudomycoides appeared the most distant species.