Discrimination of Bacillus cereus and Bacillus thuringiensis with 16S rRNA and gyrB gene based PCR primers and sequencing of their annealing sites

AIMS: To evaluate the possibility for discrimination of Bacillus cereus and B. thuringiensis using 16S rRNA and gyrB gene based PCR methods, and to obtain the sequences of the primer annealing sites so that the PCR results may be explained. METHODS AND RESULTS: Based on the sequence difference in the variable region (V1) of 16S rRNA and in the gyrB gene between B. cereus and B. thuringiensis, PCR primers specific to these Bacillus spp. were designed. When these primers were used to discriminate B. cereus and B. thuringiensis, six of 82 B. cereus strains were identified as B. thuringiensis while 67 of 73 B. thuringiensis strains were identified as B. cereus. Sequence analysis of the primer annealing sites showed that there is no clear-cut difference in the V1 region of 16S rRNA, and in the gyrB gene, between B. cereus and B. thuringiensis strains. CONCLUSIONS: Although 16S rDNA based probes and gyrB gene based PCR primers have been suggested for the discrimination of B. cereus and B. thuringiensis strains, when a large number of Bacillus strains was tested, results showed that discrimination between B. cereus and B. thuringiensis is difficult. Therefore, to distinguish B. thuringiensis from B. cereus, a single feature, such as the presence of a parasporal crystal protein or cry gene, may sometimes be reliable. SIGNIFICANCE AND IMPACT OF THE STUDY: Discrimination between B. cereus and B. thuringiensis is a challenging debate to which this paper makes a contribution.     

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.     

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.     

Markers for species-specific detection of bacteria from Bacillus cereus group

rpoB and gyr genes (and their fragments) of chromosomal DNA of bacteria from Bacillus cereus group - B. anthracis, B. cereus, and B. thuringiensis - which are the potential markers for their genotyping were sequenced and phylogenetic trees were constructed. Sets of primers for species-specific detection of B. anthracis, B. cereus, and B. thuringiensis by multiplex polymerase chain reaction were designed. Also primers sets, which allow to differentiate strains of B. anthracis with various plasmid profiles (containing both plasmids (pXO1+, pXO2+), and without one (pXO1+, pXO2- or pXO1-, pXO2+) or both plasmids (pXO1-, pXO2-), determining pathogenic characteristics of the strains, were developed. For multiplex PCR primer sets were optimized on the annealing temperature of primers and amplicon length. Itwas shown that phylogenetic tree can be applied as an indicator of reliability and accuracy of taxonomical classification of microorganisms' species and subspecies. Comparison of pXO1 and pXO2 plasmid sequences of B. anthracis showed that these plasmids contain 18 and 4 palindrome sequences respectively which can potentially form thermodynamically stable hairpin-loop structures.     

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.     

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.     

Discrimination of Bacillus anthracis and closely related microorganisms by analysis of 16S and 23S rRNA with oligonucleotide microarray

Analysis of 16S rRNA sequences is a commonly used method for the identification and discrimination of microorganisms. However, the high similarity of 16S and 23S rRNA sequences of Bacillus cereus group organisms (up to 99-100%) and repeatedly failed attempts to develop molecular typing systems that would use DNA sequences to discriminate between species within this group have resulted in several suggestions to consider B. cereus and B. thuringiensis, or these two species together with B. anthracis, as one species. Recently, we divided the B. cereus group into seven subgroups, Anthracis, Cereus A and B, Thuringiensis A and B, and Mycoides A and B, based on 16S rRNA, 23S rRNA and gyrB gene sequences and identified subgroup-specific makers in each of these three genes. Here we for the first time demonstrated discrimination of these seven subgroups, including subgroup Anthracis, with a 3D gel element microarray of oligonucleotide probes targeting 16S and 23S rRNA markers. This is the first microarray enabled identification of B. anthracis and discrimination of these seven subgroups in pure cell cultures and in environmental samples using rRNA sequences. The microarray bearing perfect match/mismatch (p/mm) probe pairs was specific enough to discriminate single nucleotide polymorphisms (SNPs) and was able to identify targeted organisms in 5min. We also demonstrated the ability of the microarray to determine subgroup affiliations for B. cereus group isolates without rRNA sequencing. Correlation of these seven subgroups with groupings based on multilocus sequence typing (MLST), fluorescent amplified fragment length polymorphism analysis (AFLP) and multilocus enzyme electrophoresis (MME) analysis of a wide spectrum of different genes, and the demonstration of subgroup-specific differences in toxin profiles, psychrotolerance, and the ability to harbor some plasmids, suggest that these seven subgroups are not based solely on neutral genomic polymorphisms, but instead reflect differences in both the genotypes and phenotypes of the B. cereus group organisms.     

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.     

Siderophores of Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis

Three Bacillus anthracis Sterne strains (USAMRIID, 7702, and 34F2) and Bacillus cereus ATCC 14579 excrete two catecholate siderophores, petrobactin (which contains 3,4-dihydroxybenzoyl moieties) and bacillibactin (which contains 2,3-dihydroxybenzoyl moieties). However, the insecticidal organism Bacillus thuringiensis ATCC 33679 makes only bacillibactin. Analyses of siderophore production by previously isolated [Cendrowski et al., Mol. Microbiol. 52 (2004) 407-417] B. anthracis mutant strains revealed that the B. anthracis bacACEBF operon codes for bacillibactin production and the asbAB gene region is required for petrobactin assembly. The two catecholate moieties also were synthesized by separate routes. PCR amplification identified both asbA and asbB genes in the petrobactin producing strains whereas B. thuringiensis ATCC 33679 retained only asbA. Petrobactin synthesis is not limited to the cluster of B. anthracis strains within the B. cereus sensu lato group (in which B. cereus, B. anthracis, and B. thuringiensis are classified), although petrobactin might be prevalent in strains with pathogenic potential for vertebrates.     

Comparative analysis of 23S ribosomal RNA gene sequences of Bacillus anthracis and emetic Bacillus cereus determined by PCR-direct sequencing

The primary structures of the 23S ribosomal RNA genes of Bacillus anthracis and an emetic strain of Bacillus cereus were determined by direct sequencing of enzymatically amplified chromosomal DNA. The 23S rRNA gene sequences of B. anthracis and B. cereus were found to be almost identical and showed only two differences (a single nucleotide change, and a single base insertion in B. cereus). The feasibility of using PCR-direct sequencing for the rapid sequence determination of large-subunit rRNA genes is demonstrated.     

Genome differences that distinguish Bacillus anthracis from Bacillus cereus and Bacillus thuringiensis

The three species of the group 1 bacilli, Bacillus anthracis, B. cereus, and B. thuringiensis, are genetically very closely related. All inhabit soil habitats but exhibit different phenotypes. B. anthracis is the causative agent of anthrax and is phylogenetically monomorphic, while B. cereus and B. thuringiensis are genetically more diverse. An amplified fragment length polymorphism analysis described here demonstrates genetic diversity among a collection of non-anthrax-causing Bacillus species, some of which show significant similarity to B. anthracis. Suppression subtractive hybridization was then used to characterize the genomic differences that distinguish three of the non-anthrax-causing bacilli from B. anthracis Ames. Ninety-three DNA sequences that were present in B. anthracis but absent from the non-anthrax-causing Bacillus genomes were isolated. Furthermore, 28 of these sequences were not found in a collection of 10 non-anthrax-causing Bacillus species but were present in all members of a representative collection of B. anthracis strains. These sequences map to distinct loci on the B. anthracis genome and can be assayed simultaneously in multiplex PCR assays for rapid and highly specific DNA-based detection of B. anthracis.     

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.     

Distribution of S-layers on the surface of Bacillus cereus strains: phylogenetic origin and ecological pressure

Bacillus anthracis , Bacillus cereus and Bacillus thuringiensis have been described as members of the Bacillus cereus group but are, in fact, one species. B. anthracis is a mammal pathogen, B. thuringiensis an entomopathogen and B. cereus a ubiquitous soil bacterium and an occasional human pathogen. In two clinical isolates of B. cereus , in some B. thuringiensis strains and in B. anthracis , an S-layer has been described. We investigated how the S-layer is distributed in B. cereus , and whether phylogeny or ecology could explain its presence on the surface of some but not all strains. We first developed a simple biochemical assay to test for the presence of the S-layer. We then used the assay with 51 strains of known genetic relationship: 26 genetically diverse B. cereus and 25 non- B. anthracis of the B. anthracis cluster. When present, the genetic organization of the S-layer locus was analysed further. It was identical in B. cereus and B. anthracis . Nineteen strains harboured an S-layer, 16 of which belonged to the B. anthracis cluster. All 19 were B. cereus clinical isolates or B. thuringiensis , except for one soil and one dairy strain. These findings suggest a common phylogenetic origin for the S-layer at the surface of B. cereus strains and, presumably, ecological pressure on its maintenance.     

Identification of self-transmissible plasmids in four Bacillus thuringiensis subspecies.

The transfer of plasmids by mating from four Bacillus thuringiensis subspecies to Bacillus anthracis and Bacillus cereus recipients was monitored by selecting transcipients which acquired plasmid pBC16 (Tcr). Transcipients also inherited a specific large plasmid from each B. thuringiensis donor at a high frequency along with a random array of smaller plasmids. The large plasmids (ca. 50 to 120 megadaltons), pXO13, pXO14, pXO15, and pXO16, originating from B. thuringiensis subsp. morrisoni, B. thuringiensis subsp. toumanoffi, B. thuringiensis subsp. alesti, and B. thuringiensis subsp. israelensis, respectively, were demonstrated to be responsible for plasmid mobilization. Transcipients containing any of the above plasmids had donor capability, while B. thuringiensis strains cured of each of them were not fertile, indicating that the plasmids confer conjugation functions. Confirmation that pXO13, pXO14, and pXO16 were self-transmissible was obtained by the isolation of fertile B. anthracis and B. cereus transcipients that contained only pBC16 and one of these plasmids. pXO14 was efficient in mobilizing the toxin and capsule plasmids, pXO1 and pXO2, respectively, from B. anthracis transcipients to plasmid-cured B. anthracis or B. cereus recipients. DNA-DNA hybridization experiments suggested that DNA homology exists among pXO13, pXO14, and the B. thuringiensis subsp. thuringiensis conjugative plasmids pXO11 and pXO12. Matings performed between strains which each contained the same conjugative plasmid demonstrated reduced efficiency of pBC16 transfer. However, in many instances when donor and recipient strains contained different conjugative plasmids, the efficiency of pBC16 transfer appeared to be enhanced.     

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.     

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 .     

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.