Multilocus Sequence Typing Scheme for Bacteria of the Bacillus cereus Group

In this study we developed a multilocus sequence typing (MLST) scheme for bacteria of the Bacillus cereus group. This group, which includes the species B. cereus, B. thuringiensis, B. weihenstephanensis, and B. anthracis, is known to be genetically very diverse. It is also very important because it comprises pathogenic organisms as well as bacteria with industrial applications. The MLST system was established by using 77 strains having various origins, including humans, animals, food, and soil. A total of 67 of these strains had been analyzed previously by multilocus enzyme electrophoresis, and they were selected to represent the genetic diversity of this group of bacteria. Primers were designed for conserved regions of housekeeping genes, and 330- to 504-bp internal fragments of seven such genes, adk, ccpA, ftsA, glpT, pyrE, recF, and sucC, were sequenced for all strains. The number of alleles at individual loci ranged from 25 to 40, and a total of 53 allelic profiles or sequence types (STs) were distinguished. Analysis of the sequence data showed that the population structure of the B. cereus group is weakly clonal. In particular, all five B. anthracis isolates analyzed had the same ST. The MLST scheme which we developed has a high level of resolution and should be an excellent tool for studying the population structure and epidemiology of the B. cereus group.     

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.     

The Bacillus cereus group: novel aspects of population structure and genome dynamics

AIMS: To provide new insights into the population and genomic structure of the Bacillus cereus group of bacteria. METHODS AND RESULTS: The genetic relatedness among B. cereus group strains was assessed by multilocus sequence typing (MLST) using an optimized scheme based on seven chromosomal housekeeping genes. A set of 48 strains from different clinical sources was included, and six clonal complexes containing several genetically similar isolates from unrelated patients were identified. Interestingly, several clonal groups contained strains that were isolated from similar human sources. Furthermore, comparative whole genome sequence analysis of 16 strains led to the discovery of novel ubiquitous genome features of the B. cereus group, such as atypical group II introns, IStrons, and hitherto uncharacterized repeated elements. CONCLUSIONS: The B. cereus group constitutes a coherent population unified by the presence of ubiquitous and specific genetic elements which do not show any pattern, either in their sequences or genomic locations, which allows to differentiate between the member species of the group. Nevertheless, the population is very dynamic, as particular lineages of clinical origin can evolve to form clonal complexes. At the genome level, the dynamic behaviour is indicated by the presence of numerous mobile and repeated elements. SIGNIFICANCE AND IMPACT OF THE STUDY: The B. cereus group of bacteria comprises species that are of medical and economic importance. The MLST data, along with the primers and protocols used, will be available in a public, web-accessible database (http://mlstoslo.uio.no).     

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.     

Phylogenetic analysis of Bacillus cereus isolates from severe systemic infections using multilocus sequence typing scheme

Bacillus cereus strains from cases of severe or lethal systemic infections, including respiratory symptoms cases, were analyzed using multilocus sequence typing scheme of B. cereus MLST database. The isolates were evenly distributed between the two main clades, and 60% of them had allele profiles new to the database. Half of the collection's strains clustered in a lineage neighboring Bacillus anthracis phylogenetic origin. Strains from lethal cases with respiratory symptoms were allocated in both main clades. This is the first report of strains causing respiratory symptoms to be identified as genetically distant from B. anthracis. The phylogenetic location of the presented here strains was compared with all previously submitted to the database isolates from systemic infections, and were found to appear in the same clusters where clinical isolates from other studies had been assigned. It seems that the pathogenic strains are forming clusters on the phylogenetic tree.     

Rapid genotypic detection of Bacillus anthracis and the Bacillus cereus group by multiplex real-time PCR melting curve analysis

Bacillus anthracis has four plasmid possible virulence genotypes: pXO1+/pXO2+, pXO1+/pXO2-, pXO1-/pXO2+ or pXO1-/pXO2-. Due to the lack of a specific chromosomal marker for B. anthracis, differentiation of the pXO1-/pXO2- form of B. anthracis from closely related Bacillus cereus group species is difficult. In this study, we evaluate the ability of sspE, pXO1 and pXO2 primers to discriminate individual B. anthracis and the B. cereus group genotypes using multiplex real-time PCR and melting curve analysis. Optimal conditions for successful multiplex assays have been established. Purified DNAs from 38 bacterial strains including 11 strains of B. anthracis and 18 B. cereus group strains were analyzed. Nine of the B. cereus group near-neighbor strains were shown by multilocus sequence typing to be phylogenetically proximate to the B. anthracis clade. We have demonstrated that the four plasmid genotypes of B. anthracis and B. cereus group near-neighbors were differentially and simultaneously discriminated by this assay.     

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.     

Multilocus sequence typing reveals that Bacillus cereus strains isolated from clinical infections have distinct phylogenetic origins

Eight strains of Bacillus cereus isolated from bacteremia and soft tissue infections were assigned to seven sequence types (STs) by multilocus sequence typing (MLST). Two strains from different locations had identical STs. The concatenated sequences of the seven STs were aligned with 65 concatenated sequences from reference STs and a neighbor-joining tree was constructed. Two strains were distantly related to all reference STs. Three strains were recovered in a clade that included Bacillus anthracis, B. cereus and rare Bacillus thuringiensis strains while the other three strains were assigned to two STs that were more closely affiliated to most of the B. thuringiensis STs. We conclude that invasive B. cereus strains do not form a single clone or clonal complex of highly virulent strains.     

Allelic diversity and population structure of Bacillus sphaericus as revealed by multilocus sequence typing

The genetic diversity of 35 Bacillus sphaericus strains was analyzed by a newly developed multilocus sequence typing (MLST) scheme, toxin gene pool survey, and mosquito bioassay. The results demonstrated that strains assigned to the same sequence type (ST) had the same occurrence of toxin genes. Further sequence analysis revealed that toxic strains presented a nearly clonal population structure, whereas nontoxic strains had a high level of heterogeneity and were significantly distinct from toxic strains.     

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.     

Comparison of minisatellite polymorphisms in the Bacillus cereus complex: a simple assay for large-scale screening and identification of strains most closely related to Bacillus anthracis

Polymorphism of five tandem repeats that are monomorphic in Bacillus anthracis was investigated in 230 isolates of the B. cereus group and in 5 sequenced B. cereus genomes in search for markers allowing identification of B. cereus and B. thuringiensis strains most closely related to B. anthracis. Using this multiple-locus variable number of tandem repeat analysis (MLVA), a cluster of 30 strains was selected for further characterization. Eventually, six of these were characterized by multilocus sequence type analysis. One of the strains is only six point mutations (of almost 3,000 bp) away from B. anthracis and was also proposed to be closest to B. anthracis by MLVA analysis. However, this strain remains separated from B. anthracis by a number of significant genetic events observed in B. anthracis, including the loss of the hemolysin activity, the presence of four prophages, and the presence of the two virulence plasmids, pXO1 and pXO2. One particular minisatellite marker provides an efficient assay to identify the subset of B. cereus and B. thuringiensis strains closely related to B. anthracis. Based on these results, a very simple assay is proposed that allows the screening of hundreds of strains from the B. cereus complex, with modest equipment and at a low cost, to eventually fill the gap with B. anthracis and better understand the origin and making of this dangerous pathogen.     

Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis--one species on the basis of genetic evidence

Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis are members of the Bacillus cereus group of bacteria, demonstrating widely different phenotypes and pathological effects. B. anthracis causes the acute fatal disease anthrax and is a potential biological weapon due to its high toxicity. B. thuringiensis produces intracellular protein crystals toxic to a wide number of insect larvae and is the most commonly used biological pesticide worldwide. B. cereus is a probably ubiquitous soil bacterium and an opportunistic pathogen that is a common cause of food poisoning. In contrast to the differences in phenotypes, we show by multilocus enzyme electrophoresis and by sequence analysis of nine chromosomal genes that B. anthracis should be considered a lineage of B. cereus. This determination is not only a formal matter of taxonomy but may also have consequences with respect to virulence and the potential of horizontal gene transfer within the B. cereus group.     

Multilocus sequence typing system for the endosymbiont Wolbachia pipientis

The eubacterial genus Wolbachia comprises one of the most abundant groups of obligate intracellular bacteria, and it has a host range that spans the phyla Arthropoda and Nematoda. Here we developed a multilocus sequence typing (MLST) scheme as a universal genotyping tool for Wolbachia. Internal fragments of five ubiquitous genes (gatB, coxA, hcpA, fbpA, and ftsZ) were chosen, and primers that amplified across the major Wolbachia supergroups found in arthropods, as well as other divergent lineages, were designed. A supplemental typing system using the hypervariable regions of the Wolbachia surface protein (WSP) was also developed. Thirty-seven strains belonging to supergroups A, B, D, and F obtained from singly infected hosts were characterized by using MLST and WSP. The number of alleles per MLST locus ranged from 25 to 31, and the average levels of genetic diversity among alleles were 6.5% to 9.2%. A total of 35 unique allelic profiles were found. The results confirmed that there is a high level of recombination in chromosomal genes. MLST was shown to be effective for detecting diversity among strains within a single host species, as well as for identifying closely related strains found in different arthropod hosts. Identical or similar allelic profiles were obtained for strains harbored by different insect species and causing distinct reproductive phenotypes. Strains with similar WSP sequences can have very different MLST allelic profiles and vice versa, indicating the importance of the MLST approach for strain identification. The MLST system provides a universal and unambiguous tool for strain typing, population genetics, and molecular evolutionary studies. The central database for storing and organizing Wolbachia bacterial and host information can be accessed at http://pubmlst.org/wolbachia/.     

Extending the Bacillus cereus group genomics to putative food-borne pathogens of different toxicity

The Bacillus cereus group represents sporulating soil bacteria containing pathogenic strains which may cause diarrheic or emetic food poisoning outbreaks. Multiple locus sequence typing revealed a presence in natural samples of these bacteria of about 30 clonal complexes. Application of genomic methods to this group was however biased due to the major interest for representatives closely related to Bacillus anthracis. Albeit the most important food-borne pathogens were not yet defined, existing data indicate that they are scattered all over the phylogenetic tree. The preliminary analysis of the sequences of three genomes discussed in this paper narrows down the gaps in our knowledge of the B. cereus group. The strain NVH391-98 is a rare but particularly severe food-borne pathogen. Sequencing revealed that the strain should be a representative of a novel bacterial species, for which the name Bacillus cytotoxis or Bacillus cytotoxicus is proposed. This strain has a reduced genome size compared to other B. cereus group strains. Genome analysis revealed absence of sigma B factor and the presence of genes encoding diarrheic Nhe toxin, not detected earlier. The strain B. cereus F837/76 represents a clonal complex close to that of B. anthracis. Including F837/76, three such B. cereus strains had been sequenced. Alignment of genomes suggests that B. anthracis is their common ancestor. Since such strains often emerge from clinical cases, they merit a special attention. The third strain, KBAB4, is a typical facultative psychrophile generally found in soil. Phylogenic studies show that in nature it is the most active group in terms of gene exchange. Genomic sequence revealed high presence of extra-chromosomal genetic material (about 530kb) that may account for this phenomenon. Genes coding Nhe-like toxin were found on a big plasmid in this strain. This may indicate a potential mechanism of toxicity spread from the psychrophile strain community. The results of this genomic work and ecological compartments of different strains incite to consider a necessity of creating prophylactic vaccines against bacteria closely related to NVH391-98 and F837/76. Presumably developing of such vaccines can be based on the properties of non-pathogenic strains such as KBAB4 or ATCC14579 reported here or earlier. By comparing the protein coding genes of strains being sequenced in this project to others we estimate the shared proteome, or core genome, in the B. cereus group to be 3000+/-200 genes and the total proteome, or pan-genome, to be 20-25,000 genes.     

Molecular characterization of Bacillus anthracis using multiplex PCR, ERIC-PCR and RAPD

AIMS: To investigate the molecular characterization of Bacillus anthracis strains by multiplex PCR, enterobacterial repetitive intergenic consensus-PCR (ERIC-PCR) and random amplification of polymorphic DNA (RAPD). METHODS AND RESULTS: Three primers were used to amplify the cya, cap and cereolysinAB genes in the multiplex PCR. Two distinct ERIC-PCR and RAPD fragments, which separated B. anthracis into two groups, were used as probes in Southern hybridization experiments. The probes hybridized only to the cya+ B. anthracis strains identified by the multiplex PCR. Nucleotide sequence analysis of the two cloned fragments showed they were from the pXO1 plasmid of B. anthracis. CONCLUSION: Multiplex PCR simultaneously identified isolates of the Bacillus cereus group and the B. anthracis virulence factors. ERIC-PCR and RAPD, combined with the Southern hybridization analyses, differentiated B. anthracis strains and separated them from the closely related B. cereus group bacteria. SIGNIFICANCE AND IMPACT OF THE STUDY: ERIC-PCR and RAPD assay could be effective in differentiating virulent from avirulent B. anthracis. Our results also show that the amplification of the large plasmids was allowed in the ERIC-PCR and RAPD assay.     

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.     

Identification and characterization of Bacillus anthracis by multiplex PCR analysis of sequences on plasmids pXO1 and pXO2 and chromosomal DNA

Abstract Bacillus anthracis can be identified on the basis of the detection of virulence factor genes located on two plasmids, pXO1 and pXO2. Thus isolates lacking both pXO1 and pXO2 are indistinguishable from closely related B. cereus group bacteria. We developed a multiplex PCR assay for characterization of B. anthracis isolates, and simultaneous confirmation of the species identity independent of plasmid content. The assay amplifies lef, cya, pag (pXO1) and cap (pXO2) genes, and a B. anthracis specific chromosomal marker, giving an easy-to-read profile. This system unambiguously identified virulent (pXO1⁺/2⁺) and avirulent (pXO1⁺/2⁻, pXO1⁻/2⁺ and pXO1⁻/2⁻) strains of B. anthracis and distinguished 'anthrax-like' strains from other B. cereus group bacteria.     

Multiple-locus sequence typing analysis of Bacillus cereus and Bacillus thuringiensis reveals separate clustering and a distinct population structure of psychrotrophic strains

We used multilocus sequence typing (MLST) to characterize phylogenetic relationships for a collection of Bacillus cereus group strains isolated from forest soil in the Paris area during a mild winter. This collection contains multiple strains isolated from the same soil sample and strains isolated from samples from different sites. We characterized 115 strains of this collection and 19 other strains based on the sequences of the clpC, dinB, gdpD, panC, purF, and yhfL loci. The number of alleles ranged from 36 to 53, and a total of 93 allelic profiles or sequence types were distinguished. We identified three major strain clusters-C, T, and W-based on the comparison of individual gene sequences or concatenated sequences. Some less representative clusters and subclusters were also distinguished. Analysis of the MLST data using the concept of clonal complexes led to the identification of two, five, and three such groups in clusters C, T, and W, respectively. Some of the forest isolates were closely related to independently isolated psychrotrophic strains. Systematic testing of the strains of this collection showed that almost all the strains that were able to grow at a low temperature (6 degrees C) belonged to cluster W. Most of these strains, including three independently isolated strains, belong to two clonal complexes and are therefore very closely related genetically. These clonal complexes represent strains corresponding to the previously identified species Bacillus weihenstephanensis. Most of the other strains of our collection, including some from the W cluster, are not psychrotrophic. B. weihenstephanensis (cluster W) strains appear to comprise an effectively sexual population, whereas Bacillus thuringiensis (cluster T) and B. cereus (cluster C) have clonal population structures.     

Fine-Scale Phylogeographic Structure of Borrelia lusitaniae Revealed by Multilocus Sequence Typing

Borrelia lusitaniae is an Old World species of the Lyme borreliosis (LB) group of tick-borne spirochetes and prevails mainly in countries around the Mediterranean Basin. Lizards of the family Lacertidae have been identified as reservoir hosts of B. lusitaniae. These reptiles are highly structured geographically, indicating limited migration. In order to examine whether host geographic structure shapes the evolution and epidemiology of B. lusitaniae, we analyzed the phylogeographic population structure of this tick-borne bacterium using a recently developed multilocus sequence typing (MLST) scheme based on chromosomal housekeeping genes. A total of 2,099 questing nymphal and adult Ixodes ricinus ticks were collected in two climatically different regions of Portugal, being ∼130 km apart. All ticks were screened for spirochetes by direct PCR. Attempts to isolate strains yielded 16 cultures of B. lusitaniae in total. Uncontaminated cultures as well as infected ticks were included in this study. The results using MLST show that the regional B. lusitaniae populations constitute genetically distinct populations. In contrast, no clear phylogeographic signals were detected in sequences of the commonly used molecular markers ospA and ospC. The pronounced population structure of B. lusitaniae over a short geographic distance as captured by MLST of the housekeeping genes suggests that the migration rates of B. lusitaniae are rather low, most likely because the distribution of mediterranean lizard populations is highly parapatric. The study underlines the importance of vertebrate hosts in the geographic spread of tick-borne microparasites.