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.     

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.     

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.     

The Ba813 chromosomal DNA sequence effectively traces the whole Bacillus anthracis community

Plasmid genes that are responsible for virulence of Bacillus anthracis are important targets for the DNA-based detection of anthrax. We evaluated the distribution of the Ba813 chromosomal DNA sequence (Ba813) within closely related Bacillus species. Ba813 was systematically identified from 47 strains or isolates of B. anthracis tested, thus indicating its reliability as a tracer for that species. From the 60 strains of closely related Bacillus spp. examined, three bona fide B. cereus and one bona fide B. thuringiensis were found to harbour Ba813. This marker was also detected in Bacillus sp. isolates that were present at high levels in soil samples collected in a place where an anthrax outbreak had occurred. The significance and the possible function of the Ba813 locus is discussed.     

两株炭疽芽胞杆菌MLST 新序列型(ST)

【目的】2株炭疽芽胞杆菌(Bacillus anthracis)17003-14和17003-32的多位点序列分型(Multilocussequence typing,MLST)研究。【方法】选取B.anthracis基因组7个常见管家基因位点glpF、gmk、ilvD、pta、pur、pycA和tpi进行PCR扩增、测序,与MLST数据库中的等位基因序列进行比对,确定菌株的序列型(sequence type,ST)。【结果】B.anthracis 17003-14和17003-32的等位基因编号分别为113、31、1、43、1、53、7和113、31、1、43、1、53、37,比对结果显示这2株细菌的等位基因编号组合未见报道。【结论】17003-14和17003-32为新ST菌株,已被MLST数据库确认,注册号(pubMLST id)分别为id-1053和id-1054。     

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.     

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.     

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.     

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.     

Investigating Burkholderia cepacia complex populations recovered from Italian maize rhizosphere by multilocus sequence typing

The Burkholderia cepacia complex (BCC) comprises at least nine closely related species of abundant environmental microorganisms. Some of these species are highly spread in the rhizosphere of several crop plants, particularly of maize; additionally, as opportunistic pathogens, strains of the BCC are capable of colonizing humans. We have developed and validated a multilocus sequence typing (MLST) scheme for the BCC. Although widely applied to understand the epidemiology of bacterial pathogens, MLST has seen limited application to the population analysis of species residing in the natural environment; we describe its novel application to BCC populations within maize rhizospheres. 115 BCC isolates were recovered from the roots of different maize cultivars from three different Italian regions over a 9-year period (1994-2002). A total of 44 sequence types (STs) were found of which 41 were novel when compared with existing MLST data which encompassed a global database of 1000 clinical and environmental strains representing nearly 400 STs. In this study of rhizosphere isolates approximately 2.5 isolates per ST was found, comparable to that found for the whole BCC population. Multilocus sequence typing also resolved inaccuracies associated with previous identification of the maize isolates based on recA gene restriction fragment length polymorphims and species-specific polymerase chain reaction. The 115 maize isolates comprised the following BCC species groups, B. ambifaria (39%), BCC6 (29%), BCC5 (10%), B. pyrrocinia (8%), B. cenocepacia IIIB (7%) and B. cepacia (6%), with BCC5 and BCC6 potentially constituting novel species groups within the complex. Closely related clonal complexes of strains were identified within B. cepacia, B. cenocepacia IIIB, BCC5 and BCC6, with one of the BCC5 clonal complexes being distributed across all three sampling sites. Overall, our analysis demonstrates that the maize rhizosphere harbours a massive diversity of novel BCC STs, so that their addition to our global MLST database increased the ST diversity by 10%.     

Evaluation of the Biolog system for the identification of Bacillus anthracis

The potential of the Biolog system for the identification of Bacillus anthracis was evaluated. In-house generated databases allowed the correct identification of 19 of 20 isolates of B. anthracis within 24 h. Five strains of the closely related B. cereus/thuringiensis group were misidentified as B. anthracis. For this reason the test could only serve as a primary screen with further testing being required to confirm identity. In addition 20% of all the strains of bacilli examined during the study gave unreadable reaction profiles due to false-positive reactions.     

Multilocus Sequence Typing as a Replacement for Serotyping in Salmonella enterica

Salmonella enterica subspecies enterica is traditionally subdivided into serovars by serological and nutritional characteristics. We used Multilocus Sequence Typing (MLST) to assign 4,257 isolates from 554 serovars to 1092 sequence types (STs). The majority of the isolates and many STs were grouped into 138 genetically closely related clusters called eBurstGroups (eBGs). Many eBGs correspond to a serovar, for example most Typhimurium are in eBG1 and most Enteritidis are in eBG4, but many eBGs contained more than one serovar. Furthermore, most serovars were polyphyletic and are distributed across multiple unrelated eBGs. Thus, serovar designations confounded genetically unrelated isolates and failed to recognize natural evolutionary groupings. An inability of serotyping to correctly group isolates was most apparent for Paratyphi B and its variant Java. Most Paratyphi B were included within a sub-cluster of STs belonging to eBG5, which also encompasses a separate sub-cluster of Java STs. However, diphasic Java variants were also found in two other eBGs and monophasic Java variants were in four other eBGs or STs, one of which is in subspecies salamae and a second of which includes isolates assigned to Enteritidis, Dublin and monophasic Paratyphi B. Similarly, Choleraesuis was found in eBG6 and is closely related to Paratyphi C, which is in eBG20. However, Choleraesuis var. Decatur consists of isolates from seven other, unrelated eBGs or STs. The serological assignment of these Decatur isolates to Choleraesuis likely reflects lateral gene transfer of flagellar genes between unrelated bacteria plus purifying selection. By confounding multiple evolutionary groups, serotyping can be misleading about the disease potential of S. enterica. Unlike serotyping, MLST recognizes evolutionary groupings and we recommend that Salmonella classification by serotyping should be replaced by MLST or its equivalents.     

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.     

Genotyping of Bacillus cereus strains by microarray-based resequencing

The ability to distinguish microbial pathogens from closely related but nonpathogenic strains is key to understanding the population biology of these organisms. In this regard, Bacillus anthracis, the bacterium that causes inhalational anthrax, is of interest because it is closely related and often difficult to distinguish from other members of the B. cereus group that can cause diverse diseases. We employed custom-designed resequencing arrays (RAs) based on the genome sequence of Bacillus anthracis to generate 422 kb of genomic sequence from a panel of 41 Bacillus cereus sensu lato strains. Here we show that RAs represent a "one reaction" genotyping technology with the ability to discriminate between highly similar B. anthracis isolates and more divergent strains of the B. cereus s.l. Clade 1. Our data show that RAs can be an efficient genotyping technology for pre-screening the genetic diversity of large strain collections to selected the best candidates for whole genome sequencing.     

Using multilocus sequence data to define the pneumococcus

We investigated the genetic relationships between serotypeable pneumococci and nonserotypeable presumptive pneumococci using multilocus sequence typing (MLST) and partial sequencing of the pneumolysin gene (ply). Among 121 nonserotypeable presumptive pneumococci from Finland, we identified isolates of three classes: those with sequence types (STs) identical to those of serotypeable pneumococci, suggesting authentic pneumococci in which capsular expression had been downregulated or lost; isolates that clustered among serotypeable pneumococci on a tree based on the concatenated sequences of the MLST loci but which had STs that differed from those of serotypeable pneumococci in the MLST database; and a more diverse collection of isolates that did not cluster with serotypeable pneumococci. The latter isolates typically had sequences at all seven MLST loci that were 5 to 10% divergent from those of authentic pneumococci and also had distinct and divergent ply alleles. These isolates are proposed to be distinct from pneumococci but cannot be resolved from them by optochin susceptibility, bile solubility, or the presence of the ply gene. Complete resolution of pneumococci from the related but distinct population is problematic, as recombination between them was evident, and a few isolates of each population possessed alleles at one or occasionally more MLST loci from the other population. However, a tree based on the concatenated sequences of the MLST loci in most cases unambiguously distinguished whether a nonserotypeable isolate was or was not a pneumococcus, and the sequence of the ply gene fragment was found to be useful to resolve difficult cases.     

Multilocus Sequence Typing of Uropathogenic ESBL-Producing <Emphasis Type="Italic">Escherichia coli</Emphasis> Isolated in a Brazilian Community

In the present study, multilocus sequence typing (MLST), pulsed-field gel electrophoresis (PFGE) and polymerase chain reaction detection of three resistance genes were combined to characterize seven uropathogenic E. coli isolated from outpatients. Selected portions of seven housekeeping and three antibiotic-resistance genes of the isolates were sequenced. The seven isolates were classified into four different sequence types (STs) by MLST and five PGFE types. Three isolates had a novel allelic profile representing a new ST designated as ST528 and showed the same PFGE and resistance genes. Two isolates, both characterized as ST359, were differentiated by PFGE and shared only one of the antibiotic-resistance genes studied. Comparison of MLST results with those of PFGE and resistance genes demonstrated that Escherichia coli had acquired different antibiotic-resistance genes and DNA rearrangements, causing alterations in PFGE patterns but maintaining the same ST. Furthermore, this article also reports the first detection of a CTX-M-2 ESBL E. coli and SHV-5 in a Brazilian community.     

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.     

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 β -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 Alu I 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.     

Multi-Locus Sequence Typing of a Geographically and Temporally Diverse Sample of the Highly Clonal Human Pathogen Bartonella quintana

Bartonella quintana is a re-emerging pathogen and the causative agent of a variety of disease manifestations in humans including trench fever. Various typing methods have been developed for B. quintana, but these tend to be limited by poor resolution and, in the case of gel-based methods, a lack of portability. Multilocus sequence typing (MLST) has been used to study the molecular epidemiology of a large number of pathogens, including B. henselae, a close relative of B. quintana. We developed a MLST scheme for B. quintana based on the 7 MLST loci employed for B. henselae with two additional loci to cover underrepresented regions of the B. quintana chromosome. A total of 16 B. quintana isolates spanning over 60 years and three continents were characterized. Allelic variation was detected in five of the nine loci. Although only 8/4270 (0.002%) of the nucleotide sites examined were variable over all loci, these polymorphisms resolved the 16 isolates into seven sequence types (STs). We also demonstrate that MLST can be applied on uncultured isolates by direct PCR from cardiac valve tissue, and suggest this method presents a promising approach for epidemiological studies in this highly clonal organism. Phylogenetic and clustering analyses suggest that two of the seven STs form a distinct lineage within the population.     

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.