A phylogenetic analysis of Bacillus thuringiensis serovars by RFLP-based ribotyping

AIMS: To determine the 23S and 5S rRNA gene fingerprints in order to reveal phylogenetic relationships among Bacillus thuringiensis strains. METHODS AND RESULTS: Eighty-six B. thuringiensis strains which include 80 serovar type strains, five intraserovar strains and a non-serotypeable strain, wuhanensis, were tested. Total DNA was digested with EcoRI and HindIII. The 23S and 5S rRNA gene restriction fragment length polymorphisms showed 82 distinctive ribopatterns. The dendrogram generated by numerical analysis showed 10 phylogenetic groups and six ungrouped serovars at the 95.5% DNA relatedness rate. A second dendrogram was constructed using a combination of the data from this study and from a previous study on 16S rRNA gene fingerprinting. It revealed eight distinct phylogenetic groups and three ungrouped serovars at the 94% DNA relatedness rate. CONCLUSION: This method permitted the classification and positioning of a wide variety of B. thuringiensis strains on a phylogenetic tree. Bacillus thuringiensis strains appear to be relatively homogeneous and to share a high degree of DNA relatedness. SIGNIFICANCE AND IMPACT OF THE STUDY: This study contributes a further step to the definition of valid taxonomic sublevels for the B. thuringiensis species.     

Phylogenetic analysis of Bacillus thuringiensis serovars based on 16S rRNA gene restriction fragment length polymorphisms

AIMS: To determine the 16S rRNA gene fingerprints of Bacillus thuringiensis strains to reveal phylogenetic relationships among them. METHODS AND RESULTS: Using 16S rRNA gene restriction fragment length polymorphisms generated by HindIII and EcoRI, 86 Bacillus thuringiensis strains were classified. This includes 80 B. thuringiensis serovars and five more strains, kurstaki HD-1, subtoxicus, dendrolimus, tenebrionis and sandiego, to assess not only interserovar DNA relatedness but also intraserovar DNA relatedness, and the non-motile strain, hence non-serotypeable, B. thuringiensis var. wuhanensis. All 86 B. thuringiensis strains tested showed distinct ribotypes. The dendrogram resulting from the numerical analysis of the distance matrix shows four distinct phylogenetic groups and two ungrouped serovars, finitimus and bolivia, at the 92.5% DNA relatedness rate. CONCLUSION: 16S rRNA gene fingerprinting cannot only be used for the classification of B. thuringiensis strains amenable or not to serotyping, but can also reveal phylogenetic relationships between strains. SIGNIFICANCE AND IMPACT OF THE STUDY: In future screening programmes, 16S rRNA gene restriction pattern analysis could be determined for novel B. thuringiensis strains, allowing them not only to be grouped but also to be positioned on the phylogenetic tree.     

Phylogenetic analysis of Bacillus thuringiensis serovars based on 16S rRNA gene restriction fragment length polymorphisms

K.-B. JOUNG AND J.-C. CÔTÉ. 2001 .
Aims: To determine the 16S rRNA gene fingerprints of Bacillus thuringiensis strains to reveal phylogenetic relationships among them.
Methods and Results: Using 16S rRNA gene restriction fragment length polymorphisms generated by Hin dIII and Eco RI, 86 Bacillus thuringiensis strains were classified. This includes 80 B. thuringiensi s serovars and five more strains, kurstaki HD-1, subtoxicus , dendrolimus , tenebrionis and sandiego , to assess not only interserovar DNA relatedness but also intraserovar DNA relatedness, and the non-motile strain, hence non-serotypeable, B. thuringiensis var. wuhanensis . All 86 B. thuringiensis strains tested showed distinct ribotypes. The dendrogram resulting from the numerical analysis of the distance matrix shows four distinct phylogenetic groups and two ungrouped serovars, finitimus and bolivia , at the 92·5% DNA relatedness rate.
Conclusions: 16S rRNA gene fingerprinting cannot only be used for the classification of B. thuringiensis strains amenable or not to serotyping, but can also reveal phylogenetic relationships between strains.
Significance and Impact of the Study: In future screening programmes, 16S rRNA gene restriction pattern analysis could be determined for novel B. thuringiensis strains, allowing them not only to be grouped but also to be positioned on the phylogenetic tree.     

The occurrence and properties of Bacillus thuringiensis isolated from free-living animals

AIMS: To assess the prevalence and properties of Bacillus thuringiensis isolated from the intestines of small mammals. METHODS AND RESULTS: Bacillus thuringiensis was found in 11% of rodents and 17% of insectivores. Using PFGE of chromosomal DNA, SDS-PAGE of whole-cell proteins and biochemical tests (API system), 12 isolates and three reference strains were classified. Numerical analysis revealed 61% and 89% similarity of protein profiles and biochemical properties of the bacilli, respectively. The results of PFGE were consistent with the outcomes of the analysis of protein profiles. CONCLUSIONS: Although B. thuringiensis is not common in the intestines of small mammals, it is heterogeneous at the genotypic and phenotypic level. SIGNIFICANCE AND IMPACT OF THE STUDY: The results presented here help to explain the diversity and the ecological significance of B. thuringiensis. Future study should focus on the toxic activity of the isolated strains.     

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.     

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.     

Physiological and phylogenetic diversity of thermophilic spore-forming hydrocarbon-oxidizing bacteria from oil fields

The distribution and population density of aerobic hydrocarbon-oxidizing bacteria in the high-temperature oil fields of Western Siberia, Kazakhstan, and China were studied. Seven strains of aerobic thermophilic spore-forming bacteria were isolated from the oil fields and studied by microbiological and molecular biological methods. Based on the 16S rRNA gene sequences, phenotypic characteristics, and the results of DNA-DNA hybridization, the taxonomic affiliation of the isolates was tentatively established. The strains were assigned to the first and fifth subgroups of the genus Bacillus on the phylogenetic branch of the gram-positive bacteria. Strains B and 421 were classified as B. licheniformis. Strains X and U, located between B. stearothermophilus and B. thermocatenulatus on the phylogenetic tree, and strains K, Sam, and 34, related but not identical to B. thermodenitrificans and B. thermoleovorans, undoubtedly represent two new species. Phylogenetically and metabolically related representatives of thermophilic bacilli were found to occur in geographically distant oil fields.     

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 heteroduplex assay for monitoring Bacillus anthracis and close relatives in environmental samples

A fluorescent heteroduplex method was developed to assess the presence of 16S rRNA gene (rDNA) sequences from Bacillus anthracis and close relatives in PCR-amplified 16S rDNA sequence mixtures from environmental samples. The method uses a single-stranded, fluorescent DNA probe, 464 nucleotides in length, derived from a B. anthracis 16S rRNA gene. The probe contains a unique, engineered deletion such that all probe-target duplexes are heteroduplexes with an unpaired G at position 343 (deltaG343). Heteroduplex profiles of sequences >/=85% similar to the probe were produced using an ABI 377 sequencer in less than 3 h. The method divides strains of the Bacillus cereus-Bacillus thuringiensis-B. anthracis group into two subgroups. Each subgroup is defined by a specific 16S rRNA gene sequence type. Sequence type A, containing one mismatch with the probe, occurs in B. anthracis and a small number of closely related clonal lineages represented mostly by food-borne pathogenic isolates of B. cereus and B. thuringiensis. Sequence type B, containing two mismatches with the probe, is found in the majority of B. cereus and B. thuringiensis strains examined to date. Sequence types A and B, when hybridized to the probe, generate two easily differentiated heteroduplexes. Thus, from heteroduplex profiles, the presence of B. cereus-B. thuringiensis-B. anthracis subgroups in environmental samples can be inferred unambiguously. The results show that fluorescent heteroduplex analysis is an effective profiling technique for detection and differentiation of sequences representing small phylogenetic or functional groups in environmental samples.     

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.     

Properties of Bacillus thuringiensis isolated from bank voles

AIMS: To assess the properties of B. thuringiensis naturally occurring in the intestines of bank voles. METHODS AND RESULTS: Seventeen Bacillus thuringiensis strains, exhibiting typical growth on selective medium for the B. cereus group and characterized by the ability to produce parasporal crystals, were isolated from bank voles trapped in the ?omza Landscape Park of the Narew River Valley (north-east Poland). All isolates were characterized by pulsed field gel electrophoresis (PFGE) of chromosomal DNA and SDS polyacrylamide gel electrophoresis (SDS-PAGE) of whole-cell proteins. Six pulsotypes were found with PFGE typing, using SmaI or NotI as restriction enzymes. Significant differences in chromosome size, ranging from 2.4 to 4.2 Mb for the B. thuringiensis strains studied, were noted. Strain heterogeneity in pulsotypes was also reflected by the similarity of whole-cell protein profiles of the strains. Environmental isolates and reference strains grouped at 71% similarity according to SDS-PAGE data and at 84% on the basis of biochemical tests. CONCLUSIONS: B. thuringiensis from intestines of bank voles demonstrated an important level of heterogeneity. The comparison of PFGE profiles and SDS-PAGE of whole-cell protein patterns may be useful to evaluate the relationship between B. thuringiensis isolates. SIGNIFICANCE AND IMPACT OF THE STUDY: The results presented in this paper may help to explain the diversity of B. thuringiensis.     

Molecular and phenotypic characterisation of Bacillus thuringiensis isolated during epizootics in Cydia pomonella L

Twelve Bacillus thuringiensis strains were isolated from intestinal tracts of Cydia pomonella larvae during epizootics in different laboratory insect culture lines. Phenotypic and genetic similarity of these isolates, together with two cultured from Leucoma salicis larvae and 14 reference B. thuringiensis strains were determined. The epizootic bacteria did not form a single group based on numerical analysis of biochemical properties. Simple RAPD method with only one primer does not allow estimating the genetic similarity of B. thuringiensis strains. We propose a novel strategy based on combining several DNA patterns obtained by RAPD technique with different primers for B. thuringiensis typing. Majority of infections in the C. pomonella culture lines were caused by bacteria with different genotypes. However, two isolates cultured from infected insects at different time (one strain was isolated in 1990 and the other in 1992) had identical DNA fingerprint that suggested a possibility of these bacteria to survive in the laboratory and to cause infections in different years. The results of SDS-PAGE of whole-cell proteins revealed a possibility to apply protein profile analysis in epidemiological investigations of infections caused by B. thuringiensis. Strains with identical DNA patterns had very similar whole-cell protein profiles.     

gyrB as a phylogenetic discriminator for members of the Bacillus anthracis-cereus-thuringiensis group

Bacillus anthracis, the causative agent of the human disease anthrax, Bacillus cereus, a food-borne pathogen capable of causing human illness, and Bacillus thuringiensis, a well-characterized insecticidal toxin producer, all cluster together within a very tight clade (B. cereus group) phylogenetically and are indistinguishable from one another via 16S rDNA sequence analysis. As new pathogens are continually emerging, it is imperative to devise a system capable of rapidly and accurately differentiating closely related, yet phenotypically distinct species. Although the gyrB gene has proven useful in discriminating closely related species, its sequence analysis has not yet been validated by DNA:DNA hybridization, the taxonomically accepted "gold standard". We phylogenetically characterized the gyrB sequences of various species and serotypes encompassed in the "B. cereus group," including lab strains and environmental isolates. Results were compared to those obtained from analyses of phenotypic characteristics, 16S rDNA sequence, DNA:DNA hybridization, and virulence factors. The gyrB gene proved more highly differential than 16S, while, at the same time, as analytical as costly and laborious DNA:DNA hybridization techniques in differentiating species within the B. cereus group.     

The phylogenetic diversity of thermophilic members of the genus Bacillus as revealed by 16S rDNA analysis

Abstract Sixteen thermophilic strains of the genus Bacillus , representing eight validly described and six invalidly described species, as well as one unassigned strain, were investigated by comparative 16S rDNA analyses and the sequences compared to the existing database for the genera Bacillus and Alicyclobacillus . The majority of strains were found to cluster in two groups represented by B. stearothermophilus and B. pallidus. Bacillus smithii, B. thermocloacae , and B. thermoruber are phylogenetically well separated and cluster within the radiation of mesophilic bacilli. The as yet undescribed taxon 'B. flavothermus' warrants species status. B. schlegelii and B. tusciae group peripherally with members of Alicyclobacillus and may be reclassified when more phenotypic data support their phylogenetic position.     

Fluorescent Amplified Fragment Length Polymorphism Analysis of Norwegian Bacillus cereus and Bacillus thuringiensis Soil Isolates

We examined 154 Norwegian B. cereus and B. thuringiensis soil isolates (collected from five different locations), 8 B. cereus and 2 B. thuringiensis reference strains, and 2 Bacillus anthracis strains by using fluorescent amplified fragment length polymorphism (AFLP). We employed a novel fragment identification approach based on a hierarchical agglomerative clustering routine that identifies fragments in an automated fashion. No method is free of error, and we identified the major sources so that experiments can be designed to minimize its effect. Phylogenetic analysis of the fluorescent AFLP results reveals five genetic groups in these group 1 bacilli. The ATCC reference strains were restricted to two of the genetic groups, clearly not representative of the diversity in these bacteria. Both B. anthracis strains analyzed were closely related and affiliated with a B. cereus milk isolate (ATCC 4342) and a B. cereus human pathogenic strain (periodontitis). Across the entire study, pathogenic strains, including B. anthracis, were more closely related to one another than to the environmental isolates. Eight strains representing the five distinct phylogenetic clusters were further analyzed by comparison of their 16S rRNA gene sequences to confirm the phylogenetic status of these groups. This analysis was consistent with the AFLP analysis, although of much lower resolution. The innovation of automated genotype analysis by using a replicated and statistical approach to fragment identification will allow very large sample analyses in the future.     

Evaluation of ribosomal RNA gene restriction patterns for the classification of Bacillus species and related genera

AIMS: To identify Bacillus species and related genera by fingerprinting based on ribosomal RNA gene restriction patterns; to compare ribosomal RNA gene restriction patterns-based phylogenetic trees with trees based on 16S rRNA gene sequences; to evaluate the usefulness of ribosomal RNA gene restriction patterns as a taxonomic tool for the classification of Bacillus species and related genera. METHODS AND RESULTS: Seventy-eight bacterial species which include 42 Bacillus species, 31 species from five newly created Bacillus-related genera, and five species from five phenotypically related genera were tested. A total of 77 distinct 16S rRNA gene hybridization banding patterns were obtained. The dendrogram resulting from UPGMA analysis showed three distinct main genetic clusters at the 75% banding pattern similarity. A total of 77 distinct 23S and 5S rRNA genes hybridization banding patterns were obtained, and the dendrogram showed four distinct genetic clusters at the 75% banding pattern similarity. A third dendrogram was constructed using a combination of the data from the 16S rRNA gene fingerprinting and the 23S and 5S rRNA genes fingerprinting. It revealed three distinct main phylogenetic clusters at the 75% banding pattern similarity. CONCLUSIONS: The Bacillus species along with the species from related genera were identified successfully and differentiated by ribosomal RNA gene restriction patterns, and most were distributed with no apparent order in various clusters on each of the three dendrograms. SIGNIFICANCE AND IMPACT OF THE STUDY: Our data indicate that ribosomal RNA gene restriction patterns can be used to reconstruct the phylogeny of the Bacillus species and derived-genera that approximates, but does not duplicate, phylogenies based on 16S rRNA gene sequences.     

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.     

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.