Identification of a region of genetic variability among Bacillus anthracis strains and related species.

The identification of a region of sequence variability among individual isolates of Bacillus anthracis as well as the two closely related species, Bacillus cereus and Bacillus mycoides, has made a sequence-based approach for the rapid differentiation among members of this group possible. We have identified this region of sequence divergence by comparison of arbitrarily primed (AP)-PCR "fingerprints" generated by an M13 bacteriophage-derived primer and sequencing the respective forms of the only polymorphic fragment observed. The 1,480-bp fragment derived from genomic DNA of the Sterne strain of B. anthracis contained four consecutive repeats of CAATATCAACAA. The same fragment from the Vollum strain was identical except that two of these repeats were deleted. The Ames strain of B. anthracis differed from the Sterne strain by a single-nucleotide deletion. More than 150 nucleotide differences separated B. cereus and B. mycoides from B. anthracis in pairwise comparisons. The nucleotide sequence of the variable fragment from each species contained one complete open reading frame (ORF) (designated vrrA, for variable region with repetitive sequence), encoding a potential 30-kDa protein located between the carboxy terminus of an upstream ORF (designated orf1) and the amino terminus of a downstream ORF (designated lytB). The sequence variation was primarily in vrrA, which was glutamine- and proline-rich (30% of total) and contained repetitive regions. A large proportion of the nucleotide substitutions between species were synonymous. vrrA has 35% identity with the microfilarial sheath protein shp2 of the parasitic worm Litomosoides carinii.     

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.     

Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis

Bacillus anthracis is one of the most genetically homogeneous pathogens described, making strain discrimination particularly difficult. In this paper, we present a novel molecular typing system based on rapidly evolving variable-number tandem repeat (VNTR) loci. Multiple-locus VNTR analysis (MLVA) uses the combined power of multiple alleles at several marker loci. In our system, fluorescently labeled PCR primers are used to produce PCR amplification products from eight VNTR regions in the B. anthracis genome. These are detected and their sizes are determined using an ABI377 automated DNA sequencer. Five of these eight loci were discovered by sequence characterization of molecular markers (vrrC(1), vrrC(2), vrrB(1), vrrB(2), and CG3), two were discovered by searching complete plasmid nucleotide sequences (pXO1-aat and pXO2-at), and one was known previously (vrrA). MLVA characterization of 426 B. anthracis isolates identified 89 distinct genotypes. VNTR markers frequently identified multiple alleles (from two to nine), with Nei's diversity values between 0.3 and 0.8. Unweighted pair-group method arithmetic average cluster analysis identified six genetically distinct groups that appear to be derived from clones. Some of these clones show worldwide distribution, while others are restricted to particular geographic regions. Human commerce doubtlessly has contributed to the dispersal of particular clones in ancient and modern times.     

Molecular evolution and diversity in Bacillus anthracis as detected by amplified fragment length polymorphism markers.

Bacillus anthracis causes anthrax and represents one of the most molecularly monomorphic bacteria known. We have used AFLP (amplified fragment length polymorphism) DNA markers to analyze 78 B. anthracis isolates and six related Bacillus species for molecular variation. AFLP markers are extremely sensitive to even small sequence variation, using PCR and high-resolution electrophoresis to examine restriction fragments. Using this approach, we examined ca. 6.3% of the Bacillus genome for length mutations and ca. 0.36% for point mutations. Extensive variation was observed among taxa, and both cladistic and phenetic analyses were used to construct a phylogeny of B. anthracis and its closest relatives. This genome-wide analysis of 357 AFLP characters (polymorphic fragments) indicates that B. cereus and B. thuringiensis are the closest taxa to B. anthracis, with B. mycoides slightly more distant. B. subtilis, B. polymyxa, and B. stearothermophilus shared few AFLP markers with B. anthracis and were used as outgroups to root the analysis. In contrast to the variation among taxa, only rare AFLP marker variation was observed within B. anthracis, which may be the most genetically uniform bacterial species known. However, AFLP markers did establish the presence or absence of the pXO1 and pXO2 plasmids and detected 31 polymorphic chromosomal regions among the 79 B. anthracis isolates. Cluster analysis identified two very distinct genetic lineages among the B. anthracis isolates. The level of variation and its geographic distribution are consistent with a historically recent African origin for this pathogenic organism. Based on AFLP marker similarity, the ongoing anthrax epidemic in Canada and the northern United States is due to a single strain introduction that has remained stable over at least 30 years and a 1,000-mile distribution.     

Molecular diversity in Bacillus anthracis

Molecular typing of Bacillus anthracis has been extremely difficult due to the lack of polymorphic DNA markers. We have identified nine novel variable number tandemly repeated loci from previously known amplified fragment length polymorphism markers or from the DNA sequence. In combination with the previously known vrrA locus, these markers provide discrimination power to genetically characterize B. anthracis isolates. The variable number tandem repeat (VNTR) loci are found in both gene coding (genic) and non-coding (non-genic) regions. The genic differences are 'in frame' and result in additions or deletion of amino acids to the predicted proteins. Due the rarity of molecular differences, the VNTR changes represent a significant portion of the genetic variation found within B. anthracis. This variation could represent an important adaptive mechanism. Marker similarity and differences among diverse isolates have identified seven major diversity groups that may represent the only world-wide B. anthracis clones. The lineages reconstructed using these data may reflect the dispersal and evolution of this pathogen.     

Molecular characterization of Korean Bacillus anthracis isolates by amplified fragment length polymorphism analysis and multilocus variable-number tandem repeat analysis

We analyzed the genetic relationships and molecular characteristics of 34 Bacillus anthracis isolates from soil and clinical samples in various regions of Korea and 17 related Bacillus species, using the amplified fragment length polymorphism (AFLP) and multilocus variable-number tandem repeat (MLVA) approaches. Triplicate AFLP profiles of these strains showed high reproducibility and identified 376 polymorphisms. AFLP phylogenetic analysis of B. anthracis isolates showed a high level of similarity, 0.93, and this monomorphic fragment profile proved to be useful to differentiate B. anthracis strains from other Bacillus species. The B. cereus group was separated from other Bacillus species at a level of similarity of 0.68. Among them, some B. cereus strains showed genetic interspersion with B. thuringiensis strains. The evolutionary pattern of nucleotide differences among B. anthracis strains with the eight MLVA markers showed nine MLVA types. Three MLVA types, M1 to M3, were pathogenic B. anthracis isolates and were assigned as new genotypes belonging to the A4 and B3 clusters, compared with 89 genotypes deduced from previous data. This indicates that differences in cluster prevalence and distribution may be influenced more by MLVA markers on two plasmids loci and human activity. Consequently, we suggest that the novel MLVA type may represent significant evidence for historic adaptation to environmental conditions of the Asian continent, particularly Korea. Therefore, MLVA techniques may be available for molecular monitoring on anthrax-release-related bioterrorism and further study is required for the continuous epidemiological study of variable anthrax collections.     

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.     

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.     

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.     

Genetic relationships of Bacillus anthracis and closely related species based on variable-number tandem repeat analysis and BOX-PCR genomic fingerprinting

Variable-number tandem repeats (VNTR) analysis and BOX-repeat-based PCR (BOX-PCR) genomic fingerprinting were performed on 25 Bacillus strains to investigate the genetic relatedness of Bacillus anthracis to the closely related species. Based on VNTR analysis, all B. anthracis strains could be assigned to (VNTR)(4), which is the most commonly found type in the world. Interestingly, a (VNTR)(2) was also observed in Bacillus cereus KCTC 1661 and with an exact match to the tandem repeats found in B. anthracis. This finding has never been reported before in the closely related species. According to the BOX-PCR, B. anthracis strains clustered together and separated reliably from the closely related species. However, B. cereus KCTC 1661 was linked to the B. anthracis cluster and showed close relationships with B. anthracis strains. These results indicated that there was a strong correlation between VNTR analysis and BOX-PCR genomic fingerprinting.     

Bacillus anthracis pXO1 Plasmid Sequence Conservation among Closely Related Bacterial Species

The complete sequencing and annotation of the 181.7-kb Bacillus anthracis virulence plasmid pXO1 predicted 143 genes but could only assign putative functions to 45. Hybridization assays, PCR amplification, and DNA sequencing were used to determine whether pXO1 open reading frame (ORF) sequences were present in other bacilli and more distantly related bacterial genera. Eighteen Bacillus species isolates and four other bacterial species were tested for the presence of 106 pXO1 ORFs. Three ORFs were conserved in most of the bacteria tested. Many of the pXO1 ORFs were detected in closely related Bacillus species, and some were detected only in B. anthracis isolates. Three isolates, Bacillus cereus D-17, B. cereus 43881, and Bacillus thuringiensis 33679, contained sequences that were similar to more than one-half of the pXO1 ORF sequences examined. The majority of the DNA fragments that were amplified by PCR from these organisms had DNA sequences between 80 and 98% similar to that of pXO1. Pulsed-field gel electrophoresis revealed large potential plasmids present in both B. cereus 43881 (341 kb) and B. thuringiensis ATCC 33679 (327 kb) that hybridized with a DNA probe composed of six pXO1 ORFs.     

Genetic diversity in the protective antigen gene of Bacillus anthracis

Bacillus anthracis is a gram-positive spore-forming bacterium that causes the disease anthrax. The anthrax toxin contains three components, including the protective antigen (PA), which binds to eucaryotic cell surface receptors and mediates the transport of toxins into the cell. In this study, the entire 2,294-nucleotide protective antigen gene (pag) was sequenced from 26 of the most diverse B. anthracis strains to identify potential variation in the toxin and to further our understanding of B. anthracis evolution. Five point mutations, three synonymous and two missense, were identified. These differences correspond to six different haploid types, which translate into three different amino acid sequences. The two amino acid changes were shown to be located in an area near a highly antigenic region critical to lethal factor binding. Nested primers were used to amplify and sequence this same region of pag from necropsy samples taken from victims of the 1979 Sverdlovsk incident. This investigation uncovered five different alleles among the strains present in the tissues, including two not seen in the 26-sample survey. One of these two alleles included a novel missense mutation, again located just adjacent to the highly antigenic region. Phylogenetic (cladistic) analysis of the pag corresponded with previous strain grouping based on chromosomal variation, suggesting that plasmid evolution in B. anthracis has occurred with little or no horizontal transfer between the different strains.     

Genetic distribution of 295 Bacillus cereus group members based on adk-screening in combination with MLST (Multilocus Sequence Typing) used for validating a primer targeting a chromosomal locus in B. anthracis

The genetic distribution of 295 Bacillus cereus group members has been investigated by using a modified Multilocus Sequence Typing method (MLST). By comparing the nucleic acid sequence of the adk gene fragment, isolates of B. cereus group members most related to B. anthracis may be easily identified. The genetic distribution, with focus on the B. anthracis close neighbours, was used to evaluate a new primer set for specific identification of B. anthracis. This primer set, BA5510-1/2, targeted the putative B. anthracis specific gene BA5510. Real-time PCR using BA5510-1/2 amplified the target fragment from all B. anthracis strains tested and only two (of 289) non-B. anthracis strains analysed. This is one of the most thoroughly validated chromosomal B. anthracis markers for real-time PCR identification, in which the screened collection contained several very closely related B. anthracis strains.     

Differential identification of Bacillus anthracis from environmental Bacillus species using microarray analysis

AIMS: To determine whether microarray analysis could be employed for the differential identification of a range of environmental Bacillus sp. from four strains of Bacillus anthracis. METHODS AND RESULTS: Oligonucleotide probes were designed that were specific to virulence factor genes of B. anthracis (pag, lef and cap), the variable number tandem repeat region of the B. anthracis vrrA gene and to the 16S-23S rRNA intergenic transcribed spacer region (ITS) and pleiotropic regulator (plcR) regions of the Bacillus cereus subgroup species. Generic probes were also designed to hybridize with conserved regions of the 16S rRNA genes of Bacillus (as a positive control), Neisseria sp., Pseudomonas sp., Streptococcus sp., Mycobacterium sp. and to all members of the Enterobacteriaceae to allow simultaneous detection of these bacteria. Identification of B. anthracis was found to rely entirely on hybridization of DNA specific to regions of the pag, lef and cap genes. Cross-reaction was observed between B. anthracis and other Bacillus species with all the other Bacillus probes tested. Results obtained using microarray hybridizations were confirmed using conventional microbiological techniques and found to have very high comparability. CONCLUSIONS: Microarray-based assays are an effective method for the identification of B. anthracis from mixed-culture environmental samples without problems of false-positivity that have been observed with conventional PCR assays. SIGNIFICANCE AND IMPACT OF THE STUDY: Identification of environmental Bacillus sp. by conventional PCR is prone to potential for reporting false-positives. This study provides a method for the exclusion of such isolates.     

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.     

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.     

Sequence and analysis of the DNA encoding protective antigen of Bacillus anthracis

The nucleotide sequence of the protective antigen (PA) gene from Bacillus anthracis and the 5' and 3' flanking sequences were determined. PA is one of three proteins comprising anthrax toxin; and its nucleotide sequence is the first to be reported from B. anthracis. The open reading frame (ORF) is 2319 bp long, of which 2205 bp encode the 735 amino acids of the secreted protein. This region is preceded by 29 codons, which appear to encode a signal peptide having characteristics in common with those of other secreted proteins. A consensus TATAAT sequence was located at the putative -10 promoter site. A Shine-Dalgarno site similar to that found in genes of other Bacillus sp. was located 7 bp upstream from the ATG start codon. The codon usage for the PA gene reflected its high A + T (69%) base composition and differed from those of genes for bacterial proteins from most other sequences examined. The TAA translation stop codon was followed by an inverted repeat forming a potential termination signal. In addition, a 192-codon ORF of unknown significance, theoretically encoding a 21.6-kDa protein, preceded the 5' end of the PA gene.     

Growth characteristics of Bacillus anthracis compared to other Bacillus spp. on the selective nutrient media Anthrax Blood Agar and Cereus Ident Agar

Anthrax Blood Agar (ABA) and Cereus Ident Agar (CEI) were evaluated as selective growth media for the isolation of Bacillus anthracis using 92 B. anthracis and 132 other Bacillus strains from 30 species. The positive predictive values for the identification of B. anthracis on ABA, CEI, and the combination of both were 72%, 71%, and 90%, respectively. Thus, less than 10% of all species were misidentified using both nutrient media. Species which might be misidentified as B. anthracis were B. cereus, B. mycoides, and B. thuringiensis. Particularly, 30% of B. weihenstephanensis strains were misidentified as B. anthracis.