Cloning and Sequencing of a 36-kb Region of the Bacillus subtilis Genome between the gnt and iol Operons

Within the framework of an international project for the sequencingof the entire Bacillus subtilis genome, a 36-kb chromosome segment,which covers the region between the gnt and iol operons, hasbeen cloned and sequenced. This region (36447 bp) contains 33complete open reading frames (ORFs; genes) including the fourgnt genes and one partial gene. A homology search for the productsof the 33 complete ORFs revealed significant homology to knownproteins in 16 of them such as tetracycline resistance protein(Clostridium perfringens), asparagine synthetase (Arabidopsisthaliana), aldehyde dehydrogenase (Pseudomonas oleovorans),2,5-dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase (P. paucimobilis),heat shock protein HtpG (Escherichia coli), galactose-protonsymporter (E. coli), auxin-induced protein (common tobacco),glucitol operon repressor (E. coli) and methylmalonate-semialdehydedehydrogenase (P. aeruginosa). Unlike the regions we sequencedso far, this region contained two short sequence multiplications:one was a tandem sequence duplication (409 and 410 bp), andthe other a triplication consisting of two highly conserved118-bp tandem sequences preceded by a less conserved similarsequence (129 bp). The reasons for the presence of these sequencemultiplications in the gnt to iol region were deduced.     

Cloning and Sequencing of a 23-kb Region of the Bacillus subtilis Genome between the iol and hut Operons

Within the framework of an international project for the sequencingof the entire Bacillus subtilis genome, a 23-kb chromosomalsegment, which covers the region between the iol and hut operons,has been cloned and sequenced, creating a 99-kb contig fromthe gnt operon to the wapA locus. This region (23351 bp) contains25 complete open reading frames (ORFs; genes) including deoR,dra, nupC and pdp and two partial ones. The region (5140 bp)containing these four genes, being also sequenced by H. H. Saxildet al., was sequenced by subjecting a long polymerase chainreaction product to random sequencing using phage M13mp19. However,we could detect no conflict, between two independently determinedsequences, which could be attributed to our sequencing method.A homology search for the 24 newly identified gene productsrevealed significant homology to known proteins in 14 of them.It was notable that three proteins, encoded by the successivegenes (yxeMNO), exhibited meaningful homology to the E. coliGlnHPQ products constituting a periplasmic ATP-dependent transportsystem for glutamine.     

Complete genome sequence of the industrial bacterium Bacillus licheniformis and comparisons with closely related Bacillus species

Background

Bacillus licheniformis is a Gram-positive, spore-forming soil bacterium that is used in the biotechnology industry to manufacture enzymes, antibiotics, biochemicals and consumer products. This species is closely related to the well studied model organism Bacillus subtilis, and produces an assortment of extracellular enzymes that may contribute to nutrient cycling in nature.

Results

We determined the complete nucleotide sequence of the B. licheniformis ATCC 14580 genome which comprises a circular chromosome of 4,222,336 base-pairs (bp) containing 4,208 predicted protein-coding genes with an average size of 873 bp, seven rRNA operons, and 72 tRNA genes. The B. licheniformis chromosome contains large regions that are colinear with the genomes of B. subtilis and Bacillus halodurans, and approximately 80% of the predicted B. licheniformis coding sequences have B. subtilis orthologs.

Conclusions

Despite the unmistakable organizational similarities between the B. licheniformis and B. subtilis genomes, there are notable differences in the numbers and locations of prophages, transposable elements and a number of extracellular enzymes and secondary metabolic pathway operons that distinguish these species. Differences include a region of more than 80 kilobases (kb) that comprises a cluster of polyketide synthase genes and a second operon of 38 kb encoding plipastatin synthase enzymes that are absent in the B. licheniformis genome. The availability of a completed genome sequence for B. licheniformis should facilitate the design and construction of improved industrial strains and allow for comparative genomics and evolutionary studies within this group of Bacillaceae.     

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.     

Sequence Analysis of the groESL-cotA Region of the Bacillus subtilis Genome, Containing the Restriction/Modification System Genes

We have determined a 35-kb sequence of the groESL-gutR-cotA(45°–52°) region of the Bacillus subtilis genome.In addition to the groESL, gutRB and cotA genes reported previously,we have newly identified 24 ORFs including gutA and fruC genes,encoding glucitol permease and fructokinase, respectively. Theinherent restriction/modification system genes, hsdMR and hsdMM,were mapped between groESL and gutRB, and we have identifiedtwo open reading frames (ORFs) encoding 5-methylcytosine formingDNA methyl transferase and an operon probably encoding a restrictionenzyme complex. The unusual genome structure of few ORFs andlower GC content around the restriction/modification genes stronglysuggests that the region originated from a bacteriophage integratedduring evolution.     

Members of the pogo superfamily of DNA-mediated transposons in the human genome

Bacillus subtilis, likeEscherichia coli, possesses several sets of genes involved in the utilization of-glucosides. InE. coli, all these genes are cryptic, including the genes forming thebgl operon, thus leading to a Bgl^– phenotype. We screened forB. subtilis chromosomal DNA fragments capable of reverting the Bgl⁺ phenotype associated with anE. coli hns mutant to the Bgl^– wild-type phenotype. OneB. subtilis chromosomal fragment having this property was selected. It contained a putative Ribonucleic AntiTerminator binding site (RAT sequence) upstream from thebglP gene. Deletion studies as well as subcloning experiments allowed us to prove that the putativeB. subtilis bglP RAT sequence was responsible for the repression of theE. coli bgl operon. We propose that this repression results from the titration of the BglG antiterminator protein ofE. coli bgl operon by our putativeB. subtilis bglP RAT sequence. Thus, we report evidence for a new cross interaction between heterologous RAT-antiterminator protein pairs.     

Sequence Analysis of the 36-kb Region Between gntZ and trnY Genes of Bacillus Subtilis Genome*

Within the framework of an international Bacillus subtilis genomesequencing project, we have determined a 36-kb sequence coveringthe region between the gntZ and trnY genes. In addition to fivegenes sequenced and characterized previously, 27 putative proteincoding sequences (open reading frame; ORF) were identified.A homology search for the newly identified ORFs revealed thatsix of them had similarities to known proteins. It is notablethat new ORFs belonging to response-regulator aspartate phosphatase(Rap) and its regulator (Phr) families, and response regulatorand sensory kinase families of two-component signal transductionsystems have been identified. Furthermore, we found that some180-bp non-coding sequence, that might be an remnant of an ancientIS element, is preserved in at least five loci of the B. subtilisgenome.     

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.     

Members of the pogo superfamily of DNA-mediated transposons in the human genome

Bacillus subtilis, likeEscherichia coli, possesses several sets of genes involved in the utilization ofβ-glucosides. InE. coli, all these genes are cryptic, including the genes forming thebgl operon, thus leading to a Bgl^? phenotype. We screened forB. subtilis chromosomal DNA fragments capable of reverting the Bgl⁺ phenotype associated with anE. coli hns mutant to the Bgl^? wild-type phenotype. OneB. subtilis chromosomal fragment having this property was selected. It contained a putative Ribonucleic AntiTerminator binding site (RAT sequence) upstream from thebglP gene. Deletion studies as well as subcloning experiments allowed us to prove that the putativeB. subtilis bglP RAT sequence was responsible for the repression of theE. coli bgl operon. We propose that this repression results from the titration of the BglG antiterminator protein ofE. coli bgl operon by our putativeB. subtilis bglP RAT sequence. Thus, we report evidence for a new cross interaction between heterologous RAT-antiterminator protein pairs.     

Sequence analysis and functional studies of a chromosomal region of alkaliphilic Bacillus firmus OF4 encoding an ABC-type transporter with similarity of sequence and Na+ exclusion capacity to the Bacillus subtilis NatAB transporter

A 14.1-kb DNA fragment was cloned from a lambda library containing inserts of DNA from alkaliphilic Bacillus firmus OF4 on the basis of its hybridization to a probe from a previously sequenced alkaliphile homolog of the natA gene from Bacillus subtilis. Sequence analysis of the entire fragment revealed that, as in B. subtilis, the natA gene was part of a putative gene locus encoding an ABC-type transporter. In the alkaliphile, the transporter involved three genes, designated natCAB, that are part of a larger operon of unknown function. This is in contrast to the two-gene natAB operon and to another homolog from B. subtilis, the yhaQP genes. Like natAB, however, the alkaliphile natCAB catalyzes Na⁺ extrusion as assessed in a mutant of Escherichia coli that is deficient in Na⁺ extrusion. The full 14.1-kb fragment of alkaliphile DNA sequenced in this study contained several probable operons as well as likely monocistronic units. Among the 17 predicted ORFs apart from natCAB were acsA, a homolog of a halobacterial gene encoding acetylCoA synthetase; sspA, a homolog of a small acid-soluble spore protein; and malK, an ATP-binding component that was unaccompanied by candidates for other mal transport genes but was able to complement a malK-deficient mutant of E. coli. No strong candidates for genes encoding a secondary Na⁺/H⁺ antiporter were found in the fragment, either from the sequence analysis or from analyses of complementation of E. coli mutants by subclones of the 14.1-kb piece. There were a total of 12 ORFs whose closest and significant homologs were genes from B. subtilis; of these, one-third were in apparently different contexts, as assessed by the sequence of the neighboring genes, than the B. subtilis homologs. Received: August 30, 1998 / Accepted: November 13, 1998     

Fermentation production of keratinase from Bacillus licheniformisPWD-1 and a recombinant B. subtilis FDB-29

Bacillus licheniformis PWD-1, the parent strain, and B. subtilis FDB-29, a recombinant strain. In both strains, keratinase was induced by proteinaceous media, and repressed by carbohydrates. A seed culture of B. licheniformis PWD-1 at early age, 6–10 h, is crucial to keratinase production during fermentation, but B. subtilis FDB-29 is insensitive to the seed culture age. During the batch fermentation by both strains, the pH changed from 7.0 to 8.5 while the keratinase activity and productivity stayed at high levels. Control of pH, therefore, is not necessary. The temperature for maximum keratinase production is 37°C for both strains, though B. licheniformis is thermophilic and grows best at 50°C. Optimal levels of dissolved oxygen are 10% and 20% for B. licheniformis and B. subtilis respectively. A scale-up procedure using constant temperature at 37°C was adopted for B. subtilis. On the other hand, a temperature-shift procedure by which an 8-h fermentation at 50°C for growth followed by a shift to 37°C for enzyme production was used for B. licheniformis to shorten the fermentation time and increase enzyme productivity. Production of keratinase by B. licheniformis increased by ten-fold following this new procedure. After respective optimization of fermentation conditions, keratinase production by B. licheniformis PWD-1 is approximately 40% higher than that by B. subtilis FDB-29. Received 16 July 1998/ Accepted in revised form 07 March 1999     

Sequence Analysis of a 25-kb Segment in the 17{degrees}-19{degrees} Region of the Bacillus subtilis Chromosome Containing ada Locus

As a part of the Bacillus subtilis genome sequencing project,we have determined a 25-kb sequence covering the 17°–19°region. This region contains 26 complete open reading frames(ORFs) including the alkA and adaA/B operon, which encode genesfor adaptive response to DNA alkylation. A homology search forthe newly identified 21 ORFs revealed that 4 of them exhibita significant similarity to known proteins, e.g., methicillin-resistantStaphylococcus aureus (MRSA) protein homolog, proteins involvedin chloramphenicol resistance, glucosamine synthase and an ABCtransporter protein. The remaining 17 ORFs did not show anysignificant sequence similarities to known gene products inthe database.     

Photo‐inactivation of Bacillus endospores: inter‐specific variability of inactivation efficiency

The aims of this work were to (a) evaluate the susceptibility of endospores of Bacillus cereus, B. licheniformis, B. sphaericus and B. subtilis to photodynamic inactivation using a tricationic porphyrin as photosensitizer, (b) assess the efficiency of adsorption of the photosensitizer in endospore material as a determinant of the susceptibility of endospores of different Bacillus species to photo‐inactivation, (c) determine the value of B. cereus as a model organism for studies of antimicrobial photodynamic inactivation of bacterial endospores. The results of irradiation experiments with endospores of four species of Bacillus showed that B. cereus was the only species for which efficient endospore photo‐inactivation (> 3 log reduction) could be achieved. Endospores of B. licheniformis, B. sphaericus and B. subtilis were virtually resistant to photo‐inactivation with tricationic porphyrin. The amount of porphyrin bound to endospore material was not significantly different between species, regardless of the presence of an exosporium or exosporium‐like outer layer. The sensitivity of endospores to photodynamic inactivation with a tricationic porphyrin is highly variable among different species of the genus Bacillus. The presence of an exosporium in endospores of B. cereus and B. sphaericus, or an exosporium‐like glycoprotein layer in endospores of B. subtilis, did not affect the amount of bound photosensitizer and did not explain the inter‐species variability in susceptibility to photodynamic inactivation. The results imply that the use of B. cereus as a more amenable surrogate of the exosporium‐producing B. anthracis must be carefully considered when testing new photosensitizers for their antimicrobial photo‐inactivation properties.     

Detailed identification of desert-originated bacteria carried by Asian dust storms to Japan

Several halotolerant bacteria were isolated from dust allowed to settle passively on saline medium in Higashi-Hiroshima, Japan during Asia dust events in 2005–2006. The primary identification, based on the sequence similarity of the 16S rRNA gene, revealed that these isolates were strains of Bacillus subtilis, B. licheniformis, Staphylococcus epidermidis, Gracillibacillus sp., and Halomonas venusta. A parallel investigation carried out on desert sand collected directly from sand dunes in Dunhuang, Gobi Desert, China resulted in the revivification of seven bacterial strains that were highly identical to the B. subtilis and B. licheniformis strains obtained in Higashi-Hiroshima (99.7 and 100% of 16S rDNA sequence similarity, respectively). A subsequent genetic analysis on the group of B. licheniformis isolates based on the universally house-keeping genes, gyrB and parE, revealed high sequence similarities in both genes among the strains of both locations (99.0–99.4%), which clustered them in a monophyletic line. Phenotype characterized by numerical taxonomy for 150 physiological tests confirmed the close relatedness between strains (similarity coefficient S _SM = 96.0%). The remarkable agreement between phenotype and genotype of the bacterial isolates allows us to conclude that there may have been an aerosolized dispersion of a Gobi Desert B. licheniformis by dust storms to Japan. This study provides evidence of microbial transport by yellow dust events in North-East Asia.     

Cloning and sequencing of a new holin-encoding gene of Bacillus licheniformis

A Bacillus licheniformis DNA fragment which exhibits homology with the upstream region of the cell-wall hydrolase-encoding gene, cwlL, was cloned into Escherichia coli (Ec). Nucleotide sequencing indicated that there are two open reading frames (tentatively designated as xpaG1 and xpaG2) which encode polypeptide of 89 and 88 amino acids (aa) (10044 and 9764 Da, respectively). Ec cells harboring two compatible plasmids (pMWB1 and pHSGKH) containing the Bacillus subtilis cell-wall hydrolase-encoding gene, cwlA, and xpaG1–G2, respectively, exhibited higher extracellular cell-wall hydrolase activity than did cells harboring pMWB1 and a control plasmid, pHSG398. The aa sequence homology of XpaG2 with other polypeptides indicated that xpaG2 is a holin-encoding gene. Moreover, Ec C600 harboring a plasmid containing xpaG1-xpaG2 led to leakage of β-galactosidase into the extracellular fraction.     

Genetic engineering of Bacillus subtilis for the commercial production of riboflavin

Recombinant DNA engineering was combined with mutant selection and fermentation improvement to develop a strain of Bacillus subtilis that produces commercially attractive levels of riboflavin. The B. subtilis riboflavin production strain contains multiple copies of a modified B. subtilis riboflavin biosynthetic operon (rib operon) integrated at two different sites in the B. subtilis chromosome. The modified rib operons are expressed constitutively from strong phage promoters located at the 5′ end and in an internal region of the operon. The engineered strain also contains purine analog-resistant mutations designed to deregulate the purine pathway (GTP is the precursor for riboflavin), and a riboflavin analog-resistant mutation in ribC that deregulates the riboflavin biosynthetic pathway. Received 22 June 1998/ Accepted in revised form 6 November 1998