The complete genome sequence of the gram-positive bacterium Bacillus subtilis

Strand compositional asymmetry in bacterial and large viral genomes

Several bacterial genomes exhibit preference for G over C on the DNA leading strand extending from the origin of replication to the ter-region in the genomes of Escherichia coli, Mycoplasma genitalium, Bacillus subtilis, and marginally in Haemophilus influenzae, Mycoplasma pneumoniae, and Helicobacter pylori. Strand compositional asymmetry is not observed in the cyanobacterium Synechocystis sp. genome nor in the archaeal genomes of Methanococcus jannaschii, Methanobacterium thermoautotrophicum, and Archaeoglobus fulgidus. A strong strand compositional asymmetry is observed in beta-type but not alpha- or gamma-type human herpesviruses featuring G > C downstream of oriL and C > G upstream of oriL. Dinucleotide relative abundances (i.e., dinucleotide representations normalized by the component nucleotide frequencies) are consonant with respect to the leading and lagging strands. Strand compositional asymmetry may reflect on differences in replication synthesis of the leading versus lagging strand, on differences between template and coding strand associated with transcription-coupled repair mechanisms, on differences in gene density between the two strands, on differences in residue and codon biases in relation to gene function, expression level, or operon organization, or on differences in single or context-dependent base mutational rates. The absence of strand asymmetry in the archaeal genomes may reflect the presence of multiple origins of replication.

Archaeal translation initiation revisited: the initiation factor 2 and eukaryotic initiation factor 2B alpha-beta-delta subunit families

As the amount of available sequence data increases, it becomes apparent that our understanding of translation initiation is far from comprehensive and that prior conclusions concerning the origin of the process are wrong. Contrary to earlier conclusions, key elements of translation initiation originated at the Universal Ancestor stage, for homologous counterparts exist in all three primary taxa. Herein, we explore the evolutionary relationships among the components of bacterial initiation factor 2 (IF-2) and eukaryotic IF-2 (eIF-2)/eIF-2B, i.e., the initiation factors involved in introducing the initiator tRNA into the translation mechanism and performing the first step in the peptide chain elongation cycle. All Archaea appear to posses a fully functional eIF-2 molecule, but they lack the associated GTP recycling function, eIF-2B (a five-subunit molecule). Yet, the Archaea do posses members of the gene family defined by the (related) eIF-2B subunits alpha, beta, and delta, although these are not specifically related to any of the three eukaryotic subunits. Additional members of this family also occur in some (but by no means all) Bacteria and even in some eukaryotes. The functional significance of the other members of this family is unclear and requires experimental resolution. Similarly, the occurrence of bacterial IF-2-like molecules in all Archaea and in some eukaryotes further complicates the picture of translation initiation. Overall, these data lend further support to the suggestion that the rudiments of translation initiation were present at the Universal Ancestor stage.

Identification and characterization of a bacterial chromosome partitioning site

We have identified a DNA site involved in chromosome partitioning in B. subtilis. This site was identified in vivo as the binding site for the chromosome partitioning protein Spo0J, a member of the ParB family of partitioning proteins. Spo0J is a site-specific DNA-binding protein that recognizes a 16 bp sequence found in spo0J. Allowing two mismatches, this sequence occurs ten times in the entire B. subtilis chromosome, all in the origin-proximal approximately 20%. Eight of the ten sequences are bound to Spo0J in vivo. The presence of a site on an otherwise unstable plasmid stabilized the plasmid in a Spo0J-dependent manner, demonstrating that this site, called parS, can function as a partitioning site. This site and Spo0J are conserved in a wide range of bacterial species.

A novel protein kinase that controls carbon catabolite repression in bacteria

HPr(Ser) kinase is the sensor in a multicomponent phosphorelay system that controls catabolite repression, sugar transport and carbon metabolism in gram-positive bacteria. Unlike most other protein kinases, it recognizes the tertiary structure in its target protein, HPr, a phosphocarrier protein of the bacterial phosphotransferase system and a transcriptional cofactor controlling the phenomenon of catabolite repression. We have identified the gene (ptsK) encoding this serine/threonine protein kinase and characterized the purified protein product. Orthologues of PtsK have been identified only in bacteria. These proteins constitute a novel family unrelated to other previously characterized protein phosphorylating enzymes. The Bacillus subtilis kinase is shown to be allosterically activated by metabolites such as fructose 1,6-bisphosphate and inhibited by inorganic phosphate. In contrast to wild-type B. subtilis, the ptsK mutant is insensitive to transcriptional regulation by catabolite repression. The reported results advance our understanding of phosphorylation-dependent carbon control mechanisms in Gram-positive bacteria.

Reconstruction of amino acid biosynthesis pathways from the complete genome sequence

The complete genome sequence of an organism contains information that has not been fully utilized in the current prediction methods of gene functions, which are based on piece-by-piece similarity searches of individual genes. We present here a method that utilizes a higher level information of molecular pathways to reconstruct a complete functional unit from a set of genes. Specifically, a genome-by-genome comparison is first made for identifying enzyme genes and assigning EC numbers, which is followed by the reconstruction of selected portions of the metabolic pathways by use of the reference biochemical knowledge. The completeness of the reconstructed pathway is an indicator of the correctness of the initial gene function assignment. This feature has become possible because of our efforts to computerize the current knowledge of metabolic pathways under the KEGG project. We found that the biosynthesis pathways of all 20 amino acids were completely reconstructed in Escherichia coli, Haemophilus influenzae, and Bacillus subtilis, and probably in Synechocystis and Saccharomyces cerevisiae as well, although it was necessary to assume wider substrate specificity for aspartate aminotransferases.

The manganese-containing ribonucleotide reductase of Corynebacterium ammoniagenes is a class Ib enzyme

Ribonucleotide reductases (RNRs) are key enzymes in living cells that provide the precursors of DNA synthesis. The three characterized classes of RNRs differ by their metal cofactor and their stable organic radical. We have purified to near homogeneity the enzymatically active Mn-containing RNR of Corynebacterium ammoniagenes, previously claimed to represent a fourth RNR class. N-terminal and internal peptide sequence analyses clearly indicate that the C. ammoniagenes RNR is a class Ib enzyme. In parallel, we have cloned a 10-kilobase pair fragment from C. ammoniagenes genomic DNA, using primers specific for the known class Ib RNR. The cloned class Ib locus contains the nrdHIEF genes typical for class Ib RNR operon. The deduced amino acid sequences of the nrdE and nrdF genes matched the peptides from the active enzyme, demonstrating that C. ammoniagenes RNR is composed of R1E and R2F components typical of class Ib. We also show that the Mn-containing RNR has a specificity for the NrdH-redoxin and a response to allosteric effectors that are typical of class Ib RNRs. Electron paramagnetic resonance and atomic absorption analyses confirm the presence of Mn as a cofactor and show, for the first time, insignificant amounts of iron and cobalt found in the other classes of RNR. Our discovery that C. ammoniagenes RNR is a class Ib enzyme and possesses all the highly conserved amino acid side chains that are known to ligate two ferric ions in other class I RNRs evokes new, challenging questions about the control of the metal site specificity in RNR. The cloning of the entire NrdHIEF locus of C. ammoniagenes will facilitate further studies along these lines.

New protein kinase and protein phosphatase families mediate signal transduction in bacterial catabolite repression

Carbon catabolite repression (CCR) is the prototype of a signal transduction mechanism. In enteric bacteria, cAMP was considered to be the second messenger in CCR by playing a role reminiscent of its actions in eukaryotic cells. However, recent results suggest that CCR in Escherichia coli is mediated mainly by an inducer exclusion mechanism. In many Gram-positive bacteria, CCR is triggered by fructose-1,6-bisphosphate, which activates HPr kinase, presumed to be one of the most ancient serine protein kinases. We here report cloning of the Bacillus subtilis hprK and hprP genes and characterization of the encoded HPr kinase and P-Ser-HPr phosphatase. P-Ser-HPr phosphatase forms a new family of phosphatases together with bacterial phosphoglycolate phosphatase, yeast glycerol-3-phosphatase, and 2-deoxyglucose-6-phosphate phosphatase whereas HPr kinase represents a new family of protein kinases on its own. It does not contain the domain structure typical for eukaryotic protein kinases. Although up to now the HPr modifying/demodifying enzymes were thought to exist only in Gram-positive bacteria, a sequence comparison revealed that they also are present in several Gram-negative pathogenic bacteria.

GeneMark.hmm: new solutions for gene finding

The number of completely sequenced bacterial genomes has been growing fast. There are computer methods available for finding genes but yet there is a need for more accurate algorithms. The GeneMark. hmm algorithm presented here was designed to improve the gene prediction quality in terms of finding exact gene boundaries. The idea was to embed the GeneMark models into naturally derived hidden Markov model framework with gene boundaries modeled as transitions between hidden states. We also used the specially derived ribosome binding site pattern to refine predictions of translation initiation codons. The algorithm was evaluated on several test sets including 10 complete bacterial genomes. It was shown that the new algorithm is significantly more accurate than GeneMark in exact gene prediction. Interestingly, the high gene finding accuracy was observed even in the case when Markov models of order zero, one and two were used. We present the analysis of false positive and false negative predictions with the caution that these categories are not precisely defined if the public database annotation is used as a control.

Use of asymmetric cell division and spoIIIE mutants to probe chromosome orientation and organization in Bacillus subtilis

Soon after the onset of sporulation in Bacillus subtilis, asymmetric cell division occurs to generate the differentiating prespore and mother cell types. Formation of the septum close to the cell pole initially bisects the nucleoid destined for the prespore, trapping only about one-third of the DNA in the small compartment. The remaining part of the chromosome is then transported through the septum. spoIIIE mutant cells fail to transfer the DNA and arrest with only partially segregated prespore chromosomes. Previous work has shown that the orientation of the chromosome at the time of septation is not random. Here, we use both physical and genetic methods to characterize the trapped DNA. The results show that the chromosome has a very specific orientation at the time of septation, consistent with the action of a centromere-like sequence near oriC. They also demonstrate that the chromosome is folded, or otherwise organized, in a highly ordered manner.

Rolling-circle plasmids from Bacillus subtilis: complete nucleotide sequences and analyses of genes of pTA1015, pTA1040, pTA1050 and pTA1060, and comparisons with related plasmids from gram-positive bacteria

Most small plasmids of Gram-positive bacteria use the rolling-circle mechanism of replication and several of these have been studied in considerable detail at the DNA level and for the function of their genes. Although most of the common laboratory Bacillus subtilis 168 strains do not contain plasmids, several industrial strains and natural soil isolates do contain rolling-circle replicating (RCR) plasmids. So far, knowledge about these plasmids was mainly limited to: (i) a classification into seven groups, based on size and restriction patterns; and (ii) DNA sequences of the replication region of a limited number of them. To increase the knowledge, also with respect to other functions specified by these plasmids, we have determined the complete DNA sequence of four plasmids, representing different groups, and performed computer-assisted and experimental analyses on the possible function of their genes. The plasmids analyzed are pTA1015 (5.8 kbp), pTA1040 (7.8 kbp), pTA1050 (8.4 kbp), and pTA1060 (8.7 kbp). These plasmids have a structural organization similar to most other known RCR plasmids. They contain highly related replication functions, both for leading and lagging strand synthesis. pTA1015 and pTA1060 contain a mobilization gene enabling their conjugative transfer. Strikingly, in addition to the conserved replication modules, these plasmids contain unique module(s) with genes which are not present on known RCR plasmids of other Gram-positive bacteria. Examples are genes encoding a type I signal peptidase and genes encoding proteins belonging to the family of response regulator aspartate phosphatases. The latter are likely to be involved in the regulation of post-exponential phase processes. The presence of these modules on plasmids may reflect an adaptation to the special conditions to which the host cells were exposed.

Arabidopsis NPH1: a protein kinase with a putative redox-sensing domain

The NPH1 (nonphototropic hypocotyl 1) gene encodes an essential component acting very early in the signal-transduction chain for phototropism. Arabidopsis NPH1 contains a serine-threonine kinase domain and LOV1 and LOV2 repeats that share similarity (36 to 56 percent) with Halobacterium salinarium Bat, Azotobacter vinelandii NIFL, Neurospora crassa White Collar-1, Escherichia coli Aer, and the Eag family of potassium-channel proteins from Drosophila and mammals. Sequence similarity with a known (NIFL) and a suspected (Aer) flavoprotein suggests that NPH1 LOV1 and LOV2 may be flavin-binding domains that regulate kinase activity in response to blue light-induced redox changes.