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.     

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.     

Strategy for identification of Bacillus cereus and Bacillus thuringiensis strains closely related to Bacillus anthracis

Bacillus cereus strains that are genetically closely related to B. anthracis can display anthrax-like virulence traits (A. R. Hoffmaster et al., Proc. Natl. Acad. Sci. USA 101:8449-8454, 2004). Hence, approaches that rapidly identify these "near neighbors" are of great interest for the study of B. anthracis virulence mechanisms, as well as to prevent the use of such strains for B. anthracis-based bioweapon development. Here, a strategy is proposed for the identification of near neighbors of B. anthracis based on single nucleotide polymorphisms (SNP) in the 16S-23S rRNA intergenic spacer (ITS) containing tRNA genes, characteristic of B. anthracis. By using restriction site insertion-PCR (RSI-PCR) the presence of two SNP typical of B. anthracis was screened in 126 B. cereus group strains of different origin. Two B. cereus strains and one B. thuringiensis strain showed RSI-PCR profiles identical to that of B. anthracis. The sequencing of the entire ITS containing tRNA genes revealed two of the strains to be identical to B. anthracis. The strict relationship with B. anthracis was confirmed by multilocus sequence typing (MLST) of four other independent loci: cerA, plcR, AC-390, and SG-749. The relationship to B. anthracis of the three strains described by MLST was comparable and even higher to that of four B. cereus strains associated with periodontitis in humans and previously reported as the closest known strains to B. anthracis. SNP in ITS containing tRNA genes combined with RSI-PCR provide a very efficient tool for the identification of strains closely related to B. anthracis.     

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.     

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.     

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 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.     

Fluorescent amplified fragment length polymorphism analysis of Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis isolates

DNA from over 300 Bacillus thuringiensis, Bacillus cereus, and Bacillus anthracis isolates was analyzed by fluorescent amplified fragment length polymorphism (AFLP). B. thuringiensis and B. cereus isolates were from diverse sources and locations, including soil, clinical isolates and food products causing diarrheal and emetic outbreaks, and type strains from the American Type Culture Collection, and over 200 B. thuringiensis isolates representing 36 serovars or subspecies were from the U.S. Department of Agriculture collection. Twenty-four diverse B. anthracis isolates were also included. Phylogenetic analysis of AFLP data revealed extensive diversity within B. thuringiensis and B. cereus compared to the monomorphic nature of B. anthracis. All of the B. anthracis strains were more closely related to each other than to any other Bacillus isolate, while B. cereus and B. thuringiensis strains populated the entire tree. Ten distinct branches were defined, with many branches containing both B. cereus and B. thuringiensis isolates. A single branch contained all the B. anthracis isolates plus an unusual B. thuringiensis isolate that is pathogenic in mice. In contrast, B. thuringiensis subsp. kurstaki (ATCC 33679) and other isolates used to prepare insecticides mapped distal to the B. anthracis isolates. The interspersion of B. cereus and B. thuringiensis isolates within the phylogenetic tree suggests that phenotypic traits used to distinguish between these two species do not reflect the genomic content of the different isolates and that horizontal gene transfer plays an important role in establishing the phenotype of each of these microbes. B. thuringiensis isolates of a particular subspecies tended to cluster together.     

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.     

Differentiation of Bacillus anthracis and other 'Bacillus cereus group' bacteria using IS231-derived sequences

Abstract Sequences based on the conserved 20 bp inverted repeat of IS 231 variants were used as polymerase chain reaction-based fingerprinting primers of the member species of the Bacillus cereus group ( B. anthracis, B. cereus, B. thuringiensis and B. mycoides ), because of their close association with transposons, principally Tn 4430 in B. thuringiensis . Fingerprints of B. anthracis were simple, and specifically allowed its identification and sub-differentiation from other members of the group. Fingerprints for B. cereus were strain-specific; those for B. thuringensis gave a 1650 bp product, characteristic of 1S 231 variants A-F. The same reaction conditions gave one or two bands for both B. anthracis and B. cereus that differed by restriction endonuclease mapping from the B. thuringiensis PCR product and established IS 231 restriction maps; this does not preclude some kind of relationship between these products and IS 231 .     

Use of long-range repetitive element polymorphism-PCR to differentiate Bacillus anthracis strains.

The genome of Bacillus anthracis is extremely monomorphic, and thus individual strains have often proven to be recalcitrant to differentiation at the molecular level. Long-range repetitive element polymorphism-PCR (LR REP-PCR) was used to differentiate various B. anthracis strains. A single PCR primer derived from a repetitive DNA element was able to amplify variable segments of a bacterial genome as large as 10 kb. We were able to characterize five genetically distinct groups by examining 105 B. anthracis strains of diverse geographical origins. All B. anthracis strains produced fingerprints comprising seven to eight bands, referred to as "skeleton" bands, while one to three "diagnostic" bands differentiated between B. anthracis strains. LR REP-PCR fingerprints of B. anthracis strains showed very little in common with those of other closely related species such as B. cereus, B. thuringiensis, and B. mycoides, suggesting relative heterogeneity among the non-B. anthracis strains. Fingerprints from transitional non-B. anthracis strains, which possessed the B. anthracis chromosomal marker Ba813, scarcely resembled those observed for any of the five distinct B. anthracis groups that we have identified. The LR REP-PCR method described in this report provides a simple means of differentiating B. anthracis strains.     

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.     

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.     

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 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.     

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.     

Generation of a specific marker to discriminate Gacillus anthracis from other bacteria of the Bacillus cereus group

Bacillus anthracis is a soil pathogen capable of causing anthrax that is closely related to several environmental species, including B. cereus, B. mycoides, and B. thuringiensis. DNA homology studies showed that B. anthracis, B. cereus, B. mycoides, and B. thuringiensis are closely related, with a high sequence homology. To establish a method to specifically detect B. anthracis in situations such as environmental contamination, we initially performed RAPD-PCR with a 10-mer random primer and confirmed the presence of specific PCR bands only in B. anthracis species. One region specific for B. anthracis was cloned and sequenced, and an internal primer set was designed to amplify a 241-bp DNA fragment within the sequenced region. The PCR system involving these specific primer sets has practical applications. Using lyses methods to prepare the samples for PCR, it was possible to quickly amplify the 241-bp DNA segment from samples containing only a few bacteria. Thus, the PCR detection method developed in this study is expected to facilitate the monitoring of environmental B. anthracis contamination.     

A Bacillus thuringiensis strain producing a polyglutamate capsule resembling that of Bacillus anthracis

Bacillus thuringiensis serovar Monterrey strain BGSC 4AJ1 produced a microscopically visible capsule that reacted with a fluorescent antibody specific for the poly-gamma-d-glutamic acid (PGA) capsule of Bacillus anthracis. PGA capsule biosynthesis genes with 75%, 81%, 72%, 65% and 63% similarity, respectively, to those of the B. anthracis capBCADE cluster were present on a plasmid (pAJ1-1). Strain BGSC 4AJ1, together with five strains of Bacillus cereus that hybridized to a PGA cap gene probe, were analyzed phylogenetically using six housekeeping genes of a B. cereus multilocus sequence typing scheme. Bacillus thuringiensis BGSC 4AJ1 shared four identical alleles with B. anthracis and was the second most closely related to this bacterium of the 674 isolates in the multilocus sequence typing database. The other cap+ strains were distributed among various lineages of Clade 1 of the B. cereus group.     

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.