Nucleotide Sequence and Features of the Bacillus licheniformis gnt Operon

Bacillus licheniformis was able to utilize gluconate as thesole carbon source as efficiently as Bacillus subtilis did.Southern analysis indicated that B. licheniformis likely possessesonly one gnt determinant. The nucleotide sequence (6278 bp)of the B. licheniformis DNA containing the gnt operon was determined,revealing the five complete open reading frames (ORF; genes).The putative product of the first gene, oug, did not show anysignificant homology to known proteins, but those of the secondto fifth genes exhibited striking homology to the gntRKPZ genesof B. subtilis, respectively, indicating that they are the correspondinggnt genes of B. licheniformis. Not only is the organizationof the gnt genes of these two Bacilli highly conserved, butso are the cis regulatory elements of their gnt operon. Sequenceanalysis of the upstream regions of these two gnt operons impliedthat a chromosome rearrangement in B. subtilis might have occurredimmediately upstream of the gnt operon during evolution, causingit to diverge from a common ancestor into B. licheniformis andB. subtilis.     

Production of the Novel Two-Peptide Lantibiotic Lichenicidin by Bacillus licheniformis DSM 13

Background

Lantibiotics are small microbial peptide antibiotics that are characterized by the presence of the thioether amino acids lanthionine and methyllanthionine. Lantibiotics possess structural genes which encode inactive prepeptides. During maturation, the prepeptide undergoes posttranslational modifications including the introduction of rare amino acids as lanthionine and methyllanthione as well as the proteolytic removal of the leader. The structural gene (lanA) as well as the other genes which are involved in lantibiotic modification (lanM, lanB, lanC, lanP), regulation (lanR, lanK), export (lanT(P)) and immunity (lanEFG) are organized in biosynthetic gene clusters.

Methodology/Principal Findings

Sequence comparisons in the NCBI database showed that Bacillus licheniformis DSM 13 harbours a putative lantibiotic gene cluster which comprises two structural genes (licA1, licA2) and two modification enzymes (licM1, licM2) in addition to 10 ORFs that show sequence similarities to proteins involved in lantibiotic production. A heat labile antimicrobial activity was detected in the culture supernatant and a heat stabile activity was present in the isopropanol cell wash extract of this strain. In agar well diffusion assays both fractions exhibited slightly different activity spectra against Gram-positive bacteria. In order to demonstrate the connection between the lantibiotic gene cluster and one of the antibacterial activities, two Bacillus licheniformis DSM 13 mutant strains harbouring insertions in the structural genes of the modification enzymes licM1 and licM2 were constructed. These strains were characterized by a loss of activity in the isopropanol extract and substractive MALDI-TOF predicted masses of 3020.6 Da and 3250.6 Da for the active peptides.

Conclusions/Significance

In conclusion, B. licheniformis DSM 13 produces an antimicrobial substance that represents the two-peptide lantibiotic lichenicidin and that shows activity against a wide range of Gram-positive bacteria including methicillin resistant Staphylococcus aureus strains.     

Regulation of Polyglutamic Acid Synthesis by Glutamate in Bacillus licheniformis and Bacillus subtilis

The synthesis of polyglutamic acid (PGA) was repressed by exogenous glutamate in strains of Bacillus licheniformis but not in strains of Bacillus subtilis, indicating a clear difference in the regulation of synthesis of capsular slime in these two species. Although extracellular γ-glutamyltranspeptidase (GGT) activity was always present in PGA-producing cultures of B. licheniformis under various growth conditions, there was no correlation between the quantity of PGA and enzyme activity. Moreover, the synthesis of PGA in the absence of detectable GGT activity in B. subtilis S317 indicated that this enzyme was not involved in PGA biosynthesis in this bacterium. Glutamate repression of PGA biosynthesis may offer a simple means of preventing unwanted slime production in industrial fermentations using B. licheniformis.     

Molecular cloning,sequence analysis,and characterization of a new cell wall hydrolase,CwIL, of Bacillus licheniformis

We have cloned a DNA fragment containing the gene for a cell wall hydrolase from Bacillus licheniformis FD0120 into Escherichia coli. Sequencing of the fragment showed the presence of an open reading frame (ORF; designated as cwlL), which is different from the B. licheniformis cell wall hydrolase gene cwlM, and encodes a polypeptide of 360 amino acids with a molecular mass of 38 994. The enzyme purified from the E. coli clone is an N-acetylmuramoyl-l-alanine amidase, which has a M_r value of 41 kDa as determined by SDS-polyacrylamide gel electrophoresis, and is able to digest B. licheniformis, B. subtilis and Micrococcus luteus cell walls. The nucleotide and deduced amino acid sequences of cwlL are very similar to those of ORF3 in the putative operon xpaL1-xpaL2-ORF3 in B. licheniformis MC14. Moreover, the amino acid sequence homology of CwlL with the B. subtilis amidase CwlA indicates two evolutionarily distinguishable regions in CwlL. The sequence homology of CwlL with other cell wall hydrolases and the regulation of cwlL are discussed.     

Hyaluronic Acid Production in Bacillus subtilis

The hasA gene from Streptococcus equisimilis, which encodes the enzyme hyaluronan synthase, has been expressed in Bacillus subtilis, resulting in the production of hyaluronic acid (HA) in the 1-MDa range. Artificial operons were assembled and tested, all of which contain the hasA gene along with one or more genes encoding enzymes involved in the synthesis of the UDP-precursor sugars that are required for HA synthesis. It was determined that the production of UDP-glucuronic acid is limiting in B. subtilis and that overexpressing the hasA gene along with the endogenous tuaD gene is sufficient for high-level production of HA. In addition, the B. subtilis-derived material was shown to be secreted and of high quality, comparable to commercially available sources of HA.     

Genome-Based Identification of Active Prophage Regions by Next Generation Sequencing in Bacillus licheniformis DSM13

Prophages are viruses, which have integrated their genomes into the genome of a bacterial host. The status of the prophage genome can vary from fully intact with the potential to form infective particles to a remnant state where only a few phage genes persist. Prophages have impact on the properties of their host and are therefore of great interest for genomic research and strain design. Here we present a genome- and next generation sequencing (NGS)-based approach for identification and activity evaluation of prophage regions. Seven prophage or prophage-like regions were identified in the genome of Bacillus licheniformis DSM13. Six of these regions show similarity to members of the Siphoviridae phage family. The remaining region encodes the B. licheniformis orthologue of the PBSX prophage from Bacillus subtilis. Analysis of isolated phage particles (induced by mitomycin C) from the wild-type strain and prophage deletion mutant strains revealed activity of the prophage regions BLi_Pp2 (PBSX-like), BLi_Pp3 and BLi_Pp6. In contrast to BLi_Pp2 and BLi_Pp3, neither phage DNA nor phage particles of BLi_Pp6 could be visualized. However, the ability of prophage BLi_Pp6 to generate particles could be confirmed by sequencing of particle-protected DNA mapping to prophage locus BLi_Pp6. The introduced NGS-based approach allows the investigation of prophage regions and their ability to form particles. Our results show that this approach increases the sensitivity of prophage activity analysis and can complement more conventional approaches such as transmission electron microscopy (TEM).     

Conversion of sub-megasized DNA to desired structures using a novel Bacillus subtilis genome vector

A novel genome vector using the 4215 kb Bacillus subtilis genome provides for precise target cloning and processing of the cloned DNA to the desired structure. Each process highly dependent on homologous recombination in the host B.subtilis is distinguished from the other cloning systems. A 120 kb mouse jumonji (jmj) genomic gene was processed in the genome vector to give a series of truncated sub-megasized DNA. One of these truncated segments containing the first intron was copied in a plasmid by a recombinational transfer method developed for B.subtilis. DNA manipulation previously considered difficult is argued with respect to DNA size and accuracy.     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

Whole Genome Complete Resequencing of Bacillus subtilis Natto by Combining Long Reads with High-Quality Short Reads

De novo microbial genome sequencing reached a turning point with third-generation sequencing (TGS) platforms, and several microbial genomes have been improved by TGS long reads. Bacillus subtilis natto is closely related to the laboratory standard strain B. subtilis Marburg 168, and it has a function in the production of the traditional Japanese fermented food “natto.” The B. subtilis natto BEST195 genome was previously sequenced with short reads, but it included some incomplete regions. We resequenced the BEST195 genome using a PacBio RS sequencer, and we successfully obtained a complete genome sequence from one scaffold without any gaps, and we also applied Illumina MiSeq short reads to enhance quality. Compared with the previous BEST195 draft genome and Marburg 168 genome, we found that incomplete regions in the previous genome sequence were attributed to GC-bias and repetitive sequences, and we also identified some novel genes that are found only in the new genome.     

Comparative Analysis of Physical Maps of Four Bacillus subtilis (natto) Genomes

The complete SfiI and I-CeuI physical maps of four Bacillus subtilis (natto) strains, which were previously isolated as natto (fermented soybean) starters, were constructed to elucidate the genome structure. Not only the similarity in genome size and organization but also the microheterogeneity of the gene context was revealed. No large-scale genome rearrangements among the four strains were indicated by mapping of the genes, including 10 rRNA operons (rrn) and relevant genes required for natto production, to the loci corresponding to those of the B. subtilis strain Marburg 168. However, restriction fragment length polymorphism and the presence or absence of strain-specific DNA sequences, such as the prophages SPβ, skin element, and PBSX, as well as the insertion element IS4Bsu1, could be used to identify one of these strains as a Marburg type and the other three strains as natto types. The genome structure and gene heterogeneity were also consistent with the type of indigenous plasmids harbored by the strains.     

The genome of Rhizobium leguminosarum has recognizable core and accessory components

Background

Rhizobium leguminosarum is an α-proteobacterial N₂-fixing symbiont of legumes that has been the subject of more than a thousand publications. Genes for the symbiotic interaction with plants are well studied, but the adaptations that allow survival and growth in the soil environment are poorly understood. We have sequenced the genome of R. leguminosarum biovar viciae strain 3841.

Results

The 7.75 Mb genome comprises a circular chromosome and six circular plasmids, with 61% G+C overall. All three rRNA operons and 52 tRNA genes are on the chromosome; essential protein-encoding genes are largely chromosomal, but most functional classes occur on plasmids as well. Of the 7,263 protein-encoding genes, 2,056 had orthologs in each of three related genomes (Agrobacterium tumefaciens, Sinorhizobium meliloti, and Mesorhizobium loti), and these genes were over-represented in the chromosome and had above average G+C. Most supported the rRNA-based phylogeny, confirming A. tumefaciens to be the closest among these relatives, but 347 genes were incompatible with this phylogeny; these were scattered throughout the genome but were over-represented on the plasmids. An unexpectedly large number of genes were shared by all three rhizobia but were missing from A. tumefaciens.

Conclusion

Overall, the genome can be considered to have two main components: a 'core', which is higher in G+C, is mostly chromosomal, is shared with related organisms, and has a consistent phylogeny; and an 'accessory' component, which is sporadic in distribution, lower in G+C, and located on the plasmids and chromosomal islands. The accessory genome has a different nucleotide composition from the core despite a long history of coexistence.     

Draft Genome Comparison of Representatives of the Three Dominant Genotype Groups of Dairy Bacillus licheniformis Strains

The spore-forming bacterium Bacillus licheniformis is a common contaminant of milk and milk products. Strains of this species isolated from dairy products can be differentiated into three major groups, namely, G, F1, and F2, using random amplification of polymorphic DNA (RAPD) analysis; however, little is known about the genomic differences between these groups and the identity of the fragments that make up their RAPD profiles. In this work we obtained high-quality draft genomes of representative strains from each of the three RAPD groups (designated strain G-1, strain F1-1, and strain F2-1) and compared them to each other and to B. licheniformis ATCC 14580 and Bacillus subtilis 168. Whole-genome comparison and multilocus sequence typing revealed that strain G-1 contains significant sequence variability and belongs to a lineage distinct from the group F strains. Strain G-1 was found to contain genes coding for a type I restriction modification system, urease production, and bacitracin synthesis, as well as the 8-kbp plasmid pFL7, and these genes were not present in strains F1-1 and F2-1. In agreement with this, all isolates of group G, but no group F isolates, were found to possess urease activity and antimicrobial activity against Micrococcus. Identification of RAPD band sequences revealed that differences in the RAPD profiles were due to differences in gene lengths, 3′ ends of predicted primer binding sites, or gene presence or absence. This work provides a greater understanding of the phylogenetic and phenotypic differences observed within the B. licheniformis species.     

Metabolic Engineering of Carotenoid Biosynthesis in Escherichia coli by Ordered Gene Assembly in Bacillus subtilis

We attempted to optimize the production of zeaxanthin in Escherichia coli by reordering five biosynthetic genes in the natural carotenoid cluster of Pantoea ananatis. Newly designed operons for zeaxanthin production were constructed by the ordered gene assembly in Bacillus subtilis (OGAB) method, which can assemble multiple genes in one step using an intrinsic B. subtilis plasmid transformation system. The highest level of production of zeaxanthin in E. coli (820 μg/g [dry weight]) was observed in the transformant with a plasmid in which the gene order corresponds to the order of the zeaxanthin metabolic pathway (crtE-crtB-crtI-crtY-crtZ), among a series of plasmids with circularly permuted gene orders. Although two of five operons using intrinsic zeaxanthin promoters failed to assemble in B. subtilis, the full set of operons was obtained by repressing operon expression during OGAB assembly with a p_R promoter-cI repressor system. This result suggests that repressing the expression of foreign genes in B. subtilis is important for their assembly by the OGAB method. For all tested operons, the abundance of mRNA decreased monotonically with the increasing distance of the gene from the promoter in E. coli, and this may influence the yield of zeaxanthin. Our results suggest that rearrangement of biosynthetic genes in the order of the metabolic pathway by the OGAB method could be a useful approach for metabolic engineering.     

Reconstitution of active 50 S ribosomal subunits from Bacillus licheniformis and Bacillus subtilis.

Active 50 S ribosomal subunits from Bacillus licheniformis and Bacillus subtilis can be reconstituted in vitro from dissociated RNA and proteins. The reconstituted 50 S sub-units are indistinguishable from native 50 S subunits in sedimentation on sucrose gradients and in protein composition. The procedure used is similar to that developed for reconstitution of Bacillus stearothermophilus 50 S subunits, though the optimal conditions are somewhat different. Hybrid ribosomes can be reconstituted with 23 S RNA and proteins from different sources (B. stearothermophilus and B. licheniformis or B. subtilis). The thermal stability of these ribosomes depends on the source of the proteins, and not on the source of 23 S RNA.     

The two authentic methionine aminopeptidase genes are differentially expressed in <Emphasis Type="Italic">Bacillus subtilis</Emphasis>

Background  

Two putative methionine aminopeptidase genes,map(essential) andyflG(non-essential), were identified in the genome sequence ofBacillus subtilis. We investigated whether they can function as methionine aminopeptidases and further explored possible reasons for their essentiality or dispensability inB. subtilis.