Lineage‐specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co‐evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large‐scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co‐evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.     

Genomics of the Bacillus cereus group of organisms

Members of the Bacillus cereus group of organisms include Bacillus cereus, Bacillus anthracis and Bacillus thuringiensis. Collectively, these organisms represent microbes of high economic, medical and biodefense importance. Given this significance, this group contains the highest number of closely related fully sequenced genomes, giving the unique opportunity for thorough comparative genomic analyses. Much of the disease and host specificity of members of this group can be attributed to their plasmids, which vary in size and number. Chromosomes exhibit a high level of synteny and protein similarity, with limited differences in gene content, questioning the speciation of the group members. Genomic data have spurred functional studies that combined microarrays and proteomics. Recent advances are reviewed in this article and highlight the advantages of genomic approaches to the study of this important group of bacteria.     

The autolytic phenotype of the Bacillus cereus group

AIM: To determine the autolytic phenotype of five species in the Bacillus cereus group. METHODS AND RESULTS: The autolytic rate of 96 strains belonging to five species in the B. cereus group was examined under starvation conditions at pH 6, 6.5 and 8.5 in different buffers. The autolytic rate was strain-dependent with a wide variability at pH 6, but higher and more uniform at pH 6.5. At pH 8.5, and respect to the extent of autolysis at pH 6.5, it was relatively low for most of the strains with the lowest values between 13 and 52% in Bacillus mycoides and Bacillus pseudomycoides. Peptidoglycan hydrolase patterns evaluated by renaturing sodium dodecyl sulfate-polyacrylamide gel electrophoresis using cells of Bacillus thuringiensis ssp. tolworthi HD125 as an indicator, revealed complex profiles with lytic bands of about 90, 63, 46, 41, 38, 32, 28 and 25 kDa in B. cereus, B. thuringiensis and Bacillus weihenstephanensis. Bacillus mycoides and B. pseudomycoides had simpler profiles with lytic bands of 63, 46 and 38 kDa. Changes in the autolytic pattern were observed for cells harvested at the stationary phase of growth (72 h) showing an increase in the intensity of the 25 kDa band in the case of B. cereus, B. thuringiensis and B. weihenstephanensis, while no changes were observed for B. mycoides. Using Micrococcus lysodeicticus and Listeria monocytogenes as indicators lytic activity was retained by proteins of 63, 46, 38, 32 and 25 kDa and a new one of about 20 kDa in B. mycoides. Growth in the different media did not affect the autolytic pattern. NaCl abolished the activity of all the peptidoglycan hydrolases except for those of B. mycoides and B. weihenstephanensis. Lytic activity was retained in the presence of MgCl(2), MnCl(2) and EDTA and increased at basic pH. CONCLUSIONS: Bacillus cereus/B. thuringiensis/B. weihenstephanensis showed a high extent of autolysis around neutral pH, even though they presented relatively complex autolysin profiles at alkaline pH. Bacillus mycoides/B. pseudomycoides had a higher extent of autolysis at acidic pH and a simpler autolysin pattern. SIGNIFICANCE AND IMPACT OF THE STUDY: Information on the autolytic phenotype expand the phenotypic characterization of the different species in the B. cereus group.     

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).     

Toxin-producing ability among Bacillus spp. outside the Bacillus cereus group

A total of 333 Bacillus spp. isolated from foods, water, and food plants were examined for the production of possible enterotoxins and emetic toxins using a cytotoxicity assay on Vero cells, the boar spermatozoa motility assay, and a liquid chromatography-mass spectrometry method. Eight strains produced detectable toxins; six strains were cytotoxic, three strains produced putative emetic toxins (different in size from cereulide), and one strain produced both cytotoxin(s) and putative emetic toxin(s). The toxin-producing strains could be assigned to four different species, B. subtilis, B. mojavensis, B. pumilus, or B. fusiformis, by using a polyphasic approach including biochemical, chemotaxonomic, and DNA-based analyses. Four of the strains produced cytotoxins that were concentrated by ammonium sulfate followed by dialysis, and two strains produced cytotoxins that were not concentrated by such a treatment. Two cultures maintained full cytotoxic activity, two cultures reduced their activity, and two cultures lost their activity after boiling. The two most cytotoxic strains (both B. mojavensis) were tested for toxin production at different temperatures. One of these strains produced cytotoxin at growth temperatures ranging from 25 to 42 degrees C, and no reduction in activity was observed even after 24 h of growth at 42 degrees C. The strains that produced putative emetic toxins were tested for the influence of time and temperature on the toxin production. It was shown that they produced putative emetic toxin faster or just as fast at 30 as at 22 degrees C. None of the cytotoxic strains produced B. cereus-like enterotoxins as tested by PCR or by immunological methods.     

Detection of large plasmids from the Bacillus cereus group

The members of the Bacillus cereus group, Bacillus anthracis, Bacillus thuringiensis, and B. cereus senso stricto, are largely defined by their content of large plasmids, which encode major virulence factors. Here we offer an easy, fast, and reliable protocol for the isolation and detection of large plasmids up to the size of at least 350kb. Furthermore, using this method, we report that Bacillus mycoides contain large plasmids.     

蜡状芽胞杆菌群代谢途径的分析和比较

微生物基因组测序和高通量分析方法获得了大量的数据和信息,利用这些信息研究代谢网络成为当前的一个新热点。文章在比较和分析重构代谢网络不同方法的基础上,利用蜡状芽胞杆菌群中已测序的9株蜡状芽胞杆菌、6株炭疽芽胞杆菌、6株苏云金芽胞杆菌基因组,对它们的碳水化合物代谢途径、氨基酸代谢途径和能量代谢途径进行比较与分析,找出它们的共性和特性。这3种菌都存在必需的糖酵解、三羧酸循环、丙氨酸代谢、组氨酸代谢及能量代谢等途径;同时它们还存在特殊的代谢途径,蜡状芽胞杆菌对单糖的利用率较高;炭疽芽胞杆菌的氨基酸降解和转运途径较丰富;苏云金芽胞杆菌中存在催化谷氨酸转化的代谢旁路等。代谢途径的分析为深入研究它们的食物毒素、炭疽毒素和杀虫毒素提供了新思路。     

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.     

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.     

Psychrotolerant species from the Bacillus cereus group are not necessarily Bacillus weihenstephanensis

Twenty-six strains of Bacillus cereus from different sources were determined to be either mesophilic or psychrotrophic by growth at 6 and 42 degrees C. The strains were also screened by two polymerase chain reaction (PCR) methods designed to discriminate between mesophilic and psychrotrophic types. Seventeen of the 26 strains were able to grow at 6 degrees C, but only four conformed to the new psychrotolerant species Bacillus weihenstephanensis. Among the 26 strains were two which caused outbreaks of food poisoning in Norway, and three others that were isolated from food suspected of causing illness. The presence of the gene components encoding production of enterotoxins Nhe, Hbl, EntT and a recently described cytotoxin K was determined by PCR. All the strains possessed genes for at least one of these toxins, and 19 of the 26 strains were cytotoxic in a Vero cell assay. We conclude that there are psychrotrophic B. cereus strains which cannot be classified as B. weihenstephanensis, and that intermediate forms between the two species exist. No correlation between cytotoxicity and the growth temperature of the strains was found.     

Comparative genome analysis of Bacillus cereus group genomes with Bacillus subtilis

Genome features of the Bacillus cereus group genomes (representative strains of Bacillus cereus, Bacillus anthracis and Bacillus thuringiensis sub spp. israelensis) were analyzed and compared with the Bacillus subtilis genome. A core set of 1381 protein families among the four Bacillus genomes, with an additional set of 933 families common to the B. cereus group, was identified. Differences in signal transduction pathways, membrane transporters, cell surface structures, cell wall, and S-layer proteins suggesting differences in their phenotype were identified. The B. cereus group has signal transduction systems including a tyrosine kinase related to two-component system histidine kinases from B. subtilis. A model for regulation of the stress responsive sigma factor sigmaB in the B. cereus group different from the well studied regulation in B. subtilis has been proposed. Despite a high degree of chromosomal synteny among these genomes, significant differences in cell wall and spore coat proteins that contribute to the survival and adaptation in specific hosts has been identified.     

Secondary structure of sphingomyelinase from Bacillus cereus

Of the total of 306 amino acids in the sequence of sphingomyelinase (SMPLC) from Bacillus cereus, almost half (150) are expected to be involved in the formation of loop or turn structure, while 65 and 73 residues may participate in the formation of alpha-helix and beta-structure, respectively. The helix content of SMPLC was calculated to be 0-5%, based on the CD spectra. The addition of divalent metal ions such as Mg2+ or both Ca2+ and Mg2+ had no effect on the CD spectra of SMPLC, although the addition of these metal ions caused the breakdown of membranous SM and specific adsorption of SMPLC onto erythrocyte membranes. A hydropathy study showed that SMPLC has hydrophobic regions at the N-terminal domain which must be responsible for the binding of the enzyme to the membranes. The partial homologies between the amino acid sequences of SMPLC and Clostridium perfringens alpha-toxin (phospholipase C) are discussed.     

Cyclic diguanylate regulation of Bacillus cereus group biofilm formation

Biofilm formation can be considered a bacterial virulence mechanism. In a range of Gram‐negatives, increased levels of the second messenger cyclic diguanylate (c‐di‐GMP) promotes biofilm formation and reduces motility. Other bacterial processes known to be regulated by c‐di‐GMP include cell division, differentiation and virulence. Among Gram‐positive bacteria, where the function of c‐di‐GMP signalling is less well characterized, c‐di‐GMP was reported to regulate swarming motility in Bacillus subtilis while having very limited or no effect on biofilm formation. In contrast, we show that in the Bacillus cereus group c‐di‐GMP signalling is linked to biofilm formation, and to several other phenotypes important to the lifestyle of these bacteria. The Bacillus thuringiensis 407 genome encodes eleven predicted proteins containing domains (GGDEF/EAL) related to c‐di‐GMP synthesis or breakdown, ten of which are conserved through the majority of clades of the B. cereus group, including Bacillus anthracis. Several of the genes were shown to affect biofilm formation, motility, enterotoxin synthesis and/or sporulation. Among these, cdgF appeared to encode a master diguanylate cyclase essential for biofilm formation in an oxygenated environment. Only two cdg genes (cdgA, cdgJ) had orthologs in B. subtilis, highlighting differences in c‐di‐GMP signalling between B. subtilis and B. cereus group bacteria.     

Phosphorylase-mediated mobilization of the amino group of adenine in Bacillus cereus

Mobilization of the ribose moiety of purine nucleosides as well as of the amino group of adenine may be realized in Bacillus cereus by the concerted action of three enzymes: adenosine phosphorylase, adenosine deaminase, and purine nucleoside phosphorylase. In this pathway, ribose-1-phosphate and inorganic phosphate act catalytically, being continuously regenerated by purine nucleoside phosphorylase and adenosine phosphorylase, respectively. As a result of such a metabolic pathway, adenine is quantitatively converted into hypoxanthine, thus overcoming the lack of adenase in B. cereus.     

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.     

Markers for species-specific detection of bacteria from Bacillus cereus group

rpoB and gyr genes (and their fragments) of chromosomal DNA of bacteria from Bacillus cereus group - B. anthracis, B. cereus, and B. thuringiensis - which are the potential markers for their genotyping were sequenced and phylogenetic trees were constructed. Sets of primers for species-specific detection of B. anthracis, B. cereus, and B. thuringiensis by multiplex polymerase chain reaction were designed. Also primers sets, which allow to differentiate strains of B. anthracis with various plasmid profiles (containing both plasmids (pXO1+, pXO2+), and without one (pXO1+, pXO2- or pXO1-, pXO2+) or both plasmids (pXO1-, pXO2-), determining pathogenic characteristics of the strains, were developed. For multiplex PCR primer sets were optimized on the annealing temperature of primers and amplicon length. Itwas shown that phylogenetic tree can be applied as an indicator of reliability and accuracy of taxonomical classification of microorganisms' species and subspecies. Comparison of pXO1 and pXO2 plasmid sequences of B. anthracis showed that these plasmids contain 18 and 4 palindrome sequences respectively which can potentially form thermodynamically stable hairpin-loop structures.