Multilocus sequence typing scheme for bacteria of the Bacillus cereus group

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.

Characterization of Bacillus cereus isolates associated with fatal pneumonias: strains are closely related to Bacillus anthracis and harbor B. anthracis virulence genes

Bacillus cereus is ubiquitous in nature, and while most isolates appear to be harmless, some are associated with food-borne illnesses, periodontal diseases, and other more serious infections. In one such infection, B. cereus G9241 was identified as the causative agent of a severe pneumonia in a Louisiana welder in 1994. This isolate was found to harbor most of the B. anthracis virulence plasmid pXO1 (13). Here we report the characterization of two clinical and one environmental B. cereus isolate collected during an investigation of two fatal pneumonia cases in Texas metal workers. Molecular subtyping revealed that the two cases were not caused by the same strain. However, one of the three isolates was indistinguishable from B. cereus G9241. PCR analysis demonstrated that both clinical isolates contained B. anthracis pXO1 toxin genes. One clinical isolate and the environmental isolate collected from that victim's worksite contained the cap A, B, and C genes required for capsule biosynthesis in B. anthracis. Both clinical isolates expressed a capsule; however, neither was composed of poly-D-glutamic acid. Although most B. cereus isolates are not opportunistic pathogens and only a limited number cause food-borne illnesses, these results demonstrate that some B. cereus strains can cause severe and even fatal infections in patients who appear to be otherwise healthy.

DNA microarray for detection of antibiotic resistance determinants in Bacillus anthracis and closely related Bacillus cereus

We developed a multiplex PCR for amplification of ten genes involved in resistance to ciprofloxacin, doxycycline, rifampin, and vancomycin in Bacillus anthracis and closely related Bacillus cereus. Enzymatic labelling of PCR products followed by hybridization to oligonucleotide probes on a DNA microarray enabled simultaneous detection of resistance genes tetK, tetL, tetM, tetO, vanA, and vanB and resistance-mediating point mutations in genes gyrA, gyrB, parC, and rpoB. The presented assay allows detection of clinically relevant antibiotic resistance determinants within 4h and can be used as a time-saving tool supporting conventional culture-based diagnostics.

Complete sequence analysis of novel plasmids from emetic and periodontal Bacillus cereus isolates reveals a common evolutionary history among the B. cereus-group plasmids, including Bacillus anthracis pXO1

The plasmids of the members of the Bacillus cereus sensu lato group of organisms are essential in defining the phenotypic traits associated with pathogenesis and ecology. For example, Bacillus anthracis contains two plasmids, pXO1 and pXO2, encoding toxin production and encapsulation, respectively, that define this species pathogenic potential, whereas the presence of a Bt toxin-encoding plasmid defines Bacillus thuringiensis isolates. In this study the plasmids from B. cereus isolates that produce emetic toxin or are linked to periodontal disease were sequenced and analyzed. Two periodontal isolates examined contained almost identical approximately 272-kb plasmids, named pPER272. The emetic toxin-producing isolate contained one approximately 270-kb plasmid, named pCER270, encoding the cereulide biosynthesis gene cluster. Comparative sequence analyses of these B. cereus plasmids revealed a high degree of sequence similarity to the B. anthracis pXO1 plasmid, especially in a putative replication region. These plasmids form a newly defined group of pXO1-like plasmids. However, these novel plasmids do not contain the pXO1 pathogenicity island, which in each instance is replaced by plasmid specific DNA. Plasmids pCER270 and pPER272 share regions that are not found in any other pXO1-like plasmids. Evolutionary studies suggest that these plasmids are more closely related to each other than to other identified B. cereus plasmids. Screening of a population of B. cereus group isolates revealed that pXO1-like plasmids are more often found in association with clinical isolates. This study demonstrates that the pXO1-like plasmids may define pathogenic B. cereus isolates in the same way that pXO1 and pXO2 define the B. anthracis species.

Analysis of the life cycle of the soil saprophyte Bacillus cereus in liquid soil extract and in soil

Bacillus is commonly isolated from soils, with organisms of Bacillus cereus sensu lato being prevalent. Knowledge of the ecology of B. cereus and other Bacillus species in soil is far from complete. While the older literature favors a model of growth on soil-associated organic matter, the current paradigm is that B. cereus sensu lato germinates and grows in association with animals or plants, resulting in either symbiotic or pathogenic interactions. An in terra approach to study soil-associated bacteria is described, using filter-sterilized soil-extracted soluble organic matter (SESOM) and artificial soil microcosms (ASM) saturated with SESOM. B. cereus ATCC 14579 displayed a life cycle, with the ability to germinate, grow, and subsequently sporulate in both the liquid SESOM extract and in ASM inserted into wells in agar medium. Cells grew in liquid SESOM without separating, forming multicellular structures that coalesced to form clumps and encasing the ensuing spores in an extracellular matrix. Bacillus was able to translocate from the point of inoculation through soil microcosms as shown by the emergence of outgrowths on the surrounding agar surface. Microscopic inspection revealed bundles of parallel chains inside the soil. The motility inhibitor L-ethionine failed to suppress outgrowth, ruling out translocation by a flagellar-mediated mechanism such as swimming or swarming. Bacillus subtilis subsp. subtilis Marburg and four Bacillus isolates taken at random from soils also displayed a life cycle in SESOM and ASM and were all able to translocate through ASM, even in presence of L-ethionine. These data indicate that B. cereus is a saprophytic bacterium that is able to grow in soil and furthermore that it is adapted to translocate by employing a multicellular mode of growth.

Characterization of Bacillus anthracis-like bacteria isolated from wild great apes from Cote d'Ivoire and Cameroon

We present the microbiological and molecular characterization of bacteria isolated from four chimpanzees and one gorilla thought to have died of an anthrax-like disease in Côte d'Ivoire and Cameroon. These isolates differed significantly from classic Bacillus anthracis by the following criteria: motility, resistance to the gamma phage, and, for isolates from Cameroon, resistance to penicillin G. A capsule was expressed not only after induction by CO(2) and bicarbonate but also under normal growth conditions. Subcultivation resulted in beta-hemolytic activity and gamma phage susceptibility in some subclones, suggesting differences in gene regulation compared to classic B. anthracis. The isolates from Côte d'Ivoire and Cameroon showed slight differences in their biochemical characteristics and MICs of different antibiotics but were identical in all molecular features and sequences analyzed. PCR and Southern blot analyses confirmed the presence of both the toxin and the capsule plasmid, with sizes corresponding to the B. anthracis virulence plasmids pXO1 and pXO2. Protective antigen was expressed and secreted into the culture supernatant. The isolates possessed variants of the Ba813 marker and the SG-749 fragment differing from that of classic B. anthracis strains. Multilocus sequence typing revealed a close relationship of our atypical isolates with both classic B. anthracis strains and two uncommonly virulent Bacillus cereus and Bacillus thuringiensis isolates. We propose that the newly discovered atypical B. anthracis strains share a common ancestor with classic B. anthracis or that they emerged recently by transfer of the B. anthracis plasmids to a strain of the B. cereus group.

Biology of Pseudomonas stutzeri

Pseudomonas stutzeri is a nonfluorescent denitrifying bacterium widely distributed in the environment, and it has also been isolated as an opportunistic pathogen from humans. Over the past 15 years, much progress has been made in elucidating the taxonomy of this diverse taxonomical group, demonstrating the clonality of its populations. The species has received much attention because of its particular metabolic properties: it has been proposed as a model organism for denitrification studies; many strains have natural transformation properties, making it relevant for study of the transfer of genes in the environment; several strains are able to fix dinitrogen; and others participate in the degradation of pollutants or interact with toxic metals. This review considers the history of the discovery, nomenclatural changes, and early studies, together with the relevant biological and ecological properties, of P. stutzeri.

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.

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.

Development of a multilocus sequence typing method for analysis of Lactobacillus plantarum strains

Lactobacillus plantarum is a species of considerable industrial and medical interest. To date, the lack of reliable molecular methods for definite identification at strain level has hindered studies of the population biology of this organism. Here, a multilocus sequence typing (MLST) system for this organism is described, which exploits the genetic variation present in six housekeeping loci to determine the genetic relationship among isolates. The MLST system was established using 16 L. plantarum strains that were also characterized by ribotyping and restriction fragment length polymorphism (RFLP) analysis of the PCR-amplified 16S-23S rDNA intergenic spacer region (ISR). Ribotyping grouped the strains into four groups; however, RFLP analysis of the ISRs showed no differences in the strains analysed. In contrast, MLST had a good discriminatory ability. The sequence analysis of the six genes showed 14 different allelic combinations, with 12 of them represented by only one strain. By using this MLST approach we were able to confirm the identity of two strains deposited in the Spanish Type Culture Collection as different strains. Phylogenetic analysis indicated a panmictic population structure of L. plantarum and split decomposition analysis indicated that recombination plays a role in creating genetic heterogeneity in L. plantarum. As MLST allows precise identification, and easy comparison and exchange of results obtained in different laboratories, the future application of this new molecular method could be useful for the identification of valuable L. plantarum strains.

Evaluation of different methods to discriminate Bacillus anthracis from other bacteria of the Bacillus cereus group

AIMS

To evaluate different methods that are useful for rapid and definitive discrimination of Bacillus anthracis from other bacteria of the Bacillus cereus group in environmental samples like letters claimed to contain anthrax spores.

METHODS AND RESULTS

Characterized strains and bacteria from environmental samples were analysed by microbiological and molecular methods (PCR and restriction analysis). Environmental isolates often shared several microbiological features with B. anthracis, e.g. lack of beta-haemolysis and phospholipase C activity, and only the gamma phage assay was specific for B. anthracis. PCR assays targeting markers from the virulence plasmids exclusively detected B. anthracis, but other PCR targets were also detected in nonanthrax isolates. Additionally, the restriction pattern in an AluI restriction analysis of the SG-749 fragment is not 100% specific. The loci used for multiple-locus variable-number tandem repeat analysis of B. anthracis are also present in other members of the B. cereus group, but amplicon sizes are usually different.

CONCLUSIONS

Environmental samples often contain borderline isolates closely related to B. anthracis both on microbiological and genetic levels. Real-time PCR targeting plasmidal and chromosomal markers should be used for rapid and definitive exclusion of a virulent strain of B. anthracis in such samples.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study gives an overview of the current microbiological and molecular methods used for identification of B. anthracis and shows that most assays have limits when borderline isolates present in environmental samples are analysed.

Multilocus sequence typing of Porphyromonas gingivalis strains from different geographic origins

Porphyromonas gingivalis is an important periodontal pathogen that can be isolated from both active and inactive periodontal lesions. Apparently, differences in virulence between P. gingivalis strains exist, but the mechanisms underlying these differences are not yet fully understood. To obtain more information about pathogenicity and virulence of P. gingivalis, it is relevant to assess the genetic population structure of the species and to examine the occurrence of putative virulence factors against the genetic background. Presently, multilocus sequence typing (MLST) is the best method for analyzing bacterial population structures. Forty P. gingivalis strains from worldwide sources were analyzed by MLST. Internal 310- to 420-bp DNA fragments of the eight ubiquitous chromosomal genes, ftsQ, hagB, gdpxJ, pepO, mcmA, recA, pga, and nah, were amplified by PCR and then sequenced. The number of alleles at individual loci ranged from 2 to 19, and a total of 33 allelic profiles, or sequence types (STs), were identified. Nucleotide variation between alleles was located at one or a few sites. Identical or similar STs were found in isolates from different geographic regions. Our results showed signs of a clonal population structure with a level of recombination not as high as that previously suggested for the species. We also found that P. gingivalis isolates from individual patients were genetically heterogeneous.

Molecular approaches to identify and differentiate Bacillus anthracis from phenotypically similar Bacillus species isolates

Background

Bacillus anthracis and Bacillus cereus can usually be distinguished by standard microbiological methods (e.g., motility, hemolysis, penicillin susceptibility and susceptibility to gamma phage) and PCR. However, we have identified 23 Bacillus spp. isolates that gave discrepant results when assayed by standard microbiological methods and PCR. We used multiple-locus variable-number tandem repeat analysis (MLVA), multiple-locus sequence typing (MLST), and phenotypic analysis to characterize these isolates, determine if they cluster phylogenetically and establish whether standard microbiological identification or PCR were associated with false positive/negative results.

Results

Six isolates were LRN real-time PCR-positive but resistant to gamma phage; MLVA data supported the identification of these isolates as gamma phage-resistant B. anthracis. Seventeen isolates were LRN real-time PCR-negative but susceptible to gamma phage lysis; these isolates appear to be a group of unusual gamma phage-susceptible B. cereus isolates that are closely related to each other and to B. anthracis. All six B. anthracis MLVA chromosomal loci were amplified from one unusual gamma phage-susceptible, motile, B. cereus isolate (although the amplicons were atypical sizes), and when analyzed phylogenetically, clustered with B. anthracis by MLST.

Conclusion

MLVA and MLST aided in the identification of these isolates when standard microbiological methods and PCR could not definitely identify or rule out B. anthracis. This study emphasized the need to perform multiple tests when attempting to identify B. anthracis since relying on a single assay remains problematic due to the diverse nature of bacteria.

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.

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.

Microarray analysis of Bacillus cereus group virulence factors

Bacillus cereus, B. thuringiensis and B. anthracis are closely related medically and economically important bacterial species that belong to the B. cereus group. Members of the B. cereus group carry genes encoding several important virulence factors, including enterotoxins, phospholipases and exotoxins. Since it is difficult to differentiate among B. cereus group members, and because Bacillus virulence factors are very important for pathogenesis, we explored the use of microarray-based detection of virulence factor genes as a tool for strain identification and for determining virulence. Our method requires an initial multiplex PCR amplification step, followed by identification of the PCR amplicons by hybridization to an oligonucleotide microarray containing genes for all three types of Bacillus virulence factors including B. anthracis virulence factors. The DNA chip described here contains 21 identical arrays used for analysis of seven samples in triplicates. Using the arrays, we found that virulence factors are present in several combinations in the strains analyzed. This work also demonstrates the potential of oligonucleotide microarrays for medical, food safety and biodefense analysis of microbial pathogens.

Genome sequence and gene expression of Bacillus anthracis bacteriophage Fah

Fah, a lytic bacteriophage of Bacillus anthracis, is used widely in the former Soviet Union to identify anthrax bacteria. Here, we present the analysis of a 37,974 bp sequence of the Fah genome and examine gene expression of the phage in a model host, Bacillus cereus. Half of the Fah genome contains genes coding for structural proteins and host lysis functions in an arrangement typical of Syphoviridae. The other half of the genome contains genes coding for enzymes of viral genome replication and for numerous predicted transcription factors that are likely to regulate viral gene expression. Primer extension, in vitro transcription assays, and gene array analysis identified temporal classes of Fah genes and allowed location of viral promoters. Fah does not execute host transcription shut-off and relies on host RNA polymerase (RNAP) sigma(A) holoenzyme for transcription of its early and late genes. In addition, Fah encodes a sigma factor, sigma(Fah), a close relative of Bacillus sporulation factor sigma(F) that directs bacterial RNAP to at least one late viral promoter. sigma(Fah) is negatively regulated by host SpoIIAB, an anti-sigma factor that controls sporulation. Thus, sigma(Fah) may link phage gene expression to sporulation of the host.

Cereulide-producing strains of Bacillus cereus show diversity

Producers of cereulide, the emetic toxin of Bacillus cereus, are known to constitute a specific subset within this species. We investigated physiological and genetic properties of 24 strains of B. cereus including two high cereulide producers (600-1,800 ng cereulide mg(-1) wet weight biomass), seven average producers (180-600 ng cereulide mg(-1) wet weight biomass), four low cereulide producers (20-160 ng cereulide mg(-1) wet weight biomass) and 11 non-producers representing isolates from food, food poisoning, human gut and environment. The 13 cereulide producers possessed 16S rRNA gene sequences identical to each other and identical to that of B. anthracis strains Ames, Sterne from GenBank and strain NC 08234-02, but showed diversity in the adk gene (two sequence types), in ribopatterns obtained with EcoRI and PvuII (three types of patterns), in tyrosin decomposition, haemolysis and lecithin hydrolysis (two phenotypes). The cereulide-producing isolates from the human gut represented two ribopatterns of which one was novel to cereulide-producing B. cereus and two phenotypes. We conclude that the cereulide-producing B. cereus are genetically and biochemically more diverse than hitherto thought.

Unusual group II introns in bacteria of the Bacillus cereus group

A combination of sequence and structure analysis and reverse transcriptase PCR experiments was used to characterize the group II introns in the complete genomes of two strains of the pathogen Bacillus cereus. While B. cereus ATCC 14579 harbors a single intron element in the chromosome, B. cereus ATCC 10987 contains three introns in the chromosome and four in its 208-kb pBc10987 plasmid. The most striking finding is the presence in B. cereus ATCC 10987 of an intron [B.c.I2(a)] located on the reverse strand of a gene encoding a putative cell surface protein which appears to be correlated to strains of clinical origin. Because of the opposite orientation of B.c.I2(a), the gene is disrupted. Even more striking is that B.c.I2(a) splices out of an RNA transcript corresponding to the opposite DNA strand. All other intragenic introns studied here are inserted in the same orientation as their host genes and splice out of the mRNA in vivo, setting the flanking exons in frame. Noticeably, B.c.I3 in B. cereus ATCC 10987 represents the first example of a group II intron entirely included within a conserved replication gene, namely, the alpha subunit of DNA polymerase III. Another striking finding is that the observed 3' splice site of B.c.I4 occurs 56 bp after the predicted end of the intron. This apparently unusual splicing mechanism may be related to structural irregularities in the 3' terminus. Finally, we also show that the intergenic introns of B. cereus ATCC 10987 are transcribed with their upstream genes and do splice in vivo.

Discrimination of pathogenic clinical isolates and laboratory strains of Bacillus cereus by NMR-based metabolomic profiling

Six different Bacillus cereus strains were selected from two different ecotypes: (1) three commonly used laboratory strains that are considered avirulent, and (2) three clinical isolates from meningitis patients. Screening of genomic DNA for the presence of genes encoding known toxins gave no candidate genes that were unambiguously able to distinguish between the two groups. However, the application of multivariate pattern-recognition methods to metabolite profiles derived from the different strains using 1H nuclear magnetic resonance spectroscopy (metabolomics) was able to classify the different profiles. The two different ecotypes were clearly separated on the basis of their metabolite profiles, showing that it is possible to use metabolomic methods to classify pathogens on the basis of their expressed physiology, even when it is not possible to infer a direct mechanistic link to specific virulence factors. This metabolomic approach could also have a wide range of possible applications in both general microbiology and microbial ecology for distinguishing and identifying different functional/physiological ecotypes of bacterial strains or species.

Rapid characterization of Bacillus spores targeting species-unique peptides produced with an atmospheric pressure matrix-assisted laser desorption/ionization source

New and improved strategies are eagerly sought for the rapid identification of microorganisms, particularly in mixtures. Mass spectrometry remains a powerful tool for this purpose. Small acid-soluble proteins (SASPs), which are relatively abundant in Bacillus spores, represent potential biomarkers for species characterization. Despite sharing extensive sequence homology, these proteins differ sufficiently in sequence for discrimination between species. This work focuses on the differences in sequence between SASPs from various Bacillus species. Compilation of SASP sequences from protein database searches, followed by in silico trypsin digestion and analysis of the resulting fragments, identified several species-specific peptides that could be targeted for analysis using mass spectrometry. This strategy was tested and found to be successful in the characterization of Bacillus spores both from individual species and in mixtures. Analysis was performed using an ion trap mass spectrometer with an atmospheric pressure MALDI source. This instrumentation offers the advantage of increased speed of analysis and accurate precursor ion selection for tandem mass spectrometric analysis compared with vacuum matrix-assisted laser desorption/ionization and time-of-flight instruments. The identification and targeting of species-specific peptides using this type of instrumentation offers a rapid, efficient strategy for the identification of Bacillus spores and can potentially be applied to different microorganisms.

Multilocus sequence typing reveals thatBacillus cereusstrains 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.

Emetic toxin formation of Bacillus cereus is restricted to a single evolutionary lineage of closely related strains

An in-depth polyphasic approach was applied to study the population structure of the human pathogen Bacillus cereus. To assess the intraspecific biodiversity of this species, which is the causative agent of gastrointestinal diseases, a total of 90 isolates from diverse geographical origin were studied by genetic [M13-PCR, random amplification of polymorphic DNA (RAPD), multilocus sequence typing (MLST)] and phenetic [Fourier transform Infrared (FTIR), protein profiling, biochemical assays] methods. The strain set included clinical strains, isolates from food remnants connected to outbreaks, as well as isolates from diverse food environments with a well documented strain history. The phenotypic and genotypic analysis of the compiled panel of strains illustrated a considerable diversity among B. cereus connected to diarrhoeal syndrome and other non-emetic food strains, but a very low diversity among emetic isolates. Using all typing methods, cluster analysis revealed a single, distinct cluster of emetic B. cereus strains. The isolates belonging to this cluster were neither able to degrade starch nor could they ferment salicin; they did not possess the genes encoding haemolysin BL (Hbl) and showed only weak or no haemolysis. In contrast, haemolytic-enterotoxin-producing B. cereus strains showed a high degree of heterogeneity and were scattered over different clusters when different typing methods were applied. These data provide evidence for a clonal population structure of cereulide-producing emetic B. cereus and indicate that emetic strains represent a highly clonal complex within a potentially panmictic or weakly clonal background population structure of the species. It may have originated only recently through acquisition of specific virulence factors such as the cereulide synthetase gene.

Specificity and polymorphism of the PlcR-PapR quorum-sensing system in the Bacillus cereus group

The expression of extracellular virulence factors in various species of the Bacillus cereus group is controlled by the plcR and papR genes, which encode a transcriptional regulator and a cell-cell signaling peptide, respectively. A processed form of PapR, presumably a pentapeptide, specifically interacts with PlcR to facilitate its binding to its DNA targets. This activating mechanism is strain specific, with this specificity being determined by the first residue of the pentapeptide. We carried out in vivo complementation assays and compared the PlcR-PapR sequences of 29 strains from the B. cereus group. Our findings suggested that the fifth amino acid of the pentapeptide is also involved in the specificity of activation. We identified four classes of PlcR-PapR pairs, defining four distinct pherotypes in the B. cereus group. We used these findings to look at the evolution of the PlcR-PapR quorum-sensing system with regard to the phylogeny of the species forming the B. cereus group.

Bacillus cereus food poisoning and its toxins

The genus Bacillus includes members that demonstrate a wide range of diversity from physiology and ecological niche to DNA sequence and gene regulation. The species of most interest tend to be known for their pathogenicity and are closely linked genetically. Bacillus anthracis causes anthrax, and Bacillus thuringiensis is widely used for its insecticidal properties but has also been associated with foodborne disease. Bacillus cereus causes two types of food poisoning, the emetic and diarrheal syndromes, and a variety of local and systemic infections. Although in this review we provide information on the genus and a variety of species, the primary focus is on the B. cereus strains and toxins that are involved in foodborne illness. B. cereus produces a large number of potential virulence factors, but for the majority of these factors their roles in specific infections have not been established. To date, only cereulide and the tripartite hemolysin BL have been identified specifically as emetic and diarrheal toxins, respectively. Nonhemolytic enterotoxin, a homolog of hemolysin BL, also has been associated with the diarrheal syndrome. Recent findings regarding these and other putative enterotoxins are discussed.

Allelic diversity and population structure in Oenococcus oeni as determined from sequence analysis of housekeeping genes

Oenococcus oeni is the organism of choice for promoting malolactic fermentation in wine. The population biology of O. oeni is poorly understood and remains unclear. For a better understanding of the mode of genetic variation within this species, we investigated by using multilocus sequence typing (MLST) with the gyrB, pgm, ddl, recP, and mleA genes the genetic diversity and genetic relationships among 18 O. oeni strains isolated in various years from wines of the United States, France, Germany, Spain, and Italy. These strains have also been characterized by ribotyping and restriction fragment length polymorphism (RFLP) analysis of the PCR-amplified 16S-23S rRNA gene intergenic spacer region (ISR). Ribotyping grouped the strains into two groups; however, the RFLP analysis of the ISRs showed no differences in the strains analyzed. In contrast, MLST in oenococci had a good discriminatory ability, and we have found a higher genetic diversity than indicated by ribotyping analysis. All sequence types were represented by a single strain, and all the strains could be distinguished from each other because they had unique combinations of alleles. Strains assumed to be identical showed the same sequence type. Phylogenetic analyses indicated a panmictic population structure in O. oeni. Sequences were analyzed for evidence of recombination by split decomposition analysis and analysis of clustered polymorphisms. All results indicated that recombination plays a major role in creating the genetic heterogeneity of O. oeni. A low standardized index of association value indicated that the O. oeni genes analyzed are close to linkage equilibrium. This study constitutes the first step in the development of an MLST method for O. oeni and the first example of the application of MLST to a nonpathogenic food production bacteria.

Microarray-based resequencing of multiple Bacillus anthracis isolates

Custom-designed resequencing arrays were used to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication and by comparison to independently generated shotgun sequence

We used custom-designed resequencing arrays to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication (discrepancy rate of 7.4 × 10^-7) and by comparison to independently generated shotgun sequence (discrepancy rate < 2.5 × 10^-6). Population genomics studies of microbial pathogens using rapid resequencing technologies such as resequencing arrays are critical for recognizing newly emerging or genetically engineered strains.

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.

The bcr1 DNA repeat element is specific to the Bacillus cereus group and exhibits mobile element characteristics

Bacillus cereus strains ATCC 10987 and ATCC 14579 harbor an approximately 155-bp repeated element, bcr1, which is conserved in B. cereus, B. anthracis, B. thuringiensis, and B. mycoides but not in B. subtilis and B. licheniformis. In this study, we show by Southern blot hybridizations that bcr1 is present in all 54 B. cereus group strains tested but absent in 11 Bacillus strains outside the group, suggesting that bcr1 may be specific and ubiquitous to the B. cereus group. By comparative analysis of the complete genome sequences of B. cereus ATCC 10987, B. cereus ATCC 14579, and B. anthracis Ames, we show that bcr1 is exclusively present in the chromosome but absent from large plasmids carried by these strains and that the numbers of full-length bcr1 repeats for these strains are 79, 54, and 12, respectively. Numerous copies of partial bcr1 elements are also present in the three genomes (91, 128, and 53, respectively). Furthermore, the genomic localization of bcr1 is not conserved between strains with respect to chromosomal position or organization of gene neighbors, as only six full-length bcr1 loci are common to at least two of the three strains. However, the intergenic sequence surrounding a specific bcr1 repeat in one of the three strains is generally strongly conserved in the other two, even in loci where bcr1 is found exclusively in one strain. This finding indicates that bcr1 either has evolved by differential deletion from a very high number of repeats in a common ancestor to the B. cereus group or is moving around the chromosome. The identification of bcr1 repeats interrupting genes in B. cereus ATCC 10987 and ATCC 14579 and the presence of a flanking TTTAT motif in each end show that bcr1 exhibits features characteristic of a mobile element.

Impact of genomics on microbial food safety

Genome sequences are now available for many of the microbes that cause food-borne diseases. The information contained in pathogen genome sequences, together with the development of themed and whole-genome DNA microarrays and improved proteomics techniques, might provide tools for the rapid detection and identification of such organisms, for assessing their biological diversity and for understanding their ability to respond to stress. The genomic information also provides insight into the metabolic capacity and versatility of microbes; for example, specific metabolic pathways might contribute to the growth and survival of pathogens in a range of niches, such as food-processing environments and the human host. New concepts are emerging about how pathogens function, both within foods and in interactions with the host. The future should bring the first practical benefits of genome sequencing to the field of microbial food safety, including strategies and tools for the identification and control of emerging pathogens.

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.

Oligonucleotide microarray for identification ofBacillus anthracisbased on intergenic transcribed spacers in ribosomal DNA

We developed a DNA microarray for identification of Bacillus anthracis and other phylogenetic groupings within the "Bacillus cereus group". Nucleotide sequences of 16S-23S ribosomal DNA internal transcribed spacers containing genes for tRNA(Ile) from 52 B. anthracis strains were found to be identical to sequences from seven strains published previously and different from all other bacteria. When 42 oligonucleotide probes targeting polymorphic sites were immobilized on glass slides and hybridized to fluorescently labeled PCR amplification products, one or more mismatches could be discriminated in all but one cases. Hence, hybridization events were highly specific and identification of B. anthracis was straightforward.

Invertebrates as a source of emerging human pathogens

Despite their importance, little is known about the origins of many emerging human pathogens. However, given the age and current predominance of invertebrates, it is likely that bacteria-invertebrate interactions are not only a present source of human pathogens but have also shaped their evolution. Pathogens of invertebrate and unicellular organisms represent an extensive reservoir of bacterial strains equipped with virulence factors that evolved to overcome the innate immune responses of their hosts. This reservoir might represent a source of new human pathogenic strains and might also foster the spread of novel virulence factors into existing human commensal or pathogenic bacteria. This article examines the available evidence for this concept by examining pairs of closely related bacteria, one of which is benign, but insect associated, and one of which is a human pathogen.

Genetic exchange and plasmid transfers in Borrelia burgdorferi sensu stricto revealed by three-way genome comparisons and multilocus sequence typing

Comparative genomics of closely related bacterial isolates is a powerful method for uncovering virulence and other important genome elements. We determined draft sequences (8-fold coverage) of the genomes of strains JD1 and N40 of Borrelia burgdorferi sensu stricto, the causative agent of Lyme disease, and we compared the predicted genes from the two genomes with those from the previously sequenced B31 genome. The three genomes are closely related and are evolutionarily approximately equidistant ( approximately 0.5% pairwise nucleotide differences on the main chromosome). We used a Poisson model of nucleotide substitution to screen for genes with elevated levels of nucleotide polymorphisms. The three-way genome comparison allowed distinction between polymorphisms introduced by mutations and those introduced by recombination using the method of phylogenetic partitioning. Tests for recombination suggested that patches of high-density nucleotide polymorphisms on the chromosome and plasmids arise by DNA exchange. The role of recombination as the main mechanism driving B. burgdorferi diversification was confirmed by multilocus sequence typing of 18 clinical isolates at 18 polymorphic loci. A strong linkage between the multilocus sequence genotypes and the major alleles of outer-surface protein C (ospC) suggested that balancing selection at ospC is a dominant force maintaining B. burgdorferi diversity in local populations. We conclude that B. burgdorferi undergoes genome-wide genetic exchange, including plasmid transfers, and previous reports of its clonality are artifacts from the use of geographically and ecological isolated samples. Frequent recombination implies a potential for rapid adaptive evolution and a possible polygenic basis of B. burgdorferi pathogenicity.

Comparative genetic diversity of Pseudomonas stutzeri genomovars, clonal structure, and phylogeny of the species

A combined phylogenetic and multilocus DNA sequence analysis of 26 Pseudomonas stutzeri strains distributed within the 9 genomovars of the species has been performed. Type strains of the two most closely related species (P. balearica, former genomovar 6, and P. mendocina), together with P. aeruginosa, as the type species of the genus, have been included in the study. The extremely high genetic diversity and the clonal structure of the species were confirmed by the sequence analysis. Clustering of strains in the consensus phylogeny inferred from the analysis of seven nucleotide sequences (16S ribosomal DNA, internally transcribed spacer region 1, gyrB, rpoD, nosZ, catA, and nahH) confirmed the monophyletic origin of the genomovars within the Pseudomonas branch and is in good agreement with earlier DNA-DNA similarity analysis, indicating that the selected genes are representative of the whole genome in members of the species.

The genome sequence of Bacillus cereus ATCC 10987 reveals metabolic adaptations and a large plasmid related to Bacillus anthracis pXO1

We sequenced the complete genome of Bacillus cereus ATCC 10987, a non-lethal dairy isolate in the same genetic subgroup as Bacillus anthracis. Comparison of the chromosomes demonstrated that B.cereus ATCC 10987 was more similar to B.anthracis Ames than B.cereus ATCC 14579, while containing a number of unique metabolic capabilities such as urease and xylose utilization and lacking the ability to utilize nitrate and nitrite. Additionally, genetic mechanisms for variation of capsule carbohydrate and flagella surface structures were identified. Bacillus cereus ATCC 10987 contains a single large plasmid (pBc10987), of approximately 208 kb, that is similar in gene content and organization to B.anthracis pXO1 but is lacking the pathogenicity-associated island containing the anthrax lethal and edema toxin complex genes. The chromosomal similarity of B.cereus ATCC 10987 to B.anthracis Ames, as well as the fact that it contains a large pXO1-like plasmid, may make it a possible model for studying B.anthracis plasmid biology and regulatory cross-talk.