|Home | About | Journals | Submit | Contact Us | Français|
Bacillus megaterium is deep-rooted in the Bacillus phylogeny, making it an evolutionarily key species and of particular importance in understanding genome evolution, dynamics, and plasticity in the bacilli. B. megaterium is a commercially available, nonpathogenic host for the biotechnological production of several substances, including vitamin B12, penicillin acylase, and amylases. Here, we report the analysis of the first complete genome sequences of two important B. megaterium strains, the plasmidless strain DSM319 and QM B1551, which harbors seven indigenous plasmids. The 5.1-Mbp chromosome carries approximately 5,300 genes, while QM B1551 plasmids represent a combined 417 kb and 523 genes, one of the largest plasmid arrays sequenced in a single bacterial strain. We have documented extensive gene transfer between the plasmids and the chromosome. Each strain carries roughly 300 strain-specific chromosomal genes that account for differences in their experimentally confirmed phenotypes. B. megaterium is able to synthesize vitamin B12 through an oxygen-independent adenosylcobalamin pathway, which together with other key energetic and metabolic pathways has now been fully reconstructed. Other novel genes include a second ftsZ gene, which may be responsible for the large cell size of members of this species, as well as genes for gas vesicles, a second β-galactosidase gene, and most but not all of the genes needed for genetic competence. Comprehensive analyses of the global Bacillus gene pool showed that only an asymmetric region around the origin of replication was syntenic across the genus. This appears to be a characteristic feature of the Bacillus spp. genome architecture and may be key to their sporulating lifestyle.
Bacillus megaterium was first described by Anton De Bary more than 1 century ago in 1884 (14). Named for its large size, a “megat(h)erium” (Greek for big animal) of 1.5 by 4 μm, this microorganism is the largest of all bacilli. Long before Bacillus subtilis was introduced as a Gram-positive model organism, B. megaterium was used for studies on biochemistry as well as bacteriophages (13). The French microbiologist Maurice Lemoigne in 1925 discovered the polyester polyhydroxybutyrate in B. megaterium as an important energy storage molecule in bacteria (32), and Andre Lwoff discovered UV induction of bacteriophage in a lysogenic B. megaterium strain (35). Due to its large cell size, B. megaterium is well-suited for research on cell morphology, such as cell wall and cytoplasmic membrane biosynthesis, sporulation, spore structure and cellular organization, DNA partitioning, and protein localization (10, 61). In the 1960s, B. megaterium was used to study sporulation, since it sporulates and germinates efficiently (18). Because of its biotechnological use in the production of several substances, the nonpathogenic B. megaterium is of general interest to industry (7). In contrast to Gram-negative organisms like Escherichia coli, B. megaterium does not produce endotoxins associated with the outer membrane, which, combined with its growth on a variety of carbon sources and simple media, has made it a workhorse in food and pharmaceutical production processes for decades (i.e., α- and β-amylases used for starch modification in the baking industry and penicillin acylases essential for the synthesis of novel β-lactam antibiotics, among others ). Further, it is one of the most efficient producers of vitamin B12 (63, 64). B. megaterium has been extensively studied genetically and is amenable to genetic manipulation (63, 66). Hundreds of auxotrophs, division mutants, antibiotic-resistant, and UV-sensitive mutants have been characterized and have been previously mapped in QM B1551 (17, 20, 31, 55, 57, 62, 65). Strain QM B1551 is second only to B. subtilis in the number of multiply marked, characterized strains that are available from the Bacillus Genetic Stock Center (Vary collection, BGSC, Ohio State University; http://www.bgsc.org/). Several biotechnological mutants have been constructed in strain DSM319 and are commercially available (38, 40, 72). Recent analysis using phylogenetic analysis based on 16S rRNA genes led to the division of the genus Bacillus into 4 different families and 37 genera (34). Even after the taxonomic reorganization, Bacillus is a diverse genus with G+C content ranging from 34 to 35% (Bacillus cereus and related pathogens) to 44% (B. subtilis) and species that differ radically in lifestyles and metabolic properties. Most of the sequenced Bacillus genomes are closely related to B. cereus or B. subtilis. Recently, Porwal et al. showed that B. megaterium (with a G+C content of 38 to 39%) is only distantly related to the B. cereus and B. subtilis groups and that it is more deeply rooted in the phylogenetic tree than previously thought (44, 51). To gain insights into the genome evolution and the metabolic versatility that facilitate biotechnological applications, we have sequenced the complete genomes of B. megaterium strains QM B1551 and DSM319. We have used this information to examine the genetic diversity, genome dynamics, and phylogenetic relationships within B. megaterium and among members of the genus in great detail. Our genomic analysis reveals new genetic traits not previously seen in Bacillus and has resulted in a refined model for genome evolution and adaptation of B. megaterium.
Freeze-dried spores of B. megaterium QM B1551 prepared in 1967 by James C. Vary were provided by Patricia S. Vary (Northern Illinois University). These spores are closest in time to the original QM B1551 (Quarter Master Bacterium 1551; ATCC 12872) isolation by Hillel Levinson at the Pioneering Research Division, Quartermaster and Engineering Center, Natick, MA. B. megaterium DSM319 is a naturally plasmidless strain isolated by Stahl and Esser (53) and obtained from DSMZ (http://www.dsmz.de/). The plasmidless QM B1551 derivative, B. megaterium strain PV361 (56), was used for high-throughput phenotype arrays. PV361 was grown on rich and minimal media and showed no differences in sporulation compared to QM B1551. Germination of PV361 was performed on rich medium because it lacks pBM700, which carries key germination genes (55). Strain PV586, a Lac− derivative of PV361 (29), was used as transformation recipient in the replicon analysis of pBM600.
Detailed materials and methods describing genome sequencing, assembly, annotation, comparative genome analysis, pangenome computation, and the different experimental assays can be found in the supplemental material.
The B. megaterium complete genome sequences (QMB1551 project ID 30165; DSM319 project ID 33377) have been deposited in the NCBI GenBank under accession numbers CP001983 (chromosome, QM B1551), CP001984 to CP001990 (pBM100 to pBM700, respectively), and CP001982 (chromosome, DSM319).
The chromosomes of B. megaterium strains QM B1551 and DSM319 are circular molecules of 5,097,129 bp and 5,097,447 bp, respectively, with an average G+C content of 38.2%. A well-defined G+C bias was observed, with the right-hand replichore significantly enriched in G+C relative to the left-hand replichore (Fig. 1 and Table 1). In addition, strain QM B1551 contains seven indigenous plasmids, pBM100 to pBM700, with sizes from 5.4 kb to over 164 kb, while strain DSM319 is a naturally plasmidless isolate (53) (Fig. 2 and Table 2; see also Fig. S1 in the supplemental material). The plasmids have significantly lower G+C contents than the chromosomes (33.0 to 36.5 versus 38.2%). The chromosome of strain QM B1551 contains 5,284 genes, and that of strain DSM319 contains 5,272 genes. The two chromosomes display a high level of genome conservation, having a nucleotide sequence identity of more than 95% over 83.3% of the length of the two chromosomes. The genomes are mostly colinear (Fig. 3), with only a single rearrangement larger than three genes: a 17-gene block in QM B1551 (Fig. 3, red arrows) is inverted and displaced by 1.7 Mbp relative to DSM319 (genes BMQ_1765 to BMQ_1781 are homologous with genes BMD_3632 to BMD_3616). The rearranged region does not appear to be a single functional unit: genes are on both strands and are involved in several different biological processes. Most of the genetic differences between the two genomes are due to insertions or deletions (indels) of single genes or small groups of genes at scattered and independent genomic locations throughout both chromosomes. We cataloged a total of 300 and 254 isolate-specific genes that are organized in 96 and 106 independent clusters for strains QM B1551 and DSM319, respectively (see Table S1 in the supplemental material). The distribution of these indels is not random (Fig. 1, circles 8 and 9): we found that the gene insertions in both chromosomes are rare in a 2-Mb region around the origin of replication (ori) (Fig. 1 and Fig. 3; see also Table S1). Compared to the function of genes common to both strains, strain-specific genes are increased in functions affecting interactions with the environment: cell envelope, transport, signal transduction, and gene regulation (Fig. 4; see also Fig. 9, below). In contrast, relatively fewer strain-specific genes are associated with basic cellular processes, such as amino acid, nucleotide, or cofactor biosynthesis, central intermediary metabolism, fatty acid metabolism, protein synthesis or degradation, or transcription. The isolate-specific genes are also more frequently annotated as conserved hypothetical genes or genes coding for enzymes of unknown specificity.
Plasmids make up 11% of the QM B1551 genome (24). This strain harbors seven indigenous plasmids with plasmid copy numbers ranging between 1 and 18 copies (Table 2). Since many other strains of B. megaterium carry multiple plasmids (63), this is a critical part of the genome analysis for this species. The three largest plasmids are shown in Fig. 2. Plasmids pBM100 and pBM200 of QM B1551 were previously sequenced and deposited in GenBank. Plasmids pBM300 and pBM400 have been analyzed (30, 52, 63) and have been resequenced to higher coverage and reannotated for this study (see Fig. S1 in the supplemental material). The plasmids carry a variety of genes, including genes for sporulation, germination, regulation, transport, and lantibiotic synthesis, as well as erythromycin and rifampin resistance. There are also genes for fatty acid metabolism, cell wall hydrolysis, sigma factors, and cell division, as well as integrons, insertion sequence (IS) elements, and transposons. Because B. megaterium is frequently found with Pseudomonas species in contaminated environments, it has long been suspected of having the capability of metabolizing unusual substrates of possible bioremedial use (30, 52, 61). To that extent, B. megaterium plasmids carry genes for heavy metal resistance, including Cu and Cd export, and genes such as styrene monooxygenase are present. Several metabolic genes on the larger plasmids are organized in what appear to be functional operon structures and may enable this strain to survive in unusual habitats (Fig. 2). The possible substrates for associated transporters cannot be deduced from genomic annotation, but the identification of these regions provides the basis for further experimental characterization.
Of special note is that all seven indigenous plasmids appear to be unique to B. megaterium compared to a broad panel of diverse Bacillus spp. plasmids. Besides similarities in insertion sequence genes among the plasmids, we found three genes similar to the pXO1 virulence plasmid of B. anthracis, including a reverse transcriptase (52) (see Fig. S1 in the supplemental material). This finding may indicate either a common phylogenetic origin or genetic exchange among these plasmid-borne genes. Plasmid pBM500 (Fig. 2) is intriguing because of the presence of three possible sigma factors (BMQ_pBM50015, BMQ_pBM50023, and BMQ_pBM50090), a collagen-like gene (BMQ_pBM50081), and an erm(B) resistance gene (BMQ_pBM50045). In addition, a cytochrome P450 gene (BMQ_pBM50008) is plasmid-borne. The P450 enzymes of B. megaterium have long been used as a model system (15, 20, 21). There is a gene cluster that appears to be responsible for the biosynthesis of a fatty acid compound (BMQ_pBM50048 to BMQ_pBM50053). Plasmid pBM600 seems to be a depository for mobile genetic elements, since it carries four transposases, three integrase genes, and a group II intron. It also has eight glycosyltransferases and other glucose carrier proteins within a 25-kb region (kb 13.9 to 39.6). Plasmid pBM700 encodes several transferases and transporters, transcriptional regulators, and bacteriocin and antibiotic synthesis operons. In addition, two two-component histidine kinase regulatory systems are present (BMQ_pBM70140/141 and BMQ_pBM70155/156), strengthening the role of this plasmid in metabolic activities and strain-specific niche adaptations. Previous studies on the pBM700 germination (ger) operons have helped identify specific amino acid residues that define receptor specificity to germinant compounds in Bacillus (11, 12). Strains of QM B1551 cured of pBM700 can no longer germinate on single germinant compounds (55), and the complete germination operon, with gerUA, gerUC, and gerUB and including the downstream monocistronic gene gerVB (BMQ_pBM70070 to BMQ_pBM70073), has been identified (10, 12).
With the exception of pBM600 (Fig. 2), all replicons have been previously identified and functionally analyzed (30, 52, 64). A region containing a putative replicon gene, BMQ_pBM60001, was identified and cloned in pJM103, a plasmid that cannot replicate in B. megaterium. The clone was successfully maintained in B. megaterium PV586, indicating that it is the pBM600 replicon. The presence of a characteristic plasmid replication motif, including six copies of a 21-bp iteron upstream from BMQ_PBM60001, indicates it is most likely a theta replicon. However, the identified replicon shows no homology with the other four theta plasmid replicons in QM B1551. We note that orthologs of QM B1551 plasmid genes are found not only on the genome of QM B1551 but also on DSM319 (Fig. 1; see also Table S2 in the supplemental material). This observation indicates that strain DSM319 is potentially capable of acquisition and maintenance of plasmids.
We have observed extensive gene exchange between QM B1551 plasmids and the main chromosomes of both strains (Fig. 1, circles 10 and 11). Of the QM B1551 499 protein-coding plasmid genes, 104 have homologs (BSR, >0.4) on either the QM B1551 or DSM319 chromosome (see Table S2). Some genes share more than 90% sequence identity over the full peptide length with their chromosomal homologs. Plasmid genes with homologs on the chromosomes are found on all plasmids except pBM100 (5.4 kb). The majority of genes are encoded on pBM400, pBM500, and pBM700, while only a few genes are found on pBM200, pBM300, and pBM600 (see Table S2). These genes can be classified into three groups. There are 19 genes present on QM B1551 plasmids that have homologs on the QM B1551 chromosome but not on the DSM319 chromosome. Conversely, 32 plasmid genes have homologs on the DSM319 chromosome but not the QM B1551 chromosome. The remaining 53 genes have homologs on both the QM B1551 and DSM319 chromosomes (see Table S2). These findings demonstrate that some of the plasmid-borne capabilities of QM B1551 are also found in the plasmidless DSM319 strain. The fact that the number of genes shared between the plasmids and QM B1551 chromosome is smaller than the number shared with the DSM319 chromosome suggests that the plasmidless state of DSM319 might be the result of a recent event. The homologous genes are found scattered in different locations, both on the plasmids and on the DSM319 and QM B1551 chromosomes. Several of the genes are found in clusters of 2 to 11 genes and appear to have been cotransferred, as they have the same relative positions and orientations on both the plasmid and chromosomes with conserved intergenic nucleotide homology. This set of genes encodes a variety of physiological functions, such as metabolism, transport, regulation, sporulation, and germination, and include 45 conserved hypothetical CDS with no assigned function.
We detected major differences in the genes involved in spore coat structure and its biosynthesis between B. megaterium and B. subtilis strain 168 (see Table S3 in the supplemental material). More than half of the B. subtilis genes involved in spore coat formation are missing and could not be identified in either B. megaterium strain, which suggests that B. megaterium possesses a different spore coat structure than B. subtilis (23).
Transformation competence has been well-studied in the Gram-positive model organism B. subtilis 168 (16). To date, no conditions have been found that induce natural competence in B. megaterium. B. megaterium genomes contain 33 competence orthologs to genes known to mediate competence in B. subtilis (see Table S4 in the supplemental material) (28) and that may allow B. megaterium to incorporate foreign genetic material. The genome data revealed that B. megaterium lacks the comQXPA gene cluster and two other competence genes, comFB and comS, that are present in B. subtilis. Having detected almost all known genes for competence in other bacilli (see Table S4), we speculate that competence would be easy to engineer in the well-studied B. megaterium strains. Furthermore, it may occur naturally in some strains of the species.
B. megaterium QM B1551 is motile but requires an oxygen-permeable coverslip for sustained motility (P. S. Vary, unpublished data). The flagellar genes of both QM B1551 and DSM319 form two clusters that are almost identical to those found in B. subtilis (Fig. 1 and and5D).5D). Cluster one corresponds to the fla/che 26-kb operon in B. subtilis (42). In this cluster, only CheC, a signal-terminating phosphatase important in chemotaxis protein methylation, is missing from the operon and not found elsewhere in B. megaterium. The second and larger cluster is more complex (Fig. 5D). The fliT, fliS, and fliD genes match those of B. subtilis, but flaG is not found in B. megaterium. Moreover, the B. subtilis gene hag, which codes for flagellin, is not part of the cluster but is found almost 2 Mb away (BMQ_1093), not near any other flagellar genes except an adjacent short flagellin domain gene (BMQ_1092), which is presumably a gene fragment that may have been created during the transposition events that moved the hag gene to its new location. In its place in the cluster are six genes that apparently mediate unrelated metabolic functions (BMQ_5103 to BMQ_5108). These genes code for a GABA permease and associated symporter, a betaine aldehyde dehydrogenase, an alcohol dehydrogenase, and a GbsR transcriptional regulator, as well as a hypothetical protein. The syntenic genes consist of flagellar genes (fliWLKS), two competence-associated genes (comFA and comFC), and the degS-degU genes, which code for a two-component system that regulates the transition from exponential to stationary growth. We note here that DegU has also been shown to regulate motility via σD in B. subtilis (36).
The genomic sequences were used to investigate the underlying molecular nature of some of the unique physiological and metabolic features of B. megaterium. The transporter inventory of both B. megaterium strains and potential substrates (http://www.membranetransport.org/index.html) has been analyzed using TransportDB (50).
B. megaterium carries genes for the glyoxylate pathway. In keeping with its known phenotype, B. megaterium uses the Embden-Meyerhof pathway for glycolysis, followed by a classical Krebs cycle. In contrast to other bacilli, B. megaterium contains an intact glyoxylate pathway (see Fig. S2A in the supplemental material). From this starting point, the B. megaterium genome contains genes for the synthesis of all 20 proteinogenic amino acids, all nucleotides and cofactors used by other enzymes, and all necessary fatty acids. Genes for the transporters and special enzymes needed to utilize a wide variety of carbon sources have been found. Aerobic respiration involves three membrane-bound dehydrogenases to transfer electrons from NADH, glycerol-3-phosphate, and succinate to menaquinone (see Fig. S2B). Menaquinone is oxidized by two menaquinone oxidases of the cytochrome bd type, by a cytochrome aa3-type oxidase, which is not found in other bacilli, or by the cytochrome bc1 complex, followed by cytochrome c and cytochrome c oxidase. B. megaterium also contains a fourth oxidase with high sequence similarity to E. coli cytochrome o ubiquinol oxidase, which is not found in other bacilli. Under anaerobic conditions, B. megaterium has all the genes required to perform a mixed acid fermentation that produces lactate, 2,3-butanediol, and acetate, as previously described for B. subtilis (see Fig. S2C) (41, 47). Unlike most other bacilli, B. megaterium is not capable of using nitrate as an electron acceptor, since it lacks a membrane-bound nitrate reductase of the Nar type. Further details about B. megaterium energy metabolism are found in Fig. S2 of the supplemental material.
Predictions for the metabolic capabilities of each strain were deduced from in silico comparative genome analysis and validated by growth experiments in vivo (Fig. 6). For strain QM B1551, several unique glucarate-degrading enzymes were found in an operon-like structure, including glucarate permease (BMQ_1892), glucarate dehydratase (BMQ_1893), and the degradative enzymes galactarate dehydratase (BMQ_1888) and 5-dehydro-4-deoxyglucarate dehydratase (BMQ_1890). The final product of this glucarate degradation process is 2,5-dioxopentanoate, which is discharged into the citrate cycle via 2-oxoglutarate. Furthermore, a complete galactitol-specific phosphotransferase system (PTS) is found to be unique in B. megaterium QM B1551 (BMQ_3201 to BMQ_3204). The genes for the corresponding enzymes galactitol-1-phosphate 5-dehydrogenase (BMQ_3200) and an alcohol dehydrogenase (BMQ_3199) are also part of this genomic island. We demonstrated experimentally that B. megaterium QM B1551 (in contrast to strain DSM319) can utilize both glucarate and galactitol as carbon sources (Fig. 6A and B). Additionally, growth experiments on pullulan showed that B. megaterium QM B1551 grows more effectively on pullulan than does DSM319 (Fig. 6C), although both strains are equipped with the same chromosomal pullulanases, AmyX (BMQ_4837) and PulA (BMQ_2037). Strain DSM319 encodes a cellulase-like glycosyl hydrolase (BMD_1113), an additional citrate two-component system (BMD_2292/2293), part of a sucrose-specific PTS (BMD_3721), sorbitol dehydrogenase GutB (BMD_3588) with its transcriptional activator GutR (BMD_3589), and a probable sorbitol transporter (BMD_3587). The sorbitol uptake system proved to be very effective in DSM319 (Fig. 6D), with optical growth densities of 6, while that for QM B1551 was only 1.2.
In its varied habitats, B. megaterium is exposed to increased osmolarity that triggers water efflux from the cell. It therefore has to adjust the osmotic potential of its cytoplasm to prevent dehydration. The physiology and genetics of osmotic adjustment to high-osmolarity surroundings has been studied intensively in B. subtilis (6). In the initial phase of adjustment to high osmolarity, B. subtilis accumulates large amounts of potassium via the KtrAB and KtrCD systems to curb the efflux of water (22). Both Ktr potassium uptake systems are also present in B. megaterium (BMQ_3903/3904 and BMQ_1339/1233). In the second phase of osmotic adjustment, B. subtilis synthesizes large quantities of the compatible solute proline and acquires various types of organic osmoprotectants (e.g., glycine betaine) from environmental sources via high-affinity transport systems. Recent studies have shown that B. megaterium DSM32 synthesizes proline when it is challenged by high salinity (9). The in silico analysis of the genome sequences of QM B1551 and DSM319 revealed a complete proline biosynthetic gene cluster (proHJA; BMQ_2287 to BMQ_2289) whose counterparts are osmotically induced in B. subtilis, B. licheniformis, and Halobacillus halophilis (26, 71). B. megaterium may also use an exogenous supply of proline for osmoprotective purposes, since it possesses a homolog (BMQ_1420) of the osmotically induced proline uptake system OpuE from B. subtilis (6). Glycine betaine is a very effective osmoprotectant, and B. subtilis can acquire it either from environmental sources or synthesize it from an exogenous supply of the precursor choline via a two-step oxidation process. Based on homology to B. subtilis, several glycine betaine import systems were identified in B. megaterium. Two ABC transporters of the OpuA type, opuACAB (BMQ_0858 to BMQ_0860) and opuAABC (BMQ_1542 to BMQ_1544), an OpuC-like glycine betaine/choline importer (BMQ_3925/3926), and an OpuD-like BCCT-type transport system (BMQ_1351) are present. In addition, a gene cluster encoding the glycine betaine biosynthetic enzymes GbsA and GbsB (BMQ_5106/5107) are present.
Besides its most common soil habitat, B. megaterium is also found in diverse environments, including rice paddies, dried food, seawater, sediments, fish, and even in bee honey (63). Analysis of B. megaterium genome sequences strongly indicates that it is well-adapted to an aquatic lifestyle, since it possesses genes for the formation of gas vesicles. These intracellular structures are often found in marine microorganisms and function as flotation devices, allowing the bacterium to float up and down by adjusting the gas pressure inside the gas vesicle (60). The genome of B. megaterium QM B1551 has two gas vesicle-forming operons (Fig. 5A). The two gvp loci share eight homologous genes (gvpBRFGLSKJ) and are syntenic except for gvpN, which is found between gvpR and gvpF in the larger operon. The smaller locus (BMQ_3224 to BMQ_3231) is not present in DSM319. The larger operon encompasses 14 genes, gvpAPQBRNFGLSKJTU (BMQ_3290 to BMQ_3303), and an araC-like regulator. This gene cluster was previously discovered in B. megaterium strain VT1660, and its functional expression in E. coli renders this bacterium buoyant (33). It should be noted that genes for gas vesicles have been found in genomes of many microorganisms that do not typically live in aquatic habitats, fostering discussions on the physiological function of these floating devices (60). The proteins involved in gas vesicle formation in B. megaterium show a phylogenetic relationship far outside the Bacillus group and are related to proteins present in aquatic Cyanobacteria but also in members of Archaea. One explanation for the gas vesicles could be related to the ability of the bacterium to leave aquatic regions of low oxygen tension.
One major advantage of B. megaterium for industrial protein production is its ability to secrete proteins directly into the growth medium (2, 4, 37, 45). B. megaterium employs Sec and Tat systems for this purpose; the relevant genes were all detected in the genomes. The Sec-dependent secretion system is especially widely used and has been investigated for industrial recombinant protein production. A look at both genome sequences of B. megaterium revealed potential genes for approximately 30 secreted proteases (see Table S5 in the supplemental material). Only a few proteases have experimentally been found to be secreted, including metalloprotease InhA, extracellular protease Vpr, aminopeptidase YwaD, and neutral protease NprM (67). The gene for the major extracellular protease NprM was deleted from production strain MS941 (69), which allowed up to 1.2 g/liter functional proteins to be recovered from the medium (4, 20, 25, 54, 66). These data clearly demonstrate the advantage of B. megaterium strains for recombinant protein production.
B. megaterium was one of the first biotechnological vitamin B12 producers described for bacteria (61, 64, 70). In agreement with the well-studied biosynthetic pathway in Salmonella enterica, which is known for its ability to synthesize vitamin B12 in the presence and absence of oxygen (5, 8, 49), the genes for oxygen-independent vitamin B12 biosynthesis genes were found in the genomes of B. megaterium. They are organized in two distinct, independent operons (Fig. 5B). The larger operon, cbiWHXJCDETLFGA-cysG-cbiY-btuR, has been thoroughly described (49); the smaller operon, cbiB-cobDUSC, is associated with an uncharacterized ATP:cob(I)alamin adenosyltransferase and codes for the enzymes involved in the final steps of the biosynthetic pathway (Fig. 7). The gene cbiP coding for adenosylcobyrinic acid synthase (EC 184.108.40.206) is isolated in the genome at Mb 2.9 as a single gene with no further B12 biosynthesis context. The genome of both B. megaterium strains now allows the full reconstruction of the vitamin B12 biosynthetic pathway (Fig. 7).
Several enzymes with dependencies on vitamin B12 have been identified in silico in the B. megaterium genome. Ethanolamine ammonia-lyase EutBC (EC 220.127.116.11) and methylmalonyl coenzyme A (CoA) mutase MutAB (EC 18.104.22.168) are two that have been shown to contribute to the ability of B. megaterium to degrade lipids. Other than B. megaterium, only Geobacillus kaustophilus and B. halodurans have been found to carry the MutAB genes. EutBC is part of a large B. megaterium species-specific ethanolamine utilization operon, eutHSPABCLEM (BMQ_3678 to BMQ_3686) (Fig. 5C), together with an upstream ethanolamine two-component response regulator system (BMQ_3687/3688) and further uncharacterized genes contributing to ethanolamine utilization (BMQ_3675 to -3677). A second EutBC can also be found as part of a small eutABC-eat cluster (BMQ_2531 to -2534), where Eat is an ethanolamine permease. It is noteworthy that a phage integrase (BMQ_2528) is found three genes downstream, which may indicate a lateral acquisition of this functional unit. Furthermore, the vitamin B12-dependent methionine synthase MetH (BMQ_1293) and a vitamin B12-dependent ribonucleoside-diphosphate reductase (BMQ_4885/4886) were identified in both strains.
B. megaterium belongs to a deeply rooted lineage within the genus Bacillus, based on its phenotypic characteristics, genome size, intermediate G+C content, and 16S rRNA phylogeny (Table 1) (1, 46). Whether B. megaterium is more closely related to B. cereus and its relatives or to B. subtilis and its relatives has been difficult to resolve based on 16S rRNA gene phylogeny alone. Trees showing all three possible relationships between these groups have been published recently (27, 44, 68). To address this issue, we performed a phylogenetic analysis on 385 orthologous genes identified using BSR analysis (see Fig. S3 in the supplemental material) (48). A neighbor-joining tree was generated, using Listeria monocytogenes as the outgroup. In this tree B. megaterium, B. cereus, and B. subtilis each anchor clades that are supported in all bootstrap samples. However, the order with which these clades join varies between bootstrap samples, with each of the three possible orders supported by 21 to 51% of the bootstrap samples. Thus, the ambiguous relationship between B. megaterium, B. cereus, and B. subtilis that is seen with 16S rRNA phylogenies is shown to also exist when a large set of orthologous genes is examined, possibly indicating a large amount of lateral gene transfer between these species.
When orthologous genes were compared between B. megaterium and other Bacillus species, a region of striking synteny was observed neighboring the ori (Fig. 1), from about the 10 o'clock to the 2 o'clock position (Fig. 8; see also Fig. S4 in the supplemental material). The Gram-positive spore former Oceanobacillus iheyensis and Listeria monocytogenes, the causative agent of listeriosis, show similar regions of synteny. Increased phylogenetic distance generally correlated with an increased number of insertions and inversions in the syntenic regions. The region of synteny is dominated by genes involved in sporulation, cell envelope biosynthesis, and the transcription and translation machinery. This phenomenon might be explained by a need for these genes to be readily accessible during germination and sporulation (73).
The availability of two complete genomes of B. megaterium (Fig. 1; see also Table S1 in the supplemental material) provides the opportunity to investigate the genomic plasticity and global gene reservoir of the Bacillus group of organisms, the Bacillus pan-genome (39) (Fig. 9). The core genome, genes expected to be present in all genomes, is estimated at 2,009 protein-coding genes (Fig. 9A). We estimated the Bacillus pan-genome to contain 10,534 unique protein-coding genes (Fig. 9C). The pan-genome defines the total number of genes found at least once among all Bacillus genomes. Our estimates of the gene discovery rate indicate that Bacillus has an open pan-genome: on average, 96 new genes would be expected from each new Bacillus genome sequenced (Fig. 9B). These results are comparable to data obtained for E. coli (59) or for Staphylococcus agalacticae and S. pneumoniae (58). As evidenced by the genomic analyses of B. megaterium QM B1551 and DSM319, there is a high degree of genetic diversity even among closely related strains. The gene discovery rate in Bacillus was significantly increased by the fact that each strain had about 300 genes not found in the other strain, with an additional 499 genes found on the QM B1551 plasmids. In a species featuring an open pan-genome, such as E. coli, the predicted pan-genome size is almost 75% larger than the average genome size of the species (59), a finding comparable to our results within the Bacillus group of organisms.
The analysis of the genome sequences of the historically and biotechnologically important strains QM B1551 and DSM319 has advanced our knowledge of the metabolic versatility and the evolution of the Gram-positive aerobic spore-forming bacteria. We have identified numerous unique genetic traits not seen in any Bacillus species previously studied, including the presence of gas vesicle proteins, a glyoxylate pathway, and vitamin B12 biosynthesis, and also a lack of homology with the spore coats of B. subtilis. Our data suggest there is considerable DNA exchange between Bacillus and other Gram-positive bacteria (including pathogens), as well as on the intraspecies level between the seven indigenous plasmids and chromosomes. The observed phenomenon of large conserved syntenic regions neighboring the chromosomal ori is key in understanding the lifestyle of B. megaterium and other related spore formers. The availability of the genome sequences of these two B. megaterium strains, the reconstruction of key metabolic and energetic pathways, and the finding that most of the known Bacillus species competence genes are present, together with well-established protocols for genetic modifications, offer a promising future for the further development of B. megaterium as a model organism for systems biotechnology (3, 19, 64).
The sequencing of B. megaterium strains QM B1551 and DSM319 was supported with federal funds from the National Science Foundation under NSF contract no. 0802327 and the German Research Foundation under contract no. SFB578.
We thank Vanessa Hering for the growth experiments, Anni Moore and Melvin Duvall for phylogeny discussions, and Stefan Münnich and students of the Bioinformatics classes at NIU for help with the annotation.
We declare that we have no competing financial interests.
†Supplemental material for this article may be found at http://jb.asm.org/.
Published ahead of print on 24 June 2011.