The structure of the major cell wall polysaccharide of Bacillus anthracis is species-specific

In this report we describe the structure of the polysaccharide released from Bacillus anthracis vegetative cell walls by aqueous hydrogen fluoride (HF). This HF-released polysaccharide (HF-PS) was isolated and structurally characterized from the Ames, Sterne, and Pasteur strains of B. anthracis. The HF-PSs were also isolated from the closely related Bacillus cereus ATCC 10987 strain, and from the B. cereus ATCC 14579 type strain and compared with those of B. anthracis. The structure of the B. anthracis HF-PS was determined by glycosyl composition and linkage analyses, matrix-assisted laser desorption-time of flight mass spectrometry, and one- and two-dimensional nuclear magnetic resonance spectroscopy. The HF-PSs from all of the B. anthracis isolates had an identical structure consisting of an amino sugar backbone of -->6)-alpha-GlcNAc-(1-->4)-beta-ManNAc-(1-->4)-beta-GlcNAc-(1-->, in which the alpha-GlcNAc residue is substituted with alpha-Gal and beta-Gal at O-3 and O-4, respectively, and the beta-GlcNAc substituted with alpha-Gal at O-3. There is some variability in the presence of two of these three Gal substitutions. Comparison with the HF-PSs from B. cereus ATCC 10987 and B. cereus ATCC 14579 showed that the B. anthracis structure was clearly different from each of these HF-PSs and, furthermore, that the B. cereus ATCC 10987 HF-PS structure was different from that of B. cereus ATCC 14579. The presence of a B. anthracis-specific polysaccharide structure in its vegetative cell wall is discussed with regard to its relationship to those of other Bacillus species.     

Secondary cell wall polysaccharides from Bacillus cereus strains G9241, 03BB87 and 03BB102 causing fatal pneumonia share similar glycosyl structures with the polysaccharides from Bacillus anthracis

Secondary cell wall polysaccharides (SCWPs) are important structural components of the Bacillus cell wall and contribute to the array of antigens presented by these organisms in both spore and vegetative forms. We previously found that antisera raised to Bacillus anthracis spore preparations cross-reacted with SCWPs isolated from several strains of pathogenic B. cereus, but did not react with other phylogenetically related but nonpathogenic Bacilli, suggesting that the SCWP from B. anthracis and pathogenic B. cereus strains share specific structural features. In this study, SCWPs from three strains of B. cereus causing severe or fatal pneumonia (G9241, 03BB87 and 03BB102) were isolated and subjected to structural analysis and their structures were compared to SCWPs from B. anthracis. Complete structural analysis was performed for the B. cereus G9241 SCWP using NMR spectroscopy, mass spectrometry and derivatization methods. The analyses show that SCWPs from B. cereus G9241 has a glycosyl backbone identical to that of B. anthracis SCWP, consisting of multiple trisaccharide repeats of: →6)-α-d-GlcpNAc-(1?→?4)-β-d-ManpNAc-(1?→?4)-β-d-GlcpNAc-(1→. Both the B. anthracis and pathogenic B. cereus SCWPs are highly substituted at all GlcNAc residues with α- and β-Gal residues, however, only the SCWPs from B. cereus G9241 and 03BB87 carry an additional α-Gal substitution at O-3 of ManNAc residues, a feature lacking in the B. anthracis SCWPs. Both the B. anthracis and B. cereus SCWPs are pyruvylated, with an approximate molecular mass of ≈12,000?Da. The implications of these findings regarding pathogenicity and cell wall structure are discussed.     

Environmental and biofilm-dependent changes in a Bacillus cereus secondary cell wall polysaccharide

Bacterial species from the Bacillus genus, including Bacillus cereus and Bacillus anthracis, synthesize secondary cell wall polymers (SCWP) covalently associated to the peptidoglycan through a phospho-diester linkage. Although such components were observed in a wide panel of B. cereus and B. anthracis strains, the effect of culture conditions or of bacterial growth state on their synthesis has never been addressed. Herein we show that B. cereus ATCC 14579 can synthesize not only one, as previously reported, but two structurally unrelated secondary cell wall polymers (SCWP) polysaccharides. The first of these SCWP, →4)[GlcNAc(β1-3)]GlcNAc(β1-6)[Glc(β1-3)][ManNAc(α1-4)]GalNAc(α1-4)ManNAc(β1→, although presenting an original sequence, fits to the already described the canonical sequence motif of SCWP. In contrast, the second polysaccharide was made up by a totally original sequence, →6)Gal(α1-2)(2-R-hydroxyglutar-5-ylamido)Fuc2NAc4N(α1-6)GlcNAc(β1→, which no equivalent has ever been identified in the Bacillus genus. In addition, we established that the syntheses of these two polysaccharides were differently regulated. The first one is constantly expressed at the surface of the bacteria, whereas the expression of the second is tightly regulated by culture conditions and growth states, planktonic, or biofilm.     

Peptidoglycan N-acetylglucosamine deacetylases from Bacillus cereus, highly conserved proteins in Bacillus anthracis

The genomes of Bacillus cereus and its closest relative Bacillus anthracis contain 10 polysaccharide deacetylase homologues. Six of these homologues have been proposed to be peptidoglycan N-acetylglucosamine deacetylases. Two of these genes, namely bc1960 and bc3618, have been cloned and expressed in Escherichia coli, and the recombinant enzymes have been purified to homogeneity and further characterized. Both enzymes were effective in deacetylating cell wall peptidoglycan from the Gram(+) Bacillus cereus and Bacillus subtilis and the Gram(-) Helicobacter pylori as well as soluble chitin substrates and N-acetylchitooligomers. However, the enzymes were not active on acetylated xylan. These results provide insight into the substrate specificity of carbohydrate esterase family 4 enzymes. It was revealed that both enzymes deacetylated only the GlcNAc residue of the synthetic muropeptide N-acetyl-D-glucosamine-(beta-1,4)-N-acetylmuramyl-L-alanine-D-isoglutamine. Analysis of the constituent muropeptides of peptidoglycan from B. subtilis and H. pylori resulting from incubation of the enzymes BC1960 and BC3618 with these polymers and subsequent hydrolysis by Cellosyl and mutanolysin, respectively, similarly revealed that both enzymes deacetylate GlcNAc residues of peptidoglycan. Kinetic analysis toward GlcNAc(2-6) revealed that GlcNAc4 was the favorable substrate for both enzymes. Identification of the sequence of N-acetychitooligosaccharides (GlcNAc(2-4)) following enzymatic deacetylation by using 1H NMR revealed that both enzymes deacetylate all GlcNAc residues of the oligomers except the reducing end ones. Enzymatic deacetylation of chemically acetylated vegetative peptidoglycan from B. cereus by BC1960 and BC3618 resulted in increased resistance to lysozyme digestion. This is the first biochemical study of bacterial peptidoglycan N-acetylglucosamine deacetylases.     

Metabolic capacity of Bacillus cereus strains ATCC 14579 and ATCC 10987 interlinked with comparative genomics

Bacillus cereus is an important food-borne pathogen and spoilage organism. In this study, numerous phenotypes and the genomes of B. cereus strains ATCC 14579 and ATCC 10987 were analysed to compare their metabolic capacity and stress resistance potential. The growth performance of the two strains was assessed for nearly 2000 phenotypes, including use of nutrient sources, performance in acid and basic environments, osmo-tolerance and antibiotic resistance. Several food-relevant phenotypic differences were found between ATCC 14579 and ATCC 10987, such as differences in utilization of carbohydrates, peptides, amino acids and ammonia. Subsequently, the genomes of both strains were analysed with INPARANOID to search for strain-specific open reading frames (ORFs). B. cereus ATCC 14579 and ATCC 10987 were found to harbour 983 and 1360 strain-specific ORFs respectively. The strain-specific phenotypic features were interlinked with corresponding genetic features and for several phenotypic differences a related strain-specific genetic feature could be identified. In conclusion, the combination of phenotypic data with strain-specific genomic differences has led to detailed insight into the performance of the two B. cereus strains, and may supply indicators for the performance of these bacteria in different environments and ecological niches.     

Characterization of type II and III restriction-modification systems from Bacillus cereus strains ATCC 10987 and ATCC 14579

The genomes of two Bacillus cereus strains (ATCC 10987 and ATCC 14579) have been sequenced. Here, we report the specificities of type II/III restriction (R) and modification (M) enzymes. Found in the ATCC 10987 strain, BceSI is a restriction endonuclease (REase) with the recognition and cut site CGAAG 24-25/27-28. BceSII is an isoschizomer of AvaII (G/GWCC). BceSIII cleaves at ACGGC 12/14. The BceSIII C terminus resembles the catalytic domains of AlwI, MlyI, and Nt.BstNBI. BceSIV is composed of two subunits and cleaves on both sides of GCWGC. BceSIV activity is strongly stimulated by the addition of cofactor ATP or GTP. The large subunit (R1) of BceSIV contains conserved motifs of NTPases and DNA helicases. The R1 subunit has no endonuclease activity by itself; it strongly stimulates REase activity when in complex with the R2 subunit. BceSIV was demonstrated to hydrolyze GTP and ATP in vitro. BceSIV is similar to CglI (GCSGC), and homologs of R1 are found in 11 sequenced bacterial genomes, where they are paired with specificity subunits. In addition, homologs of the BceSIV R1-R2 fusion are found in many sequenced microbial genomes. An orphan methylase, M.BceSV, was found to modify GCNGC, GGCC, CCGG, GGNNCC, and GCGC sites. A ParB-methylase fusion protein appears to nick DNA nonspecifically. The ATCC 14579 genome encodes an active enzyme Bce14579I (GCWGC). BceSIV and Bce14579I belong to the phospholipase D (PLD) family of endonucleases that are widely distributed among Bacteria and Archaea. A survey of type II and III restriction-modification (R-M) system genes is presented from sequenced B. cereus, Bacillus anthracis, and Bacillus thuringiensis strains.     

Physical maps of the genomes of three Bacillus cereus strains.

NotI restriction maps of the chromosomes from Bacillus cereus ATCC 10876, ATCC 11778, and the B. cereus type strain ATCC 14579 have been established and compared with the previously established map of B. cereus ATCC 10987. Between 10 and 14 NotI fragments were observed, ranging from 15 to 1,300 kb, in digests of DNA from the various strains. The sizes of the genomes varied between 5.4 and 6.3 Mb. The maps were constructed by hybridization of 42 random probes, prepared from B. cereus ATCC 10987 libraries, to fragments from partial and complete NotI digests, separated by pulsed-field gel electrophoresis. Nine probes were specific for ATCC 10987 only. Probes for five B. subtilis and five B. cereus genes were also used. The NotI restriction fragment patterns of the four strains were strikingly different.     

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.     

Siderophores of Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis

Three Bacillus anthracis Sterne strains (USAMRIID, 7702, and 34F2) and Bacillus cereus ATCC 14579 excrete two catecholate siderophores, petrobactin (which contains 3,4-dihydroxybenzoyl moieties) and bacillibactin (which contains 2,3-dihydroxybenzoyl moieties). However, the insecticidal organism Bacillus thuringiensis ATCC 33679 makes only bacillibactin. Analyses of siderophore production by previously isolated [Cendrowski et al., Mol. Microbiol. 52 (2004) 407-417] B. anthracis mutant strains revealed that the B. anthracis bacACEBF operon codes for bacillibactin production and the asbAB gene region is required for petrobactin assembly. The two catecholate moieties also were synthesized by separate routes. PCR amplification identified both asbA and asbB genes in the petrobactin producing strains whereas B. thuringiensis ATCC 33679 retained only asbA. Petrobactin synthesis is not limited to the cluster of B. anthracis strains within the B. cereus sensu lato group (in which B. cereus, B. anthracis, and B. thuringiensis are classified), although petrobactin might be prevalent in strains with pathogenic potential for vertebrates.     

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

Variable numbers of rRNA gene operons in Bacillus cereus strains

Abstract Ribosornal RNA operon organisation was analysed in two Bacillus cereus strains of different chromosome size, ATCC 10987 (5.4 Mb) and F0837/76 (2.4 Mb). We estimated that there were twelve and nine copies of the rRNA operons in these two strains, respectively. In B. cereus ATCC 10987 six rRNA operons were less than 10 kb apart, while in B. cereus F0837/76 four rRNA operons were similarly clustered. The origin of replication was located in the vicinity of a rRNA operon in both strains.     

Proposed pathway for biosynthesis of the S-layer glycoprotein of Bacillus alvei.

The outer surface of the murein sacculus of the eubacterium Bacillus alvei is covered by a surface layer (S-layer) glycoprotein. The glycan chain of this S-layer glycoprotein consists of trisaccharide repeating units with ManNAc, Gal, and Glc as constituents. From cell extracts of B. alvei, nucleotide-activated derivatives of ManNAc, Gal, Glc, and GlcNAc were isolated. Furthermore, GDP- and dolichyl-activated oligosaccharides were obtained. On the basis of the isolated putative glycoprotein precursors, a pathway for the biosynthesis of the oligosaccharide chain is proposed.     

Characterization of the chemotaxis fli Y and che A genes in Bacillus cereus

This paper describes the first identification of chemotaxis genes in Bacillus cereus. We sequenced and studied the genomic organization and the expression of the cheA and fliY genes in two different B. cereus strains, ATCC 14579 and ATCC 10987. While cheA encodes a highly conserved protein acting as the main regulator of the chemotactic response in flagellated eubacteria, fliY, which has been previously described only in B. subtilis, is one of the three genes encoding proteins of the flagellar switch complex. Although the sequences and relative position of cheA and fliY were found to be identical in the two B. cereus strains analyzed, the restriction fragment containing both genes was located differently on the physical maps of B. cereus ATCC 14579 and ATCC 10987. Evidence is shown that the genomic organization and the expression of fliY and cheA in B. cereus differ significantly from that described for B. subtilis, which is considered a model microorganism for chemotaxis in gram-positive bacteria.     

Genes of Bacillus cereus and Bacillus anthracis encoding proteins of the exosporium

The exosporium is the outermost layer of spores of Bacillus cereus and its close relatives Bacillus anthracis and Bacillus thuringiensis. For these pathogens, it represents the surface layer that makes initial contact with the host. To date, only the BclA glycoprotein has been described as a component of the exosporium; this paper defines 10 more tightly associated proteins from the exosporium of B. cereus ATCC 10876, identified by N-terminal sequencing of proteins from purified, washed exosporium. Likely coding sequences were identified from the incomplete genome sequence of B. anthracis or B. cereus ATCC 14579, and the precise corresponding sequence from B. cereus ATCC 10876 was defined by PCR and sequencing. Eight genes encode likely structural components (exsB, exsC, exsD, exsE, exsF, exsG, exsJ, and cotE). Several proteins of the exosporium are related to morphogenetic and outer spore coat proteins of B. subtilis, but most do not have homologues in B. subtilis. ExsE is processed from a larger precursor, and the CotE homologue appears to have been C-terminally truncated. ExsJ contains a domain of GXX collagen-like repeats, like the BclA exosporium protein of B. anthracis. Although most of the exosporium genes are scattered on the genome, bclA and exsF are clustered in a region flanking the rhamnose biosynthesis operon; rhamnose is part of the sugar moiety of spore glycoproteins. Two enzymes, alanine racemase and nucleoside hydrolase, are tightly adsorbed to the exosporium layer; they could metabolize small molecule germinants and may reduce the sensitivity of spores to these, limiting premature germination.     

Autoinducer 2 affects biofilm formation by Bacillus cereus

Cell-free supernatants from growing Bacillus cereus strain ATCC 10987 induced luminescence in a Photorhabdus luminescens DeltaluxS mutant, indicating the production of functional autoinducer 2 (AI-2). The exogenous addition of in vitro synthesized AI-2 had an inhibitory effect on biofilm formation by B. cereus and promoted release of the cells from a preformed biofilm.     

Localization and structural analysis of a conserved pyruvylated epitope in Bacillus anthracis secondary cell wall polysaccharides and characterization of the galactose-deficient wall polysaccharide from avirulent B. anthracis CDC 684

Bacillus anthracis CDC 684 is a naturally occurring, avirulent variant and close relative of the highly pathogenic B. anthracis Vollum. Bacillus anthracis CDC 684 contains both virulence plasmids, pXO1 and pXO2, yet is non-pathogenic in animal models, prompting closer scrutiny of the molecular basis of attenuation. We structurally characterized the secondary cell wall polysaccharide (SCWP) of B. anthracis CDC 684 (Ba684) using chemical and NMR spectroscopy analysis. The SCWP consists of a HexNAc trisaccharide backbone having identical structure as that of B. anthracis Pasteur, Sterne and Ames, →4)-β-d-ManpNAc-(1?→?4)-β-d-GlcpNAc-(1?→?6)-α-d-GlcpNAc-(1→. Remarkably, although the backbone is fully polymerized, the SCWP is the devoid of all galactosyl side residues, a feature which normally comprises 50% of the glycosyl residues on the highly galactosylated SCWPs from pathogenic strains. This observation highlights the role of defective wall assembly in virulence and indicates that polymerization occurs independently of galactose side residue attachment. Of particular interest, the polymerized Ba684 backbone retains the substoichiometric pyruvate acetal, O-acetate and amino group modifications found on SCWPs from normal B. anthracis strains, and immunofluorescence analysis confirms that SCWP expression coincides with the ability to bind the surface layer homology (SLH) domain containing S-layer protein extractable antigen-1. Pyruvate was previously demonstrated as part of a conserved epitope, mediating SLH-domain protein attachment to the underlying peptidoglycan layer. We find that a single repeating unit, located at the distal (non-reducing) end of the Ba684 SCWP, is structurally modified and that this modification is present in identical manner in the SCWPs of normal B. anthracis strains. These polysaccharides terminate in the sequence: (S)-4,6-O-(1-carboxyethylidene)-β-d-ManpNAc-(1?→?4)-[3-O-acetyl]-β-d-GlcpNAc-(1?→?6)-α-d-GlcpNH(2)-(1→.     

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.     

Purification and structure elucidation of the N-acetylbacillosamine-containing polysaccharide from Bacillus licheniformis ATCC 9945.

The exopolysaccharide of Bacillus licheniformis ATCC 9945 (formerly B. subtilis ATCC 9945) contains among other glycoses 4-acetamido-2-amino-2,4,6-trideoxy-D-glucose, termed N-acetylbacillosamine (Bac2N4NAc). A similar diamino glycose, 2-acetamido-4-amino-2,4,6-trideoxy-D-glucose, was found in a surface layer (S-layer) glycoprotein preparation of Clostridium symbiosum HB25. Electron microscopic studies, however, showed that B. licheniformis ATCC 9945 is not covered with an S-layer lattice, indicating that the N-acetylbacillosamine present in that organism might be a constituent of a cell wall-associated polymer. For elucidation of the structure of the N-acetylbacillosamine-containing polysaccharide, it was purified from a trichloroacetic acid extract of B. licheniformis ATCC 9945 cells. Using different hydrolysis protocols and a hydrolysate of the S-layer glycoprotein preparation from C. symbiosum HB25 as reference, the purified polysaccharide was found to contain 2,4-diamino-2,4,6-trideoxy-glucose, 2-acetamido-2-deoxy-glucose, 2-acetamido-2-deoxy-galactose and galactose in a molar ratio of 1 : 1 : 1 : 2. One- and two-dimensional NMR spectroscopy, including 800 MHz proton magnetic resonance measurements, in combination with chemical modification and degradation experiments, revealed that the polysaccharide consists of identical pyruvylated pentasaccharide repeating units with the structure: [-->3)-[(S)Py-(3,4)-beta-D-Galp-(1-->6)]-alpha-D-GlcpNAc-(1-->3)-beta-D-Bacp2N4NAc-(1-->3)-[(S)Py-(3,4)-beta-D-Galp-(1-->6)]-beta-D-GalpNAc-(1-->](n)     

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.